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Vorwort des Herausgebers 


Der schnelle technische Fortschritt im Turbomaschinenbau, der durch extreme technische For- 
derungen und starken internationalen Wettbewerb geprägt ist, verlangt einen effizienten Aus- 
tausch und die Diskussion von Fachwissen und Erfahrung zwischen Universitäten und indus- 
triellen Partnern. Mit der vorliegenden Reihe haben wir versucht, ein Forum zu schaffen, das 
neben unseren Publikationen in Fachzeitschriften die aktuellen Forschungsergebnisse des Insti- 
tuts für Thermische Strömungsmaschinen am Karlsruher Institut für Technologie (KIT) einem 
möglichst großen Kreis von Fachkollegen aus der Wissenschaft und vor allem auch der Praxis 
zugänglich macht und den Wissenstransfer intensiviert und beschleunigt. 


Flugtriebwerke, stationäre Gasturbinen, Turbolader und Verdichter sind im Verbund mit den 
zugehörigen Anlagen faszinierende Anwendungsbereiche. Es ist nur natürlich, dass die me- 
thodischen Lösungsansätze, die neuen Messtechniken, die Laboranlagen auch zur Lösung von 
Problemstellungen in anderen Gebieten - hier denke ich an Otto- und Dieselmotoren, elektrische 
Antriebe und zahlreiche weitere Anwendungen - genutzt werden. Die effiziente, umweltfreund- 
liche und zuverlässige Umsetzung von Energie führt zu Fragen der ein- und mehrphasigen Strö- 
mung, der Verbrennung und der Schadstoffbildung, des Wärmeübergangs sowie des Verhaltens 
metallischer und keramischer Materialien und Verbundwerkstoffe. Sie stehen im Mittelpunkt 
ausgedehnter theoretischer und experimenteller Arbeiten, die im Rahmen nationaler und inter- 
nationaler Forschungsprogramme in Kooperation mit Partnern aus Industrie, Universitäten und 
anderen Forschungseinrichtungen durchgeführt werden. 


Es sollte nicht unerwähnt bleiben, dass alle Arbeiten durch enge Kooperation innerhalb des 
Instituts geprägt sind. Nicht ohne Grund ist der Beitrag der Werkstätten, der Technik-, der 
Rechner- und Verwaltungsabteilungen besonders hervorzuheben. Diplomanden und Hilfsassis- 
tenten tragen mit ihren Ideen Wesentliches bei, und natürlich ist es der stets freundschaftlich 
fordernde wissenschaftliche Austausch zwischen den Forschergruppen des Instituts, der zur 
gleichbleibend hohen Qualität der Arbeiten entscheidend beiträgt. Dabei sind wir für die Unter- 
stützung unserer Förderer außerordentlich dankbar. 


Aufgrund steigender Leistungsdichten und zur Erreichung höherer Prozesswirkungsgrade 
nimmt die thermische Belastung insbesondere der ersten Turbinenstufen moderner Flugtrieb- 
werke stetig zu. Turbineneintrittstemperaturen, die weit über der Einsatztemperatur modernster 
Superlegierungen liegen, sind nur durch die Verwendung fortschrittlicher Kühlmethoden zu 
beherrschen. Voraussetzung für eine weitere Verbesserung der Effektivität und Effizienz der 
Kühlung ist die genaue Kenntnis der externen und internen Prozesse der Wärmeübertragung. 
Im vorliegenden Band der Schriftenreihe verfolgt der Autor den Ansatz, eine geometrisch 
ähnliche aber gegenüber der realen Anwendung vergrößerte Laufschaufel unter Einhaltung 
aller relevanten aerodynamischen und thermischen Ähnlichkeitskenngrößen auf einem sta- 
tionären Heißgasversuchsstand experimentell zu untersuchen. Im Fokus der Untersuchungen 
stand die lokal und thermisch hoch aufgelöste Ermittlung der Gesamtkühleffektivität der 
Schaufel unter wohldefinierten Randbedingungen. Es sollte dabei insbesondere ein Vergleich 
unterschiedlicher Konfigurationen der internen Konvektionskühlung im thermisch besonders 
hoch belasteten Bereich der Vorderkante der Schaufel durchgeführt werden. 


Eine spezielle Herausforderung ergab sich fiir die eingesetzte Infrarotmesstechnik zur Bestim- 
mung der Oberflächentemperatur. Diese Messtechnik wird seit vielen Jahren am Institut fiir 
Thermische Strömungsmaschinen erfolgreich zur Analyse von Kühlsystemen eingesetzt. Sie 
wurde kontinuierlich in Ihrer Genauigkeit verbessert und bezüglich ihrer Einsatzmöglichkeiten 
erweitert. Im Gegensatz zu generischen Untersuchungen an abstrahierten Konfigurationen 
ist die Zugänglichkeit im Falle der untersuchten Laufschaufel stark eingeschränkt. Die von 
der Schaufel ausgehende Wärmestrahlung muss unter einem weiten Winkelbereich auf stark 
gekrümmten Oberflächen und für einen großen geometrischen Tiefenbereich mit der Infra- 
rotthermokamera gemessen werden. Dabei ergeben sich große Variationen der Emissivität der 
Oberfläche. Zudem ist die Schärfentiefe der Kameraoptik nicht ausreichend, um alle inter- 
essierenden Bereiche der Schaufeloberfläche in einer einzigen Aufnahme scharf abzubilden. 
Darüber hinaus ist eine eindeutige Zuordnung zwischen der zweidimensionalen Abbildung 
auf dem Kamerasensor zur komplex dreidimensional geformten Oberfläche a priori nicht 
gegeben. Mit Hilfe eines neuartigen genetischen Algorithmus gelingt es dem Autor, die genaue 
Kameraposition exakt zu identifizieren und auf praktisch beliebig geformten Oberflächen deren 
Temperatur lokal optisch und thermisch hoch aufgelöst zu messen können. Auf Basis der Er- 
gebnisse der thermischen Untersuchungen gibt der Autor wichtige Anregungen für die weitere 
Verbesserung der kombinierten Konvektions- und Filmkühlung zukünftiger hochbelasteter 
Gasturbinenschaufeln. 


Karlsruhe, im März 2019 Hans-Jörg Bauer 


Vorwort des Autors 


Diese Dissertation entstand während meiner Tätigkeit am Institut für Thermische Strömungs- 
maschinen (ITS) am Karlsruher Institut für Technologie (KIT). Das Forschungsprojekt zur Un- 
tersuchung von neuartigen Kühlkonzepten wurde in enger Zusammenarbeit mit Rolls-Royce 
Deutschland definiert und durchgeführt. Die Arbeiten wurden im Rahmen des LuFo IV/V För- 
derprogrammes des BMWi gefördert. 


Großer Dank gilt Herrn Professor Hans-Jörg Bauer, zum einen für die Übernahme des Hauptre- 
ferats, zum anderen aber vor allem für die Möglichkeiten, die mir am Institut für Thermische 
Strömungsmaschinen geboten wurden. Das personelle Umfeld, die Infrastruktur, das entgegen- 
gebrachte Vertrauen sowie die Freiheiten während der Arbeit ermöglichen erst die Durchfüh- 
rung von wissenschaftlichen Arbeiten. 

Weiterhin bedanken möchte ich mich auch bei Herrn Professor Heinz-Peter Schiffer, Fachge- 
bietsleiter Gasturbinen, Luft- und Raumfahrtantriebe an der Technischen Universität Darmstadt 
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Besonders danke ich auch Herrn Dr. Achmed Schulz, dem ehemaligen Leiter der Forschungs- 
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rimenteller Natur haben einen großen Beitrag zum Erfolg dieser Dissertation geleistet. 

Ich möchte auch dem ganzen Kollegium danken, das mich während der Zeit am Institut für 
Thermische Strömungsmaschinen begleitet hat. Die sehr gute Zusammenarbeit, die diversen 
Besprechungen und gemeinsamen Unternehmungen tragen sehr zum angenehmen Arbeitskli- 
ma am Institut bei. So ließen sich auch die (unvermeidbaren) Zeiten, in denen wenig geklappt 
hatte wie es sollte, angenehm überstehen. 
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... Johannes Kneer, der die Grundlagen für die 3D-Transformation in Verbindung mit der Ther- 
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Verwendung von python nahegelegt hat und für die darauffolgende andauernde Beratung im 
Institutsalltag. 
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... Felix von Plehwe, meinen langjährigen Bürokollegen. Die ein oder andere Idee, die in den 
zahllosen Gesprächen aufgekommen ist und diskutiert wurde, ist auch in diese Dissertation 
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Abschließend möchte ich meiner Lebensgefährtin Karin und meiner Familie danken. Karin für 
die Unterstützung in den vergangenen Jahren im Privaten und Beruflichen und an dieser Stelle 
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1 Einleitung 


Die zivile Luftfahrt ist ein stetig wachsender Sektor. Trotz der schon über 100 Millionen Passa- 
giere, die allein in Deutschland pro Jahr in ein Flugzeug steigen, wird vom Institut für Flugha- 
fenwesen und Luftverkehr (2015) weiterhin ein globales Wachstum von 5 % pro Jahr prognos- 
tiziert. Neue Studien einer Arbeitsgruppe der Air Transport Action Group ATAG (2013a) gehen 
außerdem davon aus, dass nur 10 % der heutigen Weltbevölkerung für 80 % der globalen Pas- 
sagierkilometer verantwortlich sind. Mit der weiteren Industrialisierung von Schwellenländern 
wird die Nachfrage an Flugreisen weiter steigen. 

Diese positive Wachstumsprognose führt zu einer großen Herausforderung. So ist die zivile 
Luftfahrt derzeit schon für ca. 2% des weltweiten AusstoBes von CO» verantwortlich. Au- 
ßerdem kann die Bildung von Stickoxiden auch bei modernen Triebwerken nicht vermieden 
werden. Die Emission dieser beiden Stoffe wird momentan und mit hoher Wahrscheinlichkeit 
auch in Zukunft stark reguliert. Die Verwendung von elektrischen Antrieben, wie sie in anderen 
Branchen verfolgt wird, ist für die Luftfahrt nach dem aktuellen Stand der Technik und auch 
in naher Zukunft noch keine Option. Selbst theoretisch erreichbare Energiedichten von elek- 
trochemischen Speichern liegen noch Größenordnungen unter der Energiedichte von Kerosin 
(Nitta et al., 2015). Das aus dem Energiebedarf eines Langstreckenfluges resultierende Gewicht 
des Speichers läge damit weit über technisch umsetzbaren Höchstabfluggewichten. Weitere kri- 
tische Aspekte bei der Verwendung fossiler Brennstoffe wie der Massenverlust während des 
Fluges und die damit erst mögliche Landung großer Flugzeuge können mit elektrochemischen 
Speichern nicht umgesetzt werden. 

Ein klima- und gesetzkonformes Wachstum kann demnach nur durch ein verbessertes, effizi- 
enteres Komplettsystem Luftfahrt erreicht werden. Die Ziele dieser Verbesserung werden von 
führenden Herstellern in Verbund mit der ATAG konkret formuliert und in verschiedenen Be- 
richten festgehalten. So wird bis 2020 eine Kraftstoffeinsparung von 50 % und eine Reduktion 
der Stickoxidemissionen um 80 % angestrebt (Busquin et al., 2001). Dadurch soll das weitere 
Wachstum der Branche CO -neutral sein. 

Bis 2050 soll der Netto-Ausstoß von Kohlenstoffdioxid, trotz des prognostizierten Wachstums, 
den Wert von 2005 erreicht haben (ATAG, 2013b). Dafür werden von europäischen Gruppen wie 
Darecki et al. (2011) konkrete Anforderungen formuliert: Eine Reduktion des CO2-AusstoBes 
pro Passagierkilometer um 75 %, eine Reduktion der Stickoxidemissionen um 90 % muss er- 
reicht werden. 


Zum Erreichen dieser gesetzten Ziele müssen alle Teilsysteme optimiert werden. Einen ho- 
hen Anteil von etwa einem Viertel sehen führende Hersteller wie zum Beispiel Rolls-Royce 
Deutschland (2017) bei den Strahltriebwerken. Maßgeblich für die effiziente Nutzung der Pri- 
märenergie und damit Grundlage für einen niedrigen Ausstoß von Kohlenstoffdioxid ist der 
Gesamtwirkungsgrad des Triebwerkes, der den zur Verfügung gestellten Schub in Relation zur 
zugeführten Kraftstoffmenge stellt. Dieser lässt sich weiter aufteilen in das Produkt des thermi- 
schen Wirkungsgrades des Kerntriebwerkes und des Vortriebswirkungsgrades. Cazalens et al. 
(2014) zeigen die angestrebten Ziele der beiden Wirkungsgrade und die damit einhergehenden 
technischen Herausforderungen auf. 
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Abbildung 1.1: Entwicklung der Turbineneintrittstemperatur von Rolls-Royce Triebwerken in 
Abhängigkeit des Zeitpunkts der Inbetriebnahme (EIS) in Anlehnung an The 
Jet Engine, Rolls-Royce PLC (2005) und Rolls-Royce Deutschland (2017) 


Die Steigerung des Vortriebswirkungsgrades wird vor allem durch eine aerodynamische Opti- 
mierung von Fan und Bypass erreicht. Hier werden größere Fandurchmesser und höhere Neben- 
stromverhältnisse angestrebt. Als Kerntechnologie gilt ein Untersetzungsgetriebe zwischen Fan 
und Niederdruckturbine. Das Kerntriebwerk hingegen soll möglichst kompakt bleiben, muss 
aber den erhöhten Leistungsbedarf des Fans decken können. Gefordert ist eine hohe spezifi- 
sche Arbeit bei hohem thermischen Wirkungsgrad. Eine Analyse des Joule-Brayton-Prozesses 
zeigt, dass dies nur durch eine Steigerung des Prozessdruckverhältnisses II bei gleichzeitig stei- 
gender Prozessmaximaltemperatur, der Turbineneintrittstemperatur 73, möglich ist. So werden 
bei heutigen Strahltriebwerken Druckverhältnisse von II ~ 50 und Turbineneintrittstempera- 
turen von 73 > 1.900K erreicht. Zukünftige Triebwerksgenerationen werden schon jetzt mit 
Druckverhältnissen von II ~ 70 und damit auch weiter steigenden Turbineneintrittstemperatu- 
ren geplant. 

Die Entwicklung der Turbineneintrittstemperatur von Strahltriebwerken des Herstellers Rolls- 
Royce ist in Abbildung 1.1 gezeigt. Seit der Einführung von gekühlten Turbinenkomponenten 
ist ein stetiger Anstieg der Turbineneintrittstemperatur festzustellen. Die Gastemperaturen lie- 
gen dabei inzwischen mehrere hundert Kelvin über der Einsatztemperatur Tg hochentwickelter 
Nickelbasislegierungen. Dies ist nur durch extensive Kühlung der von Heißgas umströmten 
Bauteile möglich. Die Kühlluft wird dem Prozess in den hinteren Verdichterstufen entnommen, 
in der Turbine zur Kühlung verwendet und dabei dem Hauptmassenstrom wieder zugeführt. 
Hierbei ergeben sich mehrere Herausforderungen: Die Kühlluft erreicht durch die Verdichtung 
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eine erhöhte Temperatur (über 900 K bei II = 50), der entnommenen Luft wird keine Wärme 
zugefiihrt und bei der graduellen Wiedereinmischung in der Turbine kann nicht die volle Ver- 
dichterarbeit zuriickgewonnen werden. 

So müssen die Bauteile unter Verwendung von möglichst wenig Kühlluft auf Temperaturen na- 
he der Temperatur des Kühlmediums selbst gekühlt werden. Diese Kühlung muss sowohl die 
Geometrie der gekühlten Bauteile selbst als auch das externe Strömungsfeld berücksichtigen. 
Die hoch komplexen Kühlsysteme müssen exakt ausgelegt werden: Eine Abweichung der Bau- 
teiltemperatur von ~ 20K (~ 2%) vom Auslegungspunkt kann nach Schätzungen von Snowsill 
et al. (2006) eine Lebensdauerreduktion von bis zu 50 % bewirken. 

Es existieren weitreichende experimentelle und numerische Untersuchungen und abgeleitete 
Korrelationen zu einer breiten Spanne verschiedener Kühlverfahren. Diese werden üblicher- 
weise während einer Vorauslegung auf ihre Eignung überprüft, auf die konkreten Geometrien 
übertragen und dann durch umfangreiche Berechnungen mithilfe von CFD und Festkörperrech- 
nungen überprüft. Aufgrund der Komplexität der Aerothermodynamik ist diese Auslegung oft 
ungenau. Daher müssen vor Inbetriebnahme neuer Systeme Thermalfarbentests durchgeführt 
werden. Dafür werden die Prüflinge in sogenannten „Demonstrator Engines“, also Technolo- 
gieträgern, die für Testreihen zur Verwendung kommen, eingebaut. Diese Tests sind zeit- und 
kostenintensiv. 


Im Rahmen dieser Arbeit wird ein Prüfstand konzipiert und aufgebaut, der erstmals die Un- 
tersuchung moderner Kühlkonzepte für Rotorschaufeln unter realen Bedingungen ermöglicht. 
Dies bedeutet, dass reale Schaufel- wie Passagengeometrien verwendet werden und dass sowohl 
für das Heißgas als auch für die Kühlluft die wichtigen dimensionslosen Kennzahlen so einge- 
stellt werden können, wie sie aus den Betriebserfahrungen für Flugtriebwerke bekannt sind. Zur 
Untersuchung kommen verschiedene Variationen von internen Geometrien bei unverändertem 
externen Schaufelprofil sowie unveränderter Filmkühlung. 

Diese Bedingungen bedeuten eine große Herausforderung für die angewendete Messtechnik. 
Oberflächentemperaturen müssen auf gekriimmten, stark gekühlten Oberflächen lokal hoch auf- 
gelöst bestimmt werden. Hierfür wird Infrarotthermographie eingesetzt, die durch neue Verfah- 
ren zur Auswertung verbessert wird. 

Im folgenden Kapitel wird ein Überblick über den aktuellen Stand der Forschung gegeben. 
Zur Definition der Ziele dieser Arbeit sind vor allem die bisherigen Errungenschaften aus den 
Forschungsgebieten der internen Kühlmethoden und der angewandten Infrarotthermographie 
grundlegend wichtig, die dort diskutiert werden. Im darauf folgenden Kapitel werden der neue 
Prüfstand sowie aerodynamische Untersuchungen präsentiert, die die erreichte Ähnlichkeit zwi- 
schen Prüfstand und Maschine zeigen. Es folgen zwei Kapitel, die sich im Detail mit den neuen 
Verfahren zur Temperaturmessung mittels Infrarotthermographie bei komplexen Randbedin- 
gungen beschäftigen. Im Anschluss werden die durch die neuen Verfahren gewonnenen experi- 
mentellen Ergebnisse diskutiert. Ein Fazit und Ausblick schließen die Arbeit ab. 


2 Aerothermische Analyse der Stromungen im Laufgitter 


Die thermische Belastung von Rotorschaufeln ist direkt abhängig vom Strömungsfeld des 
Heißgases und dem Strömungsfeld der Kühlluft im Relativsystem der Schaufel. Im Ver- 
gleich zum davor liegenden Stator werden diese durch die zwei Scheinkraftdichten der 
Zentrifugalbeschleunigung Teeni sowie der Coriolisbeschleunigung Te 


Feen =p(@x (@xF')) sowie foor = —2p (0 x V’) (2.1) 


beeinflusst. Hier beschreibt © die Winkelgeschwindigkeit, 7’ und V’ Radius und Geschwindig- 
keit im Relativsystem. 

Die Außenströmung zeigt ähnliche Primär- und Sekundärströmungsstrukturen wie die Stator- 
strömung. Für Unterschiede ist vor allem die zusätzliche Zentrifugalkraft maßgeblich, die Co- 
rioliskraft ist gering: Der Axialanteil wax des relativen Stromungsvektors w ist kollinear mit der 
vektoriellen Winkelgeschwindigkeit ©, die aus dem Umfangsanteil wy resultierende Coriolis- 
kraft ist wegen der linearen Abhängigkeit in © wesentlich kleiner als die kollineare Zentrifugal- 
kraft. Der radiale Anteil der Relativgeschwindigkeit w, ist klein. Die resultierende Außenströ- 
mung wird in Abschnitt 2.1 beschrieben. 

Die Innenströmung der Schaufelkühlung hingegen wird sowohl von der Zentrifugal- als auch 
der Corioliskraft beeinflusst, da die Hauptströmung meist radial nach außen bzw. innen gerichtet 
ist, also ||| = Vintr: Die Zentrifugalkraft führt zu einer Pumpwirkung nach außen, die Co- 
rioliskraft zur Beschleunigung in Umfangsrichtung und damit zur Ausbildung von Ausgleichs- 
strömungen mit Wirbelform in der u,ax-Ebene. Die komplexen Formen der Kühlluftströmung 
sowie der Einfluss der Rotation werden in Abschnitt 2.2 und Abschnitt 2.3 beschrieben. 


2.1 Aerodynamik der Außenströmung im rotierenden Bezugssystem 


Die Umströmung von Turbinenschaufeln wird meist an ebenen Kaskaden oder zumindest in 
stationären Bezugssystemen untersucht. Nach Darlegung der Sekundärströmungsstrukturen im 
nicht rotierenden System werden darauf aufbauend die Unterschiede zum rotierenden Bezugs- 
system aufgezeigt. In Abbildung 2.1 ist das erwartete Sekundärströmungsfeld im Relativsystem 
des Laufgitters qualitativ dargestellt. 

Frühe Untersuchungen an ebenen Kaskaden werden von Sieverding (1985) in einer Übersicht 
zusammengefasst. Ausgehend von den theoretischen Untersuchungen von Hawthorne (1955) 
beschreibt er die entstehenden Strömungsstrukturen durch Interaktion der Grenzschichten der 
Endwände mit der Schaufel. Sieverding identifiziert ein Aufrollen der Grenzschicht der jewei- 
ligen Endwand an der Schaufelvorderkante und einen sich hieraus bildenden Hufeisenwirbel 
Vhew, der dann sowohl entlang der Druckseite als auch an der Saugseite stromab transportiert 
wird. Der druckseitige Ast wird in der Passage durch den Druckgradienten zur Saugseite um- 
gelenkt, nimmt dabei Masse aus der Grenzschicht auf und vermischt sich mit dem saugseitigen 
Ast zum größeren Passagenwirbel Vpw. Die Größe dieses Wirbels wird durch die Eindringtiefe 
(penetration length, Sharma et al., 1987), also die vom Wirbel überdeckte Höhe an der Schau- 
felhinterkante, beschrieben. Des Weiteren wird ein Ablösegebiet (laminar separation bubble) 
im Bereich großer Beschleunigungen auf der Saugseite gezeigt, dessen Intensität und Größe 
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Staupunkt 


Abbildung 2.1: Sekundärströmungsfeld im rotierenden Bezugssystem, qualitativ. Wirbellinien 
auf der Schaufeloberfläche als gestrichelte Linien. Wirbelstrukturen V. Stau- 
punkte von drei Schaufeln markiert, eine Schaufel dargestellt. Koordinatensys- 
tem der Maschine (Prüfstand in eckigen Klammern). 


mit dem Grenzschichtzustand und der Freistromturbulenz variiert. Gregory-Smith et al. (1992) 
beschreiben vergleichbare Strukturen an explizit als Rotorschaufeln gekennzeichneten Geome- 
trien. Die Experimente werden aber auch ohne Rotation durchgeführt. Wang et al. (1997) be- 
stätigen diese Strukturen experimentell und zeigen zusätzlich kleinere Eckenwirbel zwischen 
Endwand und Schaufel sowie einen durch den Passagenwirbel induzierten Wandwirbel (in der 
Abbildung nicht dargestellt). 

Untersuchungen im rotierenden System zeigen weitestgehend ähnliche Strukturen, die wichti- 
gen Unterschiede werden im Folgenden erläutert. Durch die beschriebenen Sekundärströmun- 
gen im vor dem Rotor liegenden Stator ergibt sich eine w-förmige Verteilung der Axialge- 
schwindigkeit in radialer Richtung. Joslyn et al. (1992) identifizieren die Passagenwirbel und 
deren Einfluss auf die Umlenkung an der Statorhinterkante als Ursache. Nach Joslyn et al. 
(1992) ist der Effekt an der Gehäuseendwand am stärksten. Dies wurde in neueren Messungen 
von Gallus et al. (1995) und Ristic et al. (1999) widerlegt, hier wird die maximale Axialge- 
schwindigkeit nahe der nabenseitigen Endwand bestimmt. Diese Verteilung der Axialgeschwin- 
digkeit führt zu einer Verstärkung der nabennahen Sekundärströmungen. 
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Gallus et al. (1995) beschreiben außerdem den Transport des nabenseitigen Passagenwirbels in 
Richtung der Kanalmitte durch radiale Beschleunigung im rotierenden System. Dies führt zu 
einer höheren Eindringtiefe des nabenseitigen Passagenwirbels im Vergleich zum gehäuseseiti- 
gen Passagenwirbel. Ristic et al. (1999) bestätigen diese Ergebnisse durch LDA-Messungen im 
rotierenden System. Auch sie beschreiben das nabenseitige Wirbelsystem relativ zum nicht ro- 
tierenden System vergrößert und hoch stabil, während das gehäuseseitige Wirbelsystem kleiner 
und fragmentierter ausfällt. 

Auf Basis der von Wang et al. (1997) gezeigten graphischen Darstellung und der erläuterten 
Phänomene wird das erwartete Sekundärströmungsfeld in der Rotorpassage in Abbildung 2.1 
gezeigt. Dort sind die wichtigsten Wirbelsysteme dargestellt: Der Hufeisenwirbel mit beiden 
Ästen, wobei der saugseitige Ast vom Passagenwirbel überdeckt wird (Änderung der Drehrich- 
tung) und das analoge Wirbelsystem an der gehäuseseitigen Endwand, das erst weiter stromab 
in die Kanalmitte wandert. Dadurch kann sich die laminar separation bubble stabilisieren. Auf 
der Saugseite am Schaufelfuß können so zwei Wirbelstaulinien erwartet werden, während an 
der Schaufelspitze, je nach genauen Strömungsbedingungen, keine explizite Trennung der zwei 
Wirbel möglich ist. 


2.2 Kühlung von hoch belasteten Rotorschaufeln 


Bauteile im Turbinenbereich von Gasturbinen sind extrem hohen thermischen Belastungen un- 
terworfen. Vor allem die Endwände und Schaufeln sind dem Heißgas, das bereits heutzutage 
Temperaturen mehrere hundert Kelvin über den Einsatztemperaturen der verwendeten Werk- 
stoffe erreicht, direkt ausgesetzt. Eine besondere Stellung nehmen die Schaufeln der ersten 
Laufreihe ein. Da dem Prozessmedium bis zum Eintritt in den Rotor noch keine Arbeit entzogen 
wird, ist die für die Kühlung maßgebliche relative Recovery-Temperatur Trec,1,reı des Heißgases 
hoch. Diese ergibt sich unter Verwendung des Recovery-Faktors r, der Geschwindigkeit c und 
der spezifischen Wärmekapazität cp (bzw. Machzahl Ma sowie Isentropenexponent k) zu 


2 1 k—1 2 g 
c +r*5-Ma . vVPr laminare GS 
Te=T+r— =R-—— — mi r=f, (2.2) 
2Cp 14-7. Ma Pr turbulente GS 


und ist in Abbildung 2.2 als die für den Wärmeübergang maßgebliche Temperatur dargestellt. 
Zusätzlich zur hohen thermischen Belastung wirken, bedingt durch die hohen Drehzahlen der 
Hochdruckturbine, starke Zentripetalkräfte. Diese führen zu Radialspannungen insbesondere im 
Schaufelfuß und am Übergang zwischen Plattform und Schaufel. Diese kombinierte Belastung 
begünstigt plastische Schädigungsmechanismen wie Kriechen. 
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Abbildung 2.2: Temperaturverläufe an der gekühlten Wand, (a) konvektiv gekühlt, (b) konvek- 
tiv mit keramischer Dämmschicht, (c) konvektiv mit keramischer Dämmschicht 
und Filmkühlung 


Trotz der Verwendung modernster hochwarmfester Materialien auf Nickelbasis! sowie komple- 
xer Fertigungsmethoden wie dem Einkristallguss, kann die geforderte Lebensdauer bei gegebe- 
nen Temperaturen nicht erreicht werden. Die Rotorschaufeln müssen daher intensiv gekühlt 
werden. Hierbei haben sich drei Verfahren etabliert, die meist kombiniert eingesetzt werden. 


2.2.1 Konvektive Kühlung 


Bei der konvektiven Kühlung, dargestellt in Abbildung 2.2(a), wird die zu kühlende Oberfläche 
auf der dem Heißgas abgewandten Seite von Kühlluft überströmt. Der Temperaturverlauf ergibt 
sich aus den zwei Grenzschichten mit den Wärmeübergangskoeffizienten hy und hx sowie 
dem linearen Verlauf im Festkörper unter Annahme Fourierscher Wärmeleitung und konstanter 
Wärmeleitfähigkeit As. Die für die jeweiligen Wärmeübergänge an Innen- und Außenseite 
maßgeblichen Temperaturen sind die Recovery-Temperaturen der Fluide nach Gleichung 2.2. 
Für die meisten Anwendungen kann aufgrund der geringen Geschwindigkeiten auf der Seite 
der Kühlluft Trec.k = Tk ~ Tk angenommen werden. 

Üblicherweise ist die heißgasseitige Nußeltzahl Nuy durch die externe Strömung vorgegeben 
und aufgrund hoher Reynolds- und Machzahlen um bis zu eine Größenordnung höher als die 
interne Nußeltzahl Nug. Dies führt dazu, dass auf der Kühlluftseite sowohl die Nußeltzahl 
als auch die Fläche vergrößert werden müssen, um eine effektive Kühlung zu erreichen. Die 
verwendeten Geometrien sind vielfältig und basieren meist auf einer partiellen Versperrung des 
Kühlluftkanals. So wird die Oberfläche vergrößert und Turbulenz erzeugt. Mögliche Ausfüh- 
rungen sind Rippen (ribs), orthogonal oder angestellt, kreisförmige Dellen oder Erhöhungen 


! Einen Einsatz bei noch höheren Temperaturen erlauben nur keramische Werkstoffe. Diese sind aufgrund der 
Sprödigkeit nicht als Material für Schaufeln geeignet. Neue Entwicklungen von sogenannten Ceramic Matrix 
Composite-Werkstoffen, also Keramik-Keramik-Verbunde, sind zwar Gegenstand der Forschung, ein Einsatz 
für Rotorschaufeln ist noch nicht absehbar, der amerikanische Hersteller GE baut jedoch inzwischen Prototypen 
(Angrand, 2018). 
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Abbildung 2.3: Geometrische Anderungen der Wand und entstehende Rezirkulationsgebiete. 
(i) gerade Rippe, (ii) angestellte Rippe, (iii) Dellen und Erhohungen, (iv) Zylin- 
der 


(dimple and protrusion), zylindrische Einbauten (cylinders oder pin fins) sowie Kombinationen 
derselben. Diese geometrischen Ausführungen sind beispielhaft in Abbildung 2.3 dargestellt. 
Im Rahmen dieser Arbeit werden neben den genannten konventionellen Geometrien Drall- 
konzepte für die Vorderkante untersucht. Hier wird der Vorderkantenkanal der Schaufel ohne 
interne Einbauten ausgeführt. Ein im Schaufelfuß implementierter Drallerzeuger erzeugt eine 
Wirbelströmung im Kanal, die sich positiv auf den Wärmeübergang auswirken kann. Ein 
detaillierter Überblick über die zugrunde liegenden physikalischen Effekte und die zum Einsatz 
kommenden Verfahren wird in Abschnitt 2.3 gegeben. 

Als weitere Methode zur Erhöhung des internen Wärmeübergangs, speziell im Bereich 
der Schaufelvorderkante, kann Prallkühlung eingesetzt werden. Hierbei wird eine gelochte 
Zwischenwand? im Kühlluftkanal eingebaut. Uber diese erste Zwischenwand wird ein hoher 
Druckabfall realisiert, der zur Ausbildung von Kiihlluftstrahlen durch die diskreten Bohrungen 
führt. Diese Strahlen treffen mit hohem Impuls auf die Innenseite der zu kühlenden Schau- 
felwand. Die sich einstellende Verteilung des Wärmeübergangskoeffizienten ist komplex. 
Typische Merkmale sind lokal sehr hohe Wärmeübergangskoeffizienten an den Staupunkten 
der Kühlluftstrahlen und vergleichsweise niedrige Wärmeübergangskoeffizienten dazwischen. 
Durch Verwendung von Prallkühlung mit enger Staffelung lässt sich eine hohe Kühlwirkung 
erreichen. Nachteilig ist ein sehr hoher notwendiger Eintrittsdruck der Kühlluft sowie die 
lokale Inhomogenität der Verteilung des Wärmeübergangskoeffizienten, die zu hohen Tem- 


? Das Verfahren wurde ursprünglich in den mechanisch geringer belasteten Schaufeln des Leitgitters verwendet. 
Hier wird die Zwischenwand meist durch Einsetzen von Blechen nach dem Gussverfahren realisiert. Durch die 
verbesserten Feingussverfahren in jüngster Vergangenheit kann diese Zwischenwand auch direkt im Gussver- 
fahren implementiert werden. Durch die so gesteigerte strukturelle Integrität kann das Kühlverfahren auch in 
Rotorschaufeln verwendet werden. 
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peraturgradienten im Material führen kann. Aufgrund dieser Komplexität, des sehr hohen 
Druckbedarfs und des Aufbaus, der sich deutlich von dem der in den späteren Versuchen 
eingesetzten Geometrien unterscheidet, wird im Folgenden nicht detaillierter auf dieses 
Kühlkonzept eingegangen. 


2.2.2 Filmkühlung 


Bei Verwendung von Filmkühlung wird ein Teil der Kühlluft durch diskrete Bohrungen an die 
Heißgasseite der zu kühlenden Wand geführt. Das sich einstellende Temperaturprofil ist in Ab- 
bildung 2.2(c) beispielhaft für eine zylindrische angestellte Bohrung gezeigt. Durch den Kühl- 
film wird sowohl der treibende Temperaturgradient als auch der Wärmeübergangskoeffizient 
beeinflusst. Nach Einführung der adiabaten Wandtemperatur Taw, dem Wärmeübergangskoef- 
fizienten an der filmgekühlten Wand h; und eines Ersatzwärmeübergangskoeffizienten h lassen 
sich zwei Wandwärmeströme definieren 


dw = hr (Taw = Ty) (2.3a) 
qw = h (Trec,H = Tw) (2.3b) 


wobei h dann auch abhängig von der unbekannten, lokalen Temperatur des Kühlfilms ist. Diese 
zwei Definitionen lassen sich nach Choe et al. (1974) unter Annahme eines stationären Pro- 
zesses vergleichen. Zur experimentellen Bestimmung der neuen Unbekannten Tw, Taw sowie hr 
müssen mindestens zwei Versuche mit veränderten thermischen aber konstanten aerodynami- 
schen Randbedingungen durchgeführt werden. Die Ergebnisse sind jeweils Partikulärlösungen 
der Energiegleichung der Grenzschicht und können mittels dem Superpositionsprinzip der Film- 
kühlung (siehe Metzger et al. (1971) und Choe et al. (1974)) zur Lösung des mathematischen 
Problems verwendet werden. 

Die Effektivität der Filmkühlung lässt sich dann durch die dimensionslose adiabate Filmkühl- 
effektivität (x T, wy), die dimensionslose Wandtemperatur (x Tw) sowie das dimensionslose Ver- 
hältnis der Wärmeübergänge mit und ohne Filmkühlung (x hf) vollständig beschreiben. Diese 
Größen werden von vielen Faktoren beeinflusst und sind daher auch selbst Thema wissenschaft- 
licher Untersuchungen. Hervorzuheben sind diejenigen Faktoren, die maßgeblich durch kom- 
plexe? Innenströmungen beeinflusst werden: 


e Das treibende Druckverhältnis IT = p.k/py zwischen Kühlluft und Heißgas. Aufgrund von 
Druckverlusten und Sekundärströmungen gilt zu beachten, dass sich der Totaldruck der 
Kühlluftströmung entlang des Strömungskanals ändert. 


e Die Geschwindigkeit und Richtung der internen Strömung relativ zur Bohrung, beschrie- 
ben durch den Winkel «æ des (lokalen) Strömungsvektors zur Bohrungshauptachse sowie 
der Machzahl Max der Kühlluft. 


3 komplex beschreibt Bedingungen, die von der Annahme der Ausströmung aus einem Plenum mit ruhendem 
Fluid abweichen. 
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Diese Faktoren wirken sich vor allem auf den Durchflusskoeffizienten Ca, das Dichteverhältnis 
DR, die Ausblaserate M und das Impulsverhältnis / 


N 2 
C=, DR=%, m=" und 1= PKK (2.4) 
Mid PH PHCH PHC 


aus, die untereinander gekoppelt sind und den Massendurchsatz durch eine Bohrung sowie die 
Interaktion des austretenden Strahls mit dem Heißgas beschreiben. 

Ein weiterer wichtiger Faktor ist die Geometrie der Bohrung. Diese beeinflusst das Durchfluss- 
und Kühlverhalten in starkem Maße. Durchgesetzt haben sich zylindrische Bohrungen aufgrund 
der einfachen Fertigung trotz mäßiger Strömungseigenschaften, und Diffusorbohrungen, soweit 
diese fertigungstechnisch möglich sind. Hier wird meist die laidback fanshaped-Konfiguration 
verwendet: Eine zylindrische Startbohrung, die in einen in zwei Achsen (lateral und vertikal) 
aufweitenden Diffusor mündet. 


Zum Durchfluss- und Kühlverhalten von Filmkühlbohrungen wurden vor allem am Institut für 
Thermische Strömungsmaschinen weitreichende Untersuchungen durchgeführt. Gritsch et al. 
(1998a,b,c) analysieren verschiedene Bohrungsgeometrien. Sie variieren das Druckverhältnis 
II und die interne Machzahl Max bei interner Parallelanströmung (a = 0). Aus umfangreichen 
Ergebnissen leiten sie Korrelationen für zylindrische Bohrungen ab. Sie zeigen außerdem, dass 
Diffusorbohrungen höhere Durchflusskoeffizienten erreichen, sowie weniger Variation dersel- 
ben bei veränderten Randbedingungen zeigen. 

Saumweber et al. (2001) führen Untersuchungen zum Durchflussverhalten und zur adiabaten 
Filmkühleffektivität bei verschiedenen Anströmwinkeln der Kühlluftströmung durch. Sie zei- 
gen, dass die Diffusorgeometrie zwar hohe Durchflusskoeffizienten erreicht, bei senkrechter 
Innenanströmung (a = 7/2) aufgrund stark gestörter Strömung im Diffusor der Bohrung jedoch 
eine inhomogene Verteilung der Filmkühleffektivität entsteht. Ein vergleichbares Experiment 
wird von Gritsch et al. (2003) durchgeführt. Auch hier wird die komplexe Wechselwirkung 
zwischen Bohrungsgeometrie, Anströmung und Ausblaserate beschrieben. Gritsch et al. (2003) 
bewerten zylindrische Bohrungen bei senkrechter Anströmung ebenfalls besser als Diffusorgeo- 
metrien. Auch aktuelle Untersuchungen von Fraas et al. (2017, 2018) an skalierten Bohrungen 
unter Variation des Anströmwinkels æ sowie der Machzahl Max bestätigen die bisherigen Er- 
gebnisse weitestgehend: Die verwendeten laidback fanshaped-Bohrungen zeigen eine starke 
Variation der adiabaten Filmkühleffektivität bei veränderter Anströmung. 


Als weiteren wichtigen Parameter identifizieren Bunker (2005) und Bogard et al. (2006) den so- 
genannten compound-Winkel ß, also den Winkel zwischen Heißgasströmung und Bohrungsach- 
se. In einer umfassenden Literaturrecherche fassen sie aerodynamische (Bunker, 2005) und ther- 
mische (Bogard et al., 2006) Untersuchungen der heißgasseitigen Strömungsstrukturen stromab 
von Filmkühlbohrungen zusammen. Die Ergebnisse zeigen wieder, dass die Wechselwirkung 
zwischen den geometrischen und den aerodynamischen Parametern hochkomplex und schwer 
vorherzusagen ist. Auch Lim et al. (2013) untersuchen den Einfluss der verschiedenen Anstell- 
winkel und stellen fest, dass ein Winkel zwischen Heißgas und Bohrungsachse für zylindrische 
Bohrungen für die adiabate Filmkühleffektivität vorteilhaft sein kann. 

Hervorzuheben ist außerdem eine Studie von zum Heißgas angewinkelten Bohrungen (zylin- 
drisch wie auch mit Diffusor mit ß > 0) von Ganzert et al. (2000), die mittels Ölfarben die 
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lokalen Strömungsstrukturen hinter den Bohrungen untersuchen. Die Ergebnisse, speziell ei- 
ne beobachtete Aufteilung des Kühlluftstrahles in zwei feine Einzelstrahlen bei zylindrischen 
Bohrungen, werden in ähnlicher Form in Kapitel 6 beschrieben. 


Zusammenfassend ist festzuhalten, dass eine genaue Prognose der Effektivität der Filmkühlung 
schwierig bleibt. Durch die inzwischen übliche Verwendung komplexer Innengeometrien in 
Verbindung mit dreidimensional konturierten Schaufeln und Seitenwänden (siehe auch Kneer 
et al., 2016) sind hohe und lokal variierende Anströmwinkel œ sowie compound-Winkel ß die 
Regel. Experimentelle Untersuchungen müssen sich meist auf die Variation weniger Parameter 
beschränken. Zwar existiert ein großer Datensatz zur Auslegung der Filmkühlung, eine Validie- 
rung eines konkreten Schaufelkühlsystems kann aber nur im Ganzen erfolgen. 


2.2.3 Keramische Schichten 


Falls konvektive oder kombinierte Methoden die Materialtemperatur nicht weit genug senken, 
können zusätzlich keramische Dämmschichten aufgebracht werden. Die Aufbringung einer sol- 
chen Schicht ist in Abbildung 2.2(b) dargestellt. Auf die zu kühlende Oberfläche wird eine ther- 
misch hoch resistente, keramische (meist ZrO>) Schicht aufgetragen. Aufgrund der geringen 
Wärmeleitfähigkeit Azro, der Schicht kann trotz geringer Schichtdicke eine deutliche Tempera- 
turabsenkung erreicht werden. 

Die komplexe Berechnung der Wärmeübergänge an Heißgas- und Kühlluftseite muss wie in 
Unterabschnitt 2.2.1 bzw. Unterabschnitt 2.2.2 durchgeführt werden. Der Temperaturverlauf im 
Festkörper ist bei gegebenem Wärmestrom und konstanter Wärmeleitfähigkeit der Materialien 
als abschnittsweise stetige Funktion mit den Steigungen « Er und «Ag : gegeben. 


2.3 Ubersicht aktueller Untersuchungen zur Erhéhung des internen War- 
meübergangs 


Das Ziel einer Variation der Innengeometrie ist die Erhöhung des absoluten, an die Kühlluft 
übertragenen, Wärmestroms. Dies kann bei gegebenen Temperaturniveaus durch eine Vergrö- 
Berung der Fläche sowie die Erhöhung des Wärmeübergangskoeffizienten h (bzw., in dimensi- 
onsloser Darstellung, der Nußeltzahl) erreicht werden. Dem Ziel einer großen Fläche und hohen 
Nußeltzahl steht der Druckverlust im Kühlsystem gegenüber, der möglichst gering sein soll. 
Aufgrund dieser Überlegungen werden die Ergebnisse von Untersuchungen meist als Kombi- 
nation aus relativer Erhöhung der Nußeltzahl Nu/Nuy und relativer Erhöhung des Druckverlust- 
faktors f/f) angegeben. Die folgenden drei Abschnitte geben einen Überblick über die bisher 
erzielten Verbesserungen durch Einsatz verschiedener Verfahren. Die Trennung kann nur grob 
erfolgen, da der Übergang zwischen den Verfahren fließend ist und willkürlich von den jewei- 
ligen Autoren festgelegt wird. Die Leistungsfähigkeit einiger ausgewählter Konzepte ist ab- 
schließend in Abbildung 2.5, Seite 24 gezeigt. Die Konzepte, welche im Rahmen dieser Arbeit 
untersucht werden, können nicht in diesem Diagramm eingefügt werden, da die Druckverlus- 
te nur für das gesamte interne System bestimmt werden können und damit deutlich über den 
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gezeigten Untersuchungen der Teilkonfigurationen liegen. Dies ist ein weiteres Zeichen für die 
Unzulänglichkeit der bisher gesammelten Daten für die Bestimmung der Effektivität moderner 
Konzepte. 


Rippen 


Rippen, orthogonal oder angestellt in Bezug auf die Hauptströmungsrichtung, erzeugen 
Rückströmgebiete, generieren Turbulenz und fördern das Mischverhalten der Strömung. 
Zusätzlich kann eine starke Flächenvergrößerung erreicht werden. Dies führt auch dazu, dass 
Rippen, je nach genauer Konfiguration, häufig hohe Druckverluste erzeugen. Aufgrund zahl- 
reicher geometrischer Parameter (Rippenhöhe bzw. Versperrung, Rippenbreite, Rippenabstand, 
Winkel, Form) lassen sich Rippen sehr gut an die gegebenen Anforderungen bezüglich f/fo 
bzw. Nu/nu, anpassen und finden so an verschiedensten Positionen Verwendung. 

Taslim et al. (1994) untersuchen eine Vielzahl von verschiedenen Geometrien in einem 
rechteckigen Kanal mit Rippen an zwei Seiten. Sie variieren die Rippenbreite, die Rippenhöhe, 
die Rippenform (orthogonal, angewinkelt, V-förmig) und den Versatz zwischen den Rippen auf 
zwei gegenüberliegenden Seiten. Der Rippenabstand wird in Bezug auf die Rippenhöhe kon- 
stant gehalten. Sie zeigen, dass angewinkelte Rippen, trotz geringer Versperrung, aufgrund der 
zusätzlichen Umlenkung der Strömung die höchste thermische Leistungsfähigkeit erreichen. 
V-förmige Rippen werden ebenfalls positiv bewertet. Orthogonale Rippen hingegen erzielen 
bei relativ hohen Druckverlusten eine nur geringe Steigerung des Wärmeübergangs. 

Mit der gleichen experimentellen Anlage untersuchen Taslim et al. (1997, 1999) außerdem den 
Effekt des Rippenabstandes sowie weitere Variationen von Rippenhöhen und Versperrungen. 
Sie zeigen, dass vor allem der Wärmeaustausch an den Rippen selbst, und nicht zwischen 
den Rippen, maßgeblich ist. Dies kann zu inhomogenen Wandtemperaturverteilungen führen. 
Außerdem steigt die Nußeltzahl aufgrund steigenden Anteils an Sekundärströmungen mit 
zunehmender Lauflänge stark an. Dieser Anstieg ist wieder abhängig von der Versperrung. 
Eine ideale Geometrie kann in keinem Experiment gefunden werden. So interagieren alle 
geometrischen Parameter, und selbst bei einfacher Kanalgeometrie ist eine Vorhersage der 
thermischen Leistungsfähigkeit schwierig. In der Anwendung nimmt die Kühllufttemperatur 
mit der Lauflänge ebenfalls zu, was dem Nutzen durch die steigende Nußeltzahl entgegenwirkt. 
Eine Vielzahl von weiteren geometrischen Variationen werden von Cho et al. (2001) betrachtet. 
In einem rechteckigen Kanal bewerten sie die Leistungsfähigkeit von ein- und zweiseitig 
angebrachten Rippen. Zusätzlich zu den üblichen Konfigurationen mit durchgängigen Rippen 
werden auch solche mit Unterbrechungen zwischen angrenzender Wand und Rippe (sog. 
diskrete Rippe) oder in der Rippenmitte betrachtet. Auch sie zeigen eine breite Spanne von 
möglichen Auslegungspunkten und empfehlen unterbrochene Rippen. Aufgrund komplexer, 
vom rechteckigen Querschnitt abweichender Kanäle in realistischen Schaufeln werden Rippen 
aktuell meist schon als diskrete Rippen ausgeführt. 

Neuere Studien von Salameh et al. (2012), durchgeführt in einem rechteckigen U-Kanal als 
erweitertes Modell einer internen Passage, bestätigen die bisherigen Studien weitestgehend. 
Wieder wird der maximale Wärmeübergang an der Rippe selbst gezeigt. Zwischen den Rippen 
stellt sich ein komplexes, von Sekundärströmungen geprägtes Nußeltzahlfeld ein. Es werden 
nur orthogonale Rippen untersucht. 
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Domaschke et al. (2012) bestimmen den Wärmeübergangskoeffizienten und den Druckverlust 
in glatten und berippten Kanälen. Der Querschnitt des Kanals ist einer realen Geometrie 
nachempfunden, ansonsten ist der Kanal prismatisch aufgebaut. Die Messung der internen 
Oberflächentemperaur erfolgt mittels kalibrierter, thermosensitiver Farbe. Domaschke et al. 
zeigen eine Erhöhung des Wärmeübergangskoeffizienten durch Rippen. Als kritischen Bereich 
mit niedrigem Wärmeübergangskoeffizienten identifizieren sie die Spitze des Kanals rückseitig 
des Staupunktes des Heißgases. 

Eine experimentelle und numerische Studie von Waidmann et al. (2013) untersucht ebenfalls 
berippte Kanäle mit realistischen Querschnitten. Zusätzlich wird ein U-Kanal eingesetzt. Die 
Bestimmung der Zielgrößen erfolgt wieder durch Einsatz von thermosensitiver Farbe und einer 
transienten Methode zur Ableitung des Wärmeübergangskoeffizienten. Sie zeigen eine lokal 
stark variierende Verteilung des Wärmeübergangskoeffizienten an den mit Rippen versehenen 
Wänden. Im Mittel kann eine Erhöhung erreicht werden. 

Untersuchungen von Rippen in Kanälen mit Geometrien ähnlich der in Schaufelvorder- 
kanten werden weiterhin von Saxer-Felici et al. (2014) durchgeführt. Sie empfehlen eine 
verhältnismäßig hohe (10%Dy), angestellte Rippe für einen optimalen Kompromiss aus 
Druckverlusterhöhung und Wärmeübergangserhöhung. Durch umfangreiche Variation von 
geometrischen Parametern können Saxer-Felici et al. außerdem Zusammenhänge zwischen die- 
sen und den Zielgrößen Wärmeübergangserhöhung und Druckverlusterhöhung identifizieren. 
Eine vergleichbare Studie von Naik et al. (2014) ergänzt die Ergebnisse von Saxer-Felici et al. 
mit Daten für den Hinterkantenbereich von Turbinenschaufel. Diese Kanäle zeichnen sich 
meist durch ein hohes Seitenverhältnis aus. Naik et al. zeigen, dass hier die Geometrie der 
Rippen zwar die lokalen Strömungsphänomene beeinflusst, die flachengemittelten Kenngrößen 
für verschiedene Rippengeometrien aber vergleichbar sind. 

Auch ausschließlich numerische Untersuchungen werden in jüngerer Vergangenheit immer 
öfter durchgeführt. Amano et al. (2010) führen beispielsweise Berechnungen mit orthogonalen 
wie auch angestellten Rippen durch. Aufgrund der komplexen Geometrien müssen Turbulenz- 
modelle verwendet werden. Die gewonnenen numerischen Daten weichen weit von den als 
Vergleich gezeigten experimentellen Daten ab. Zusätzlich führt die Wahl des Turbulenzmodells 
zu unsystematischen Abweichungen in der Größenordnung der Messgröße selbst. 


In modernen Studien wird zusätzlich der Effekt der Rotation mit berücksichtigt, da vor allem 
für die interne Strömung eine Beeinflussung zu erwarten ist. Zum Einsatz kommen rotierende 
U-Kanäle mit verschiedenen Querschnitten, da diese eine gute Repräsentation einer internen 
Kühlluftführung darstellen: Durch die zwei Hauptkanäle kann die Strömung radial nach au- 
ßen wie auch die Strömung radial nach innen beurteilt werden. Aufgrund der Berechnung der 
wirkenden Scheinkräfte nach Gleichung 2.1 werden in den jeweiligen Kanälen gegenläufige Ef- 
fekte vermutet. Diese Strömungsstrukturen werden von Elfert et al. (2008) im rotierenden und 
nicht rotierenden System mit und ohne Rippen in einem rechteckigen Kanal experimentell un- 
tersucht. Sie zeigen die Ausbildung von Corioliswirbeln bei Rotation, also jeweils einem Wirbel 
an der linken sowie rechten Innenwand (links-rechts-Richtung senkrecht auf Verbindungslinie 
Saugseite-Druckseite), und ein dadurch induziertes, weiteres, schwächeres, um 90° versetztes 
Wirbelpaar. Durch die Umlenkung am maximalen Radius bilden sich Dean-Wirbel, die sich 
mit den nun in umgekehrter Drehrichtung entstehenden Corioliswirbeln zu einem komplexen 
Sekundärströmungsfeld überlagern. Durch die eingebrachten Rippen wird eine Abschwächung 
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der Corioliswirbel beobachtet. Des Weiteren bildet sich ein asymmetrisches Profil der Radial- 
geschwindigkeit aus: Bei radialer Strömung nach außen steigt die Massenstromdichte an der 
druckseitigen Wand, bei radialer Strömung nach innen an der saugseitigen Wand. 

Fu et al. (2006) untersuchen 45°-angestellte Rippen in einem rotierenden U-Kanal. In einer 
umfangreichen Studie werden sowohl das Verhältnis von Kanalbreite zu Kanalhöhe (aspect 
ratio), als auch verschiedene Winkel zwischen den Hauptachsen des Kanalquerschnittes und 
der Rotationsrichtung analysiert. So lassen sich die Ergebnisse auf viele komplexe Kanäle in 
Turbinenschaufeln übertragen, da das Aspektverhältnis maßgeblich für den Einfluss der Corio- 
liswirbel ist. Fu et al. (2006) zeigen, dass die Rotation vor allem den lokalen Wärmeaustausch 
an den Wänden mit hoher lokaler Radialgeschwindigkeit der Kühlluft, wie von Elfert et al. 
(2008) beschrieben, fördert, während die Wände mit niedriger lokaler Radialgeschwindigkeit 
der Kühlluft geringere Nußeltzahlen zeigen. Dieser Effekt verstärkt sich mit steigender Rota- 
tionszahl Ro. Global zeigen sie eine leichte Erhöhung der thermischen Leistungsfähigkeit bei 
Rotation für glatte sowie verrippte Kanäle. 

Experimente mit hohen Rotationszahlen in Kanälen mit dreieckigem Querschnitt und Rippen 
wurden von Liu et al. (2010) durchgeführt. Diese Kanalform eignet sich, da sie gut in die spitz 
zulaufende Geometrie der Schaufelvorderkante zu integrieren ist. Auch sie zeigen eine generel- 
le Erhöhung der Nußeltzahl bei Rotation, wobei der Effekt bei verschiedenen Rippengeometrien 
unterschiedlich stark ausgeprägt ist. 

Numerische Simulationen mit modernen Methoden sowie komplementäre experimentelle Stu- 
dien mit Rotation wurden von Singh et al. (2017) durchgeführt. Hier zeigt sich ebenfalls eine 
Abweichung zwischen Experiment und numerischer Vorhersage, wobei Trends richtig vorher- 
gesagt werden. Sie beschreiben wieder vergleichbare, durch Rotation induzierte Strömungs- 
strukturen und den generellen Anstieg der Nußeltzahl. 


Dellen und Erhöhungen 


Dellen und Erhöhungen (meist in Kreis- oder Tropfenform) erzeugen zwei gegenläufig 
rotierende, kleinere Wirbel, die advektiv in Strömungsrichtung transportiert werden und 
dissipieren. Ligrani et al. (2001) beschreiben diese Wirbelentstehung als periodisch und mit 
starkem Mischverhalten, was sich positiv auf den Wärmeübergang auswirkt. Aufgrund der 
meist geringen Versperrung ist der Druckverlust gering. Ein typischer Einsatzbereich sind 
daher oft die hinteren Turbinenstufen, bei denen ein niedriger Kühlluftdruck zur Verfügung 
steht. Ein weiterer Vorteil von Dellen gegenüber Erhöhungen (und anderen hier vorgestellten 
Verfahren) ist die Materialentfernung im Gegensatz zu einer Materialaufbringung. Dies kann 
leichtere Bauteile ermöglichen. 

Chyu et al. (1997) untersuchen solche Dellen mit Kreis- und Tränenform. Sie zeigen sowohl 
lokal aufgelöste als auch gemittelte Nußeltzahlverteilungen für Dellen bei verschiedenen 
Reynoldszahlen. Die beschriebenen Konfigurationen sind Kombinationen aus entweder 
kreisförmigen oder tropfenförmigen Dellen an zwei gegenüberliegenden Wänden in einem 
rechteckigen Kanal. Chyu et al. zeigen, dass die von ihnen verwendeten Konfigurationen eine 
mit Rippen vergleichbare Erhöhung des Wärmeübergangs erzielen können, dabei aber deutlich 
niedrigere Druckverluste (~ Faktor zwei) erzeugen. 
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Mahmood et al. (2001) führen vergleichbare experimentelle Untersuchungen durch. In einem 
rechteckigen Kanal werden auf einer Seite kreisförmige Dellen eingebracht, während auf der 
gegenüberliegenden Seite, an gleicher Position, versetzt, oder keine kreisförmigen Erhöhun- 
gen eingebracht werden. Sie können den vergleichsweise niedrigen Druckverlust bei modera- 
ter Erhöhung der Nußeltzahl bestätigen. In weiteren Versuchen untersuchen Mahmood et al. 
(2002) anhand der selben Geometrien die Wechselwirkungen von Aspektverhältnis, Tempera- 
turverhältnis zwischen Fluid und Wand und der Reynoldszahl. Sie betonen die Effektivität von 
Dellen bei flachen Kanälen und eine Steigerung der thermischen Leistungsfähigkeit bei Tem- 
peraturverhältnissen, die wesentlich kleiner als eins sind. 

Geometrische Variationen von Dellen werden von S. W. Moon et al. (2002) untersucht. Neben 
der meist verwendeten halbkreisförmigen Ausführung kommen auch zylindrische Aussparun- 
gen zum Einsatz, die sich durch scharfe rück- und vorspringende Stufen auszeichnen. In dem 
von S. W. Moon et al. (2002) realisierten Bereich der Reynoldszahl (turbulent, aber vergleichs- 
weise niedrig) bewerten sie die zylindrischen Strukturen als leistungsfähiger gegenüber den 
kreisförmigen Dellen. Eine Fertigung von scharfen, feinen Strukturen ist aber im für Schaufeln 
meist verwendeten Feingussprozess problematisch. 


Zylinder, Pins und ähnliche Einbauten 


Zylinder sowie Pins sind senkrecht zur Strömung stehende zylindrische Strukturen. Alternativ 
können auch diamant- oder würfelförmige Hindernisse verwendet werden. Neben der eigentli- 
chen Ausformung der Körper ist auch deren Anordnung entscheidend, hier wird zwischen einer 
gleichmäßigen Anordnung im Gitter und einer versetzten Anordnung unterschieden. Je nach 
Versperrung werden meist hohe Wärmeübergänge bei hohen Druckverlusten und einer starken 
Flächenvergrößerung erreicht. Typische Einsatzgebiete ergeben sich daher zur Kühlung von 
Schaufelhinterkanten. 

Die meisten Untersuchungen beschränken sich auf die thermische Leistungsfähigkeit, das aero- 
dynamische Verhalten wird wenig beschrieben. Ligrani et al. (2003) fassen die aerodynami- 
schen Strukturen zusammen: So bilden sich an jeder Vorderkante Hufeisenwirbel, die um den 
Zylinder transportiert werden. Hinter jedem Zylinder entstehen Wirbelstraßen, die dann mitein- 
ander interagieren. Die Ausbildung einer jeden Wirbelstraße ist abhängig von den Strömungsbe- 
dingungen und der Lage des jeweiligen Zylinders im Feld. So ergebenen sich mit zunehmender 
Lauflänge beliebig komplexe Sekundärströmungen, die aufgrund der hohen Versperrung zu ho- 
hem Druckverlust und hohem Wärmeübergang führen. 

Chyu (1990) untersucht solche Felder mit Zylindern bei verschiedenen Reynoldszahlen. Au- 
ßerdem variiert er die Anordnung und die Form: Durch Einbringen einer Verrundung am Zylin- 
derfuß kann der sich ausbildende Hufeisenwirbel gezielt beeinflusst werden. Er zeigt, dass die 
Einbringung der Verrundung geringen Einfluss hat, wobei die Anordnung des Feldes hingegen 
großen Einfluss auf die Nußeltzahl hat. Chyu leitet außerdem ab, dass der Wärmeübergangs- 
koeffizient auf den Zylindern selbst wie auch auf den Wänden vergleichbar ist. So kann die 
Flächenvergrößerung durch die Zylinder voll genutzt werden. Die erzeugten Druckverluste sind 
hoch. 

In neueren Untersuchungen ergänzen Chyu et al. (1998) die vorher gewonnenen Daten durch 
neue Formen der Einbauten. Sie verwenden Würfel sowie Diamanten (um 45° rotierte Würfel) 
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in regelmäßiger und versetzter Anordnung. Im Vergleich zu den vorher untersuchten Zylindern 
erreichen Würfel ähnliche Werte. Diamanten hingegen erreichen schon bei gleichmäßiger An- 
ordnung ähnliche Druckverlust- und Nußeltzahlerhöhungen wie die versetzt angeordneten Zy- 
linder. Eine versetzte Anordnung kann beide Kenngrößen weiter steigern, wobei die Erhöhung 
der Nußeltzahl deutlich geringer als die Erhöhung des Druckverlustes ausfällt. Abschließend 
bewerten sie die gleichmäßig angeordneten Würfel als am besten geeignet. 

Hwang et al. (2001) untersuchen die Kühlung einer Turbinenhinterkante durch versetzte Ein- 
bauten mit zylindrischer, Würfel- und Diamantform bei verschiedenen Reynoldszahlen in einem 
trapezförmigen Kanal. Im Vergleich zu anderen Untersuchungen realisieren sie einen graduel- 
len Masseverlust (bleed) durch Abnahme an einer Seite, wie er in internen Kühlkanälen durch 
Filmkühlung meistens vorliegt. Der Einfluss der so entstehenden Querströmung wird quantifi- 
ziert. Sie zeigen ein Minimum der Nußeltzahl, wenn die durch die Querströmung entstehende 
Strömungshauptrichtung nahezu der Versatzrichtung der Elemente gleicht. Des Weiteren kön- 
nen sie quantitativ ähnliche Resultate wie die vorhergehenden Studien zeigen. 
Untersuchungen an Turbinenschaufelhinterkanten mit Einbauten wurden auch am Institut für 
Thermische Strömungsmaschinen durchgeführt. Martini et al. verwenden sowohl experimen- 
telle (Martini et al., 2005) als auch numerische Methoden (Martini et al., 2006) zur Bestim- 
mung des Wärmeübergangskoeffizienten und der adiabaten Filmkühleffektivität an zurückge- 
schnittenen Hinterkanten mit Rippen sowie Zylinderfeldern. Sie konzentrieren sich auf den Be- 
reich nach der Ausblasung aus dem Spalt an der Schaufelhinterkante. Martini et al. können 
zeigen, dass die den Wärmeübergang dominierende turbulente Mischung am Rückschnitt durch 
Einbauten nahe der Ausblasung positiv beeinflusst werden kann. 

In neueren Untersuchungen versuchen Siw et al. (2012) die Form der Elemente selbst weiter zu 
optimieren. Durch Verwendung von Halbzylindern und Dreiecken (,,Halb-Diamanten“), die sie 
mit Vollzylindern vergleichen, wollen sie das Gewicht der Konfigurationen reduzieren. Die ex- 
perimentellen Studien werden mit verschiedenen versetzten Anordnungen der drei Geometrien 
bei verschiedenen Reynoldszahlen durchgeführt. Siw et al. (2012) zeigen eine hohe thermische 
Leistungsfähigkeit für die dreieckigen Elemente. Auch die halbzylindrischen Elemente können 
gut bewertet werden: Zwar fällt die erreichte Erhöhung des Wärmeübergangs geringer aus, die 
hier untersuchten Elemente benötigen jedoch auch eine wesentlich kleinere Fläche auf der End- 
wand. So kann eine dichtere Besetzung und damit stärkere Wärmeübergangserhöhung erreicht 
werden, sofern der zusätzliche Druckverlust in Kauf genommen wird. 


2.3.1 Drallkühlung 


Bei der Verwendung von Drallkühlung wird, meist im Schaufelfuß, ein Drallerzeuger in Vor- 
derkantennähe eingesetzt. Die Drallerzeugung erfolgt durch tangentiale Einblasung oder axiale 
Leitgitter (siehe auch Abbildung 3.3, Seite 42). Diese Drallerzeuger generieren die in Abbil- 
dung 2.4(a) gezeigte Strömungsstruktur mit positiven Eigenschaften für den internen Wärme- 
übergang. Die Erhöhung der Nußeltzahl wird maßgeblich durch zwei Gegebenheiten erreicht. 
Bei stationärem Drall gilt das Gleichgewicht zwischen Druckkräften und Zentrifugalkräften 


2P pen 
or 


(2.5) 
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das zu einem niedrigen Druck im Kern der Strömung führt. Dieser Druckabfall ist direkt nach 
dem Drallerzeuger bei höchster Umfangsgeschwindigkeit maximal und nimmt in Strömungs- 
richtung ab. Dieser positive Druckgradient im Kern bedingt ein Rückströmgebiet und damit, 
bei gleichem durchgesetzten Massenstrom, eine erhöhte Axialgeschwindigkeit in Wandnähe. 
In Verbindung mit der hohen Umfangskomponente der Geschwindigkeit ergibt sich eine be- 
tragsmäßig hohe Geschwindigkeit nahe der Wand, was Mischung, Turbulenz und damit den 
Wärmeübergang fördert. Eine mögliche Version eines als Gitter ausgeführten Drallerzeugers 
mit entstehendem Rückströmgebiet ist in Abbildung 2.4(a) gezeigt. Dieses Rückströmgebiet ist 
anhand des Verlaufes der Axialgeschwindigkeit über den Kanalradius auch in Abbildung 2.4(b) 
links dargestellt. In Abbildung 2.4(b) rechts ist der Verlauf der Umfangskomponente der Ge- 
schwindigkeit abgebildet. Dieser ergibt sich als Annäherung an einen Rankine-Wirbel (gestri- 
chelter Verlauf) mit Starrkörperwirbelzone (I) und Potentialwirbelzone (II). Durch die sich aus- 
bildende Grenzschicht an der Außenwand fällt die Umfangsgeschwindigkeit plötzlich ab. Dies 
führt zu lokalen Druckgradienten in die Grenzschicht und damit zum Aufplatzen dieser. Diese 
Durchmischung der Grenzschicht wirkt sich positiv auf den Wärmeübergang aus. 

Wie auch bei den klassischen Verfahren lässt sich die Druckverlusterhöhung als Aufwand des 
Kühlverfahrens definieren. Diese ist vor allem lokal auf die drallerzeugende Geometrie be- 
schränkt. Da aber bei den klassischen Verfahren auch die integrale Erhöhung des Druckverlus- 
tes angegeben wird, bleiben die auf das Gesamtsystem bezogenen Ergebnisse vergleichbar. Als 
Nutzen wird wieder die Erhöhung der Nußeltzahl verwendet. Drallkonzepte, die in der Literatur 
untersucht wurden und von denen diese zwei Parameter zur Verfügung gestellt werden, sind in 
Abbildung 2.5 mit größeren Symbolen dargestellt. 


Das aerodynamische Verhalten von Drallströmungen wird schon früh von Hall (1972) be- 
schrieben. Er konzentriert sich vor allem auf die Ausbildung und das Verhalten des axialen 
Rückströmgebietes, zu erkennen in Abbildung 2.4(b), linke Seite. Maßgeblich für dessen Ent- 
stehung ist demnach ein hoher Drallwinkel 9 = arctan(“/u.x) von 9 Z 40° direkt vor dem 
Entstehungsort und ein positiver Druckgradient in axialer Richtung. Aufgrund der progres- 
siven Instabilität der Drallströmung kann ein Rückströmgebiet auch bei nicht vorhandenem 
oder sogar leicht negativem Druckgradienten auftreten. Die Versperrung durch dieses sorgt 
für einen Druckanstieg und dadurch eine Repositionierung und Stabilisierung des Rückström- 
gebietes. Ein stark negativer Druckgradient, z.B. durch konvergente Kanäle, muss aber ver- 
mieden werden. Jacobs et al. (1975) führen erstmals die Drallzahl* Y = M/r.K ein. Diese be- 
schreibt die Drallintensität durch das integrale Verhältnis zwischen Drehmoment des Impuls- 
flusses M = je Püyltaxr- dr sowie dem Impulsstrom K = Hh pu2,rdr und wird mithilfe des Ra- 
dius r entdimensioniert. Die Autoren beschreiben einen Bereich von einer minimal nötigen 
Drallzahl von Y = 0,3...0,6 zur Ausbildung eines stabilen Rückströmgebietes. Kleine Werte 
sind nur ausreichend bei Strömungen mit starkem positivem Druckgradient. Bei Drallzahlen 
Y > 0,6 beobachten sie außerdem periodische Störungen in Form von umlaufenden Wirbelfä- 
den um den Hauptkern. Diese Wirbelfäden unterstützen aufgrund ihrer hohen Turbulenzgrade 
das Mischverhalten der Strömung. 


* orig. Drallparameter S, aus Gründen der Konsistenz mit späteren Veröffentlichungen wird die gängigere Be- 
zeichnung der Drallzahl verwendet 
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Abbildung 2.4: Erhöhung des internen Wärmeübergangs mittels Drallkühlung 


Chang et al. (1994) untersuchen das lokal aufgelöste Stromungsfeld einer Drallströmung (Drall- 
erzeugung durch tangentiale Einblasung) und zeigen die mögliche Aufteilung des entstehenden 
Wirbels in die zwei Zonen Starrkörperwirbel I und Potentialwirbel II. Des Weiteren stellen 
sie erste Korrelationen für den Abfall der Drallzahl mit der Lauflänge (swirl decay) auf. Da 
der erwartete Wärmeübergang stark von der Drallzahl abhängt, können diese Korrelationen für 
Vorauslegungen genutzt werden. 

Mit weit ausgereifter Messtechnik können Grundmann et al. (2012) die vorangegangenen Un- 
tersuchungen anhand eines Rohres mit tangentialer Einströmung zur Drallerzeugung bestätigen. 
Sie zeigen die erhöhte Axialgeschwindigkeit an der Wand und deren inhomogene Verteilung. 
So wandern Bereiche stark erhöhter Axialgeschwindigkeit helixförmig in Strömungsrichtung 
(sog. wall-jets). Eine generelle Erhöhung der wandnahen Axialgeschwindigkeit kann ebenfalls 
beobachtet werden (laut Quelle werden bis zu = 98% der Masse in Wandnähe transportiert). 
Ebenso können die Wirbelfäden gezeigt werden. Sie entstehen durch Interaktion der beiden 
Wirbeltypen I und II (siehe Abbildung 2.4(b)). Diese Strukturen wandern in einer Helix um 
den Hauptkern der Drallströmung und zeichnen sich durch eine positive axiale Geschwindigkeit 
aus. Sie werden von Grundmann et al. (2012) erstmals experimentell gezeigt und als core-jet 
bezeichnet. 

Diese Untersuchungen werden von Wassermann et al. (2012) ergänzt: Das einfache Rohr wird 
durch eine U-Rohr-Geometrie ersetzt, die näher an einer einsatzfähigen Geometrie liegt. Sie 
zeigen, dass durch Austrittsblenden und gezieltes Anpassen der 180°-Umlenkung die Wirbel- 
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strukturen und damit der Wärmeübergang gezielt beeinflusst werden können. Dies bezieht sich 
nur auf den ersten Teil des Kanals: Nach der Umlenkung bricht das Rezirkulationsgebiet wei- 
testgehend zusammen, ein Abfall des Wärmeübergangs wird vermutet. 

Bruschewski et al. (2015) ergänzen diese Untersuchungen mit weiteren geometrischen Variatio- 
nen des Kanals. So werden verschiedene Austrittsblenden wie auch ein 180°-Austrittskrümmer 
bei verschiedenen Drall- und Reynoldszahlen verwendet. Sie zeigen eine Veränderung der Strö- 
mungsstrukturen bei Variation der Austrittsgeometrie, können aber das für die Erhöhung des 
Wärmeübergangs wichtige Rückströmgebiet unabhängig davon zeigen. Weiterhin legen sie dar, 
dass die Drallströmung auch nach der 180°-Umlenkung stabil bleiben kann. Bruschewski et 
al. definieren einen dafür notwendigen minimalen Radius der Umlenkung, der sowohl von der 
Drallzahl als auch von der Reynoldszahl abhängt. 

Biegger et al. (2013) untersuchen eine durch tangentiale Einblasung erzeugte Drallströmung 
in einem Rohr mit kreisförmigem Querschnitt mittels 3D-PIV. Auch sie bestätigen das kom- 
plexe Strömungsfeld mit den in Abbildung 2.4 gezeigten Geschwindigkeitsverteilungen und 
Wirbelzonen I und II. Weiterhin können Biegger et al. durch instationäre Datenerfassung die 
instantanen Strömungsstrukturen auflösen und vermuten, dass die hoch turbulente Strömung in 
Wandnähe einen hohen Wärmeübergang fördert. 

Die durch die radiale Beschleunigung induzierten Grenzschichtinstabilitäten können aufgrund 
der zu geringen Auflösung der experimentellen Messungen in der Grenzschicht in diesen mo- 
dernen Untersuchungen nicht gezeigt werden. Saric (1994) fasst bisherige Untersuchungen zu 
dieser Instabilität zusammen. Ausgehend vom Rayleigh circulation criterion? , das diese Insta- 
bilitäten bei reibungsfreier Strömung beschreibt und notwendiges Kriterium auch bei reibungs- 
behafteten Strömungen ist, zeigt Saric die theoretischen Zusammenhänge der Grenzschichtmi- 
schung. An konkaven Wänden entstehen durch das Aufbrechen und Mischen der Grenzschicht 
gegenläufig rotierende Wirbelzöpfe, die sich in Umfangsrichtung erstrecken. Diese werden als 
Görtlerwirbel bezeichnet (nach Görtler (1941)). Sowohl in der Zusammenfassung von Saric 
(1994) wie auch in einer weiteren Veröffentlichung von Peerhossaini et al. (1988) können diese 
Strukturen experimentell nachgewiesen werden. Letztere zeigen, dass diese Strukturen auch bei 
turbulenten Freiströmungen stabil sind. 


Die beschriebenen aerodynamischen Strukturen legen nahe, dass stark verdrallte Strömungen 
positive Auswirkungen auf den Wärmeübergang an der Kanalwand haben. So zeigen Chang 
et al. (1995), anknüpfend an die durchgeführten aerodynamischen Messungen, thermische 
Ergebnisse. Sie leiten neue Korrelationen zur Bestimmung der Wärmeübergangserhöhung 
in Abhängigkeit von der Drallzahl ab. Diese Erhöhungen liegen in ähnlichen Bereichen 
wie bei klassischen Verfahren. Auch Glezer et al. (1996) verwenden eine Konfiguration mit 
tangentialer Einblasung, streben aber erstmals die Verwendung als Kühlverfahren für Turbi- 
nenschaufeln an. Erste Untersuchungen an Kanälen ähnlich denen in einer Turbinenschaufel 
zeigen vielversprechende Resultate mit signifikanten Wärmeübergangserhöhungen (bis zu 
Nu/Nuo < 6). Als Besonderheit untersuchen sie eine Konfiguration, bei der die gesamte Kühlluft 
über Filmkühlbohrungen abgegeben wird. Trotz der graduellen Massenentnahme kann ein 
stabiler Drall erreicht werden. Die Wärmeübergangserhöhungen fallen aufgrund der fehlenden 


5 Rayleigh postuliert, dass bei einer Strömung in Zylinderkoordinaten mit Z = (0, cy, 0),.o,; Instabilitäten auftre- 
ten, sobald ir: < 0 gilt. Hierbei ist die Zirkulation I" definiert als I = rcu. 
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Axialströmung geringer aus. Aussagen über den Druckverlust werden für keine Konfiguration 
getätigt. 

Hedlund et al. (1999, 2000) untersuchen eine Konfiguration mit zweifacher tangentialer Ein- 
blasung bei hohen Drallzahlen (Y > 10). Sie können sowohl die aerodynamischen Strukturen 
der Görtlerwirbel als auch lokal aufgelöste Nußeltzahlverteilungen zeigen. Global gemittelte 
Daten sind in guter Übereinstimmung mit den vorher durchgeführten Untersuchungen. 
Khalatov et al. (2000) führen experimentelle Studien an durch Leitgitter erzeugten Drall- 
stromungen zu Druckverlust und Nußeltzahlerhöhung durch. Sie zeigen hohe erreichbare 
Steigerungen des Wärmeübergangs bei einer starken Erhöhung des Druckverlustes. Trotzdem 
sprechen sie eine Empfehlung für das Konzept aus, da bei im Vergleich zu konventio- 
nellen Konzepten vereinfachten Geometrien ähnliche Wärmeübergänge erreicht werden. 
Khalatov et al. (2001) untersuchen außerdem Drall-Jet-Verfahren in einem Kanal, wie er 
geometrisch auch in einer Schaufelvorderkante implementiert werden kann. Gleichzeitig 
werden maschinenähnliche Parameter eingestellt (stat. Druck, Reynoldszahl, Wandtemperatur- 
zu-Kühllufttemperaturverhältnis). Sie können die vorher gewonnenen Ergebnisse trotz der 
komplexeren Einsatzbedingungen bestätigen. In einer neueren Studie von Khalatov et al. 
(2011) untersuchen sie erneut eine Drallkonfiguration mit tangentialer Einblasung in einem 
Rohr mit einer 90°-Krümmung vor dem Auslass. Sie zeigen, dass die Krümmung des Kanals 
auf den Drall stromauf wenig Einfluss hat. Sie bewerten die untersuchten Konfigurationen 
positiv: Zwar ist die erreichte thermische Leistungsfähigkeit ähnlich der von Rippen, die 
Fertigung der Kanäle ist aber deutlich einfacher. Die Thematik der Fertigung wird auch von 
Bunker (2009) behandelt. In einer Sensitivitätsstudie zeigt er, dass schon kleine Veränderungen 
der internen Geometrien, wie zum Beispiel eine Änderung des Rippenwinkels um +5° oder der 
Rippenhöhe um 20%, den Wärmeübergang maßgeblich beeinflussen. Diese Abweichungen 
sind aber nicht zu vermeiden. Eine Vereinfachung der Kanalgeometrie durch Verwendung der 
Drallkühlung kann hier Vorteile haben. 

Ahmadvand et al. (2010) verwenden axiale Leitgitter mit verschiedenen geometrischen 
Drallwinkeln zur Drallerzeugung in einem geraden Rohr. Auch sie zeigen eine Steigerung 
des Wärmeübergangs und des Druckverlustes, die maßgeblich von der Drallzahl, aber kaum 
von der Reynoldszahl abhängig ist. Die erreichten Nußeltzahlverhältnisse liegen bei ~ 2, 
während Druckverlustverhältnisse bis ~ 6 beobachtet werden. Diese nicht ausreichende 
Leistungsfähigkeit ist auf die von Ahmadvand et al. (2010) verwendete, einfache Ausführung 
der Drallerzeuger mit nicht konturierten Leitblechen zurückzuführen. 

Rao et al. (2017) führen experimentelle und numerische Untersuchungen an einem geraden 
Rohr mit einzelner und wiederholter tangentialer Einblasung durch. Sie zeigen, dass durch 
lokale Wirbelstrukturen in Wandnähe hohe Gradienten in der Verteilung der Nusseltzahl 
auftreten können. Nach Mittelung über die jeweiligen Umfänge an jedem Querschnitt verbleibt 
eine Steigerung der Nusseltzahl. Bei gleichem Massenstrom erreicht das Rohr mit einzelner 
Einblasung höhere gemittelte Nusseltzahlen als das Rohr mit wiederholter Einblasung. 
Dagegen ist bei einzelner Einblasung die Verteilung der Nusseltzahl inhomogener und der 
Druckverlust deutlich höher. 


Der Effekt der Rotation muss auch bei der Verwendung von Drallkühlung berücksichtigt wer- 
den. Winter et al. (2009) können zeigen, dass die Drallkühlung wesentlich stabiler bezüglich 
des Rotationseinflusses ist als eine rein axiale Rohrströmung. So zeigt sich die Ausbildung der 
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Corioliswirbel und der damit verbundenen saug- und druckseitigen Veränderung der Nußelt- 
zahlen (siehe z. B. Abschnitt Rippen) erst bei hohen Rotationszahlen. 

Hier ist mit einer Verschiebung der Lage der lokalen Extrema der Nusseltzahlverteilung zu 
rechnen. Diese Ergebnisse können durch weitere Untersuchungen von Winter (2010) bestätigt 
werden. 

Dontenville (2015) führt numerische Berechnungen der im Rahmen dieser Arbeit verwendeten 
Geometrien im rotierenden System durch. Er zeigt ebenfalls, dass die Rotation einen gerin- 
gen Einfluss auf den Wärmeübergang bei der Verwendung von Drallkühlung hat. Bei hohen 
Rotationszahlen bilden sich Corioliswirbel und damit lokale Verschiebungen der Nußeltzahl- 
verteilung, die gesamte übertragene Wärme des Systems wird wenig beeinflusst. 

Somit kann die Drallkühlung unter dem Aspekt der Rotation, also der Verwendung in Rotor- 
schaufeln, positiv bewertet werden. Es ergeben sich weniger Inhomogenitäten als bei klassi- 
schen Verfahren. Die Gesamtkühlwirkung wird nicht beeinflusst, eine Prognose der Leistungs- 
fähigkeit des Kühlsystems wird vereinfacht. 


Wechselwirkung zwischen Drall- und Filmkühlung 


Durch den Einsatz der Drallkühlung ergeben sich lokal hohe Umfangsgeschwindigkeiten in 
den internen Kanälen. Diese können zu einer Beeinflussung der Filmkühlung, wie in Unterab- 
schnitt 2.2.2 beschrieben, führen. Zu diesem komplexen Sachverhalt gibt es wenige Untersu- 
chungen, nennenswert sind die Arbeiten von Lerch et al. (2011a), Lerch et al. (2011b) und die 
dazugehörige, umfangreiche Dissertationsschrift von Lerch (2012). 

In einer ersten Studie (Lerch et al., 2011a) wird ein Kanal mit kreisförmigem Querschnitt in 
die Vorderkante eines symmetrischen NACA-Profils integriert. In diesem Kanal kann entweder 
eine rein axiale Strömung oder eine Drallströmung durch tangentiale Einblasung im Schaufel- 
fuß realisiert werden. Die Drallrichtung ist positiv um die Hochachse. Die Filmkühlung wird 
aus diesem Kanal gespeist, die Bohrungsachsen schließen einen spitzen Winkel mit der Rich- 
tung der lokalen Umfangsgeschwindigkeit ein. Mutmaßlich führt diese Konfiguration zu den 
niedrigsten Durchflusskoeffizienten Cp aufgrund der Richtungsumkehr der Luft bei Einströ- 
mung in die Bohrung. Darauf aufbauend (Lerch et al., 2011b) werden weitere geometrische 
Variationen untersucht: Es wird eine weitere Reihe von Filmkühlbohrungen an der gegenüber- 
liegenden Wand implementiert, deren Einströmwinkel stumpf ist. Außerdem wird mittels über 
die Schaufelhöhe verteilten Einblasebohrungen in den Vorderkantenkanal eine weitere Metho- 
de zur Drallerzeugung realisiert und so eine große Anzahl von Konfigurationen untersucht. Des 
Weiteren werden beide Drallrichtungen untersucht. Die Ergebnisse der Untersuchungen von 
Lerch et al. können wie folgt zusammengefasst werden: 


e Das für den Wandwärmeübergang wichtige Rezirkulationsgebiet bleibt bei allen Konfi- 
gurationen trotz Massenabnahme durch die Filmkühlbohrungen stabil. 


e Die Drallzahl steigt über die Schaufelhöhe an, da die Axialgeschwindigkeit durch Mas- 
senabfuhr abnimmt. Dies ist durch die Dichteänderung bei diabaten Messungen nicht 
übertragbar auf die späteren Versuche. 
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e Unterschiedliche Drallrichtungen beeinflussen das Strömungsfeld der Filmkühlung bis 
auf eine nötige Spiegelung nicht. 


e Die positive Drallrichtung bei einer Bohrungsreihe mit stumpfen Winkeln (Lerch et al., 
2011b) wird positiv bewertet (vergleichbare adiabate Filmkühleffektivität wie ohne 
Drall), bei negativer Richtung und damit spitzen Einströmwinkeln wird ein starker 
Einbruch der adiabaten Filmkühleffektivität beobachtet. 


e Durch Versuche mit gegenüberliegenden Bohrungsreihen stellen Lerch et al. eine Mas- 
senstromverschiebung zu den Bohrungen mit positivem Drallwinkel fest. Dies wiederum 
führt zu einer Verbesserung der adiabaten Filmkühleffektivität der negativ angeströmten 
Bohrungen, da durch die verringerte Ausblaserate die Interaktion und Mischung mit der 
Freiströmung reduziert wird. 


Es ergibt sich demnach eine komplexe Wechselwirkung zwischen Drall- und Filmkühlung, wel- 
che die Gesamteffektivität des Kühlsystems bestimmt. Die wichtigsten Kriterien sind hier der 
lokale Einströmwinkel (Beeinflussung von Cp), der Druckgradient über die Bohrungen (bei 
konstantem Kanaldruck und nicht symmetrischem Profil höher für die saugseitigen Bohrungen) 
sowie daraus resultierend das lokale Impulsverhältnis / für jede Bohrung. Durch geschicktes 
Abstimmen sehen Lerch et al. Potenzial für ein hoch effektives Kühlsystem. Sie merken aber 
an, dass aufgrund der Komplexität des Systems weitere Untersuchungen notwendig sind, die 
reale Geometrien der Kanäle, Schaufelprofile und Filmkühlbohrungen berücksichtigen. Weitere 
Untersuchungen, die sowohl diese realen Geometrien als auch realistische thermische Randbe- 
dingungen berücksichtigen, sind zum bisherigen Zeitpunkt nicht bekannt. 


Drallkühlkonzepte eignen sich demnach vor allem für die Vorderkante von Turbinenschaufeln 
der ersten oder zweiten Stufe, speziell von Flugzeugtriebwerken. So muss für die Erzeugung 
der kinetischen Energie der Drallströmung ein hohes Maß an innerer Energie aufgewendet wer- 
den, der Druck der Kühlluft muss hoch sein. Aufgrund starker Umlenkung ist die Erzeugung 
einer hohen Drallzahl weiterhin mit hohen Totaldruckverlusten verbunden. So können solche 
Drallerzeuger nur eingesetzt werden, wenn der Totaldruck der Kühlluft maßgeblich höher als 
der statische Druck des Heißgases ist, was gerade bei Flugzeugtriebwerken für die erste Rotor- 
stufe mit niedrigem Reaktionsgrad gilt. 

Die bisherigen Konzepte erreichen ähnliche thermische Leistungsfähigkeiten wie konventio- 
nelle Konzepte. Jedoch werden in den vorliegenden Studien Abstriche bei der Komplexität der 
drallerzeugenden Komponenten gemacht. So wurden bisher nur tangentiale Einblasungen oder 
einfachste, nicht konturierte, axiale Leitgitter untersucht. Dies steht im Widerspruch zu der 
komplexen, hoch optimierten Gestaltung der Leitgitter der Heißgasströmung. Durch einen Wis- 
senstransfer und die dreidimensionale Optimierung der Drallerzeuger für die interne Kühlung 
ist mit einer deutlichen Verbesserung der thermischen Leistungsfähigkeit zu rechnen. 
Aufgrund der umfangreichen Arbeiten auf dem Gebiet der internen Kühlverfahren kann die 
hier gegebene Zusammenfassung nur als Überblick dienen. Ein weiterer, umfangreicher Über- 
blick der verwendeten Verfahren ist in Ligrani et al. (2003) und als Erweiterung dazu in Li- 
grani (2013) zu finden. Während die erste Zusammenfassung klassische Verfahren und erste 
Drallkühlverfahren bis 2003 untersucht, werden in der zweiten Zusammenfassung klassische 
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Zyl. gleichmaBig, Chyu (1990) 

abger. Zyl. gleichmaBig, Chyu (1990) 
Tropfenform beidseitig, Chyu et al. (1997) 
Halbkreisform beidseitig, Chyu et al. (1997) 
Tropfen- Glatt, Chyu et al. (1997) 
Halbkreis- Glatt, Chyu et al. (1997) 

Pins versetzt, Chyu et al. (1998) 

Pins gleichmäßig, Chyu et al. (1998) 
Diamanten versetzt, Chyu et al. (1998) 
Würfel versetzt, Chyu et al. (1998) 
Diamanten gleichmäßig, Chyu et al. (1998) 
Würfel gleichmäßig, Chyu et al. (1998) 
Glatt, Mahmood et al. (2001) 


Dellen-Erhöhungen gegenüber versetzt, Mahmood et al. (2001) 
Dellen-Erhöhungen gegenüber ausgerichtet, Mahmood et al. (2001) 


Delle konkav klein, S. W. Moon et al. (2002) 
Delle konkav groß, S. W. Moon et al. (2002) 
Delle zyl. flach, S. W. Moon et al. (2002) 
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Op, 


Delle zyl. tief, S. W. Moon et al. (2002) 

Rippen orthogonal, Taslim et al. (1994) 

Rippen 45°, Taslim et al. (1994) 

V-Rippe zugewandt, Taslim et al. (1994) 
V-Rippe abgewandt, Taslim et al. (1994) 

Tang. Drallwinkel ca. 60°, Khalatov et al. (2000) 
2. Pass. mit 2. Drallerzeuger, Khalatov et al. (2000) 
Axiales Leitgitter 30°, Anmadvand et al. (2010) 
Axiales Leitgitter 60°, Anmadvand et al. (2010) 
Tang. Drallerzeuger div., Khalatov et al. (2011) 
Zyl. Dellen, Xie et al. (2010) 

Zyl. Erhéhungen, Xie et al. (2010) 

Rippe voll, Salameh et al. (2012) 

Rippe beschn. 46, Salameh et al. (2012) 
Rippengitter Ro = 0, Acharya et al. (2005) 
Rippengitter Ro = 0,1, Acharya et al. (2005) 
hohe V-Rippe einseitig, Maurer et al. (2007) 
hohe V-Rippe beidseitig, Maurer et al. (2007) 
Versch. Rippen, Cho et al. (2001) 


Abbildung 2.5: Übersicht der Literaturdaten zur Erhöhung des internen Wärmeübergangs 


Verfahren wie auch Drallkühlverfahren evaluiert. Zusätzlich wird der Effekt der Rotation bei 
Einsatz in Laufschaufeln thematisiert. Die Autoren betonen die hohe Komplexität sowie die 
notwendige Berücksichtigung von vielen Einflussfaktoren. Diese Wechselwirkungen werden in 
den bisherigen Untersuchungen nicht ausreichend berücksichtigt. Des Weiteren ist es von hoher 
Wichtigkeit, nicht nur globale Parameter zu berücksichtigen, sondern lokal aufgelöste Tempe- 
raturfelder zu bestimmen, da lokale Temperaturgradienten die Lebensdauer von Bauteilen stark 
beeinflussen. 

Der Effekt der Rotation wird generell als gering bewertet. Je komplexer die Innengeometrie, de- 
sto mehr überwiegt der Einfluss der Reynoldszahl. Die Rotationszahl hat wenig Einfluss auf die 
gesamte thermische Leistungsfähigkeit. Einzig eine Verschiebung der Nußeltzahlen zwischen 
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Saug- und Druckseite ist bei hohen Rotationszahlen in Hinblick auf lokale Temperaturvertei- 
lungen zu berücksichtigen. 

Eine umfangreiche Zusammenstellung zu den aerodynamischen und thermodynamischen Ef- 
fekten von Drallströmungen ist in Mitrofanova (2003) zu finden. Neben den diskutierten grund- 
legenden Effekten sowie den Vor- und Nachteilen verschiedener Geometrien wird außerdem 
ein Überblick über die Herausforderungen bei der numerischen Berechnung von diesen Strö- 
mungen gegeben. Mitrofanova betont ebenfalls nochmals die Wichtigkeit von experimentellen 
Untersuchungen. 

Auch Han et al. (2010) fassen den Kenntnisstand zu internen Kühlverfahren mit und ohne Rota- 
tion zusammen. Sie thematisieren die Komplexität des Wärmeübergangsproblems, speziell an 
der Schaufelvorderkante. Han et al. sehen die Weiter- und Neuentwicklung von Kühlverfahren 
als wichtiges Forschungsthema, das neue Triebwerkskonzepte ermöglichen kann. 


2.3.2 Kombinierte Untersuchungen zur Gesamtkühleffektivität 


Die bisher gezeigten Studien beschränken sich auf die Bewertung der einzelnen Teilsysteme 
Innenkühlung und Filmkühlung. Zwar wird die Wechselwirkung der Teilsysteme thematisiert, 
Untersuchungen von vollständigen Kühlsystemen wurden aber noch nicht vorgestellt. Die für 
eine solche Beurteilung verwendete Zielgröße ist meist die Totalkühleffektivität TCE 


Ta — Ts(X) 


TCE= 
Tu-Tk ` 


(2.6) 
welche die im Allgemeinen ortsabhängige Oberflächentemperatur Ts unter Verwendung der 
Freistromtemperatur des Heißgases Ty sowie einer Referenztemperatur der Kühlluft Tkg entdi- 
mensioniert. Eine perfekte Kühleffektivität TCE = 1 bedeutet das Erreichen der Kühllufttempe- 
ratur, ein Wert von TCE = 0 das Erreichen der Heißgastemperatur und damit keine Kühlwirkung 
der Oberfläche. Bei gegeben Randbedingungen in modernen Maschinen und verwendbaren Ma- 
terialien sollte diese Kühleffektivität über TCE Z 45% liegen (erste Rotorstufe, Rolls-Royce 
Deutschland, 2017) 


H. K. Moon et al. (2000) analysieren diese Kühleffektivität an einer realistischen (orig. engi- 
ne hardware) Leitschaufelreihe mittels Infrarotthermographie. Aufgrund rein konvektiver Küh- 
lung unter Einhaltung von Reynolds- und Machzahlähnlichkeit (global an Ein- und Auslass, 
über lokale Schaufeldruckverteilungen wird keine Aussage getroffen) ermitteln sie das Tem- 
peraturverhältnis I” zwischen Heißgas und Kühlluft als wichtigste Einflussgröße. Messdaten 
können nur für die Saugseite erfasst werden: Die Druckseite ist nicht sichtbar, die Sichtwinkel 
auf die Vorderkante sind zu hoch, um eine genaue Messung zu garantieren. Aufgrund hoher 
Unsicherheiten bei Messungen mit hohen Freistromtemperaturen (Oberflächentemperatur bis 
T, = 225°C) verwenden sie ein reversibles Verfahren, bei der die Innenströmung erhitzt wird 
und so parasitäre Strahlung von heißen Kanalwänden auf die Messfläche vermieden wird. Durch 
diese Einschränkungen kann der Messfehler der Totalkühleffektivität auf +10 % reduziert wer- 
den. 

Nathan et al. (2014) untersuchen die Vorderkante einer Leitschaufel und beschreiben detail- 
liert die für eine Übertragbarkeit der Daten einzuhaltenden Parameter. Diese sind, neben der 
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aerodynamischen Ähnlichkeit, die Biotzahl gebildet mit der Wärmeleitfähigkeit und Geometrie 
des Prüflings, das Verhältnis von internen zu externen Wärmeübergangskoeffizienten und der 
Aufwärmfaktor® (wobei die letzten beiden Ähnlichkeiten bei korrektem Temperaturverhältnis 
T zwischen Ty und Tx und aerodynamischer Ähnlichkeit implizit gegeben sind). Unter Ein- 
haltung dieser Ähnlichkeiten verwenden Nathan et al. stark hochskalierte Geometrien in einer 
linearen Kaskade. 

Ihre Ergebnisse für die Vorderkante zeigen Totalkühleffektivitäten zwischen TCE ~ 0,4 (kei- 
ne Filmkühlung) und TCE = 0,7 (starke Filmkühlung). Zusätzlich können sie durch adiabate 
Messungen auch die adiabate Filmkühleffektivität isolieren. Durch den Vergleich der beiden 
Messgrößen zeigen sie, dass die Wärmeleitung im Material sowie die interne Kühlung dazu 
führen, dass die Totalkühleffektivität höher und homogener verteilt ist als die adiabate Film- 
kühleffektivität. 

Auch in diesen aktuellen Untersuchungen werden wieder Kompromisse eingegangen: Der geo- 
metrische Messbereich ist verhältnismäßig klein, das Dichteverhältnis zwischen dem Heißgas 
und der Kühlluft kann aufgrund technischer Umstände nicht eingehalten werden. Außerdem 
wird keine Aussage über die erreichte Machzahlähnlichkeit getroffen. 

Kirollos et al. (2016) untersuchen unter Einhaltung der Reynoldszahl, Machzahl und Biotzahl 
die Oberflächentemperatur und Kühleffektivität einer unskalierten Leitschaufel (Innen- und 
Filmkühlung) unter Verwendung von Infrarotthermographie. Die Einhaltung des Temperatur- 
verhältnisses T zwischen Kühlluft und Heißgas wird nicht realisiert. Uber das resultierende 
Dichteverhältnis der Filmkühlung wird keine Aussage getroffen. Ziel der Arbeit bleibt ein re- 
lativer Vergleich zwischen zwei Herstellungsverfahren. Die Übertragbarkeit der Daten muss 
in Frage gestellt werden. Die thermographische Messung der Druckseite wird aus optisch un- 
günstiger Position mit hoher Objektstreckung in Tiefenrichtung und unter hohem Sichtwinkel 
durchgeführt. Eventuelle Einflüsse werden vernachlässigt. 

Ramachandran et al. (2015) betonen nochmals die Wichtigkeit der Einhaltung der vollen aero- 
thermischen Ähnlichkeit sowie die Notwendigkeit der Analyse von kompletten Kühlsystemen. 
Diese Feststellungen stützen sie nur mit numerischen Daten, die anhand von vereinfachten Geo- 
metrien erzeugt wurden. Ohne die Einhaltung dieser Ähnlichkeit sei keine genaue Prognose 
möglich. 


6 Die dimensionslose, graduelle Erhöhung der Kühllufttemperatur durch kontinuierliche Wärmeaufnahme im 
Kühlluftkanal. 
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2.4 Bestimmung der Oberflächentemperaturen mittels Infrarotthermo- 
graphie 


Die Infrarotthermographie eignet sich, um lokal und physikalisch hoch aufgelöste Strahlungs- 
bilder von Oberflächen aufzunehmen. Aufgrund komplexer Strahlungssituationen in heißen 
Prüfständen kann nur durch eine extensive Kalibrierung hieraus ein Temperaturbild mit ge- 
ringem Messfehler abgeleitet werden. Im Folgenden sollen die nötigen Grundlagen kurz zu- 
sammengefasst werden. Des Weiteren wird ein Überblick über die bisherigen Errungenschaften 
des Instituts für Thermische Strömungsmaschinen gegeben, die diese hoch genaue Temperatur- 
erfassung ermöglichen. 


Jeder Körper, dessen thermodynamische Temperatur von OK abweicht, sendet Strahlung aus. 
Als Infrarot wird hierbei der Wellenlängenbereich von ca. A = 0,78um bis ca. A = 1.000um 
bezeichnet. Die ideale Ausstrahlung eines Körpers wird durch das Plancksche Strahlungsgesetz 


hc? 
beschrieben. Hierbei repräsentieren h,c und k die physikalischen Konstanten Plancksches Wir- 
kungsquantum, Lichtgeschwindigkeit sowie die Boltzmann-Konstante. Die maximale Ausstrah- 
lung ist sowohl wellenlängen- als auch temperaturabhängig. Die Lage dieses Maximums wird 
durch das Wiensche Verschiebungsgesetz vereinfacht mit der Wienschen Konstante b 


b _2,898-10""m/K 
Ts Ts 
beschrieben. Hieraus ist zu erkennen, dass bei gegebenen Einsatztemperaturen Strahlungsmaxi- 
ma im mittleren Infrarotbereich von ca. A = 3 um bis ca. A = 8 um erwartet werden. Eingesetzt 
wird daher das Kamerasystem FLIR SC6000 mit gekühltem Indium-Antimonid-Detektor, der 
im Bereich von A = 3um bis A = 5 um sensitiv ist. Der Detektor ist als FPA mit 512 - 640 Pi- 
xeln ausgeführt. Mit dem verwendeten optischen Aufbau (Objektiv f = 50mm, Zwischenringe) 
wird so eine Auflösung von ca. 6 px/mm bis 10 px/mm erreicht. Diese ist proportional zum Ab- 
bildungsmaßstab und damit bei fester Brennweite auch direkt abhängig vom Objektabstand. 
Werden nun im Experiment reale Oberflächen betrachtet, ergibt sich die Ausstrahlung einer 
solchen Oberfläche mit der Emissivität e zu 


Wx(A,Ts) = (2.7) 


AW,max = (2.8) 


Ws(A, 0,Ts) = €(A,0,Ts) -Wex(A, Ts). (2.9) 


Jegliche Hindernisse schwächen diese Ausstrahlung ab, bevor sie den Detektor der Kamera 
erreicht. Hier sind sowohl die Atmosphäre wie auch alle Einbauten im optischen System (Fens- 
ter, Linsen und Filter) zu berücksichtigen. Diese können vereinfacht in der totalen Transmission 
Tor zusammengefasst werden, die sich als Produkt aller Einzeltransmissionen 7; ergibt. Die auf 
dem Detektor einfallende, von der betrachteten Oberfläche ausgesendete Strahlung, Es ergibt 
sich dann zu 

Es(A,0, Ts, Tj) œ tor(A, Ti) “Ws (A, 0, Ts). (2.10) 


Dieser Strahlung werden noch weitere Anteile tiberlagert. Dies sind Eigenstrahlungen von op- 
tischen Komponenten, Reflexionen, die an optischen Grenzflächen entstehen und, als größter 
Anteil, die von der zu untersuchenden Oberfläche reflektierte Strahlung. Letztere, sogenannte 
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Offsetstrahlung Eor, ist vor allem bei geschlossenen Umgebungen mit heißen Wänden maßgeb- 
lich. Diese Situation, die bei der Untersuchung von gekühlten Bauteilen in heißer Umgebung 
zwingend vorliegt, erschwert die Kalibrierung und kann als Hauptgrund der kritisierten Ska- 
lierungen der in Unterabschnitt 2.3.2 präsentierten Untersuchungen gesehen werden. Die am 
Detektor einfallende Strahlung kann dann durch die Überlagerung von Oberflächenstrahlung 
Es und Offsetstrahlung Eoff modelliert werden: 


Ep(A, 0, Ts, T;) = Es(A,0,T3,T;)+ Eotr(A, Ti). (2.11) 


Die komplexe Strahlungssituation kann im Allgemeinen nicht a priori bestimmt werden. Daher 
wird eine in-situ Kalibrierung notwendig. Martiny et al. (1996) schlagen eine semi-empirische, 
skalare Kalibrierfunktion vor, die dem Planckschen Gesetz nach Gleichung 2.7 nachempfunden 
ist (vereinfacht ohne Argumente dargestellt): 


b 
Up bzw. Ts = 


F 
elf In(’/Up + f) 
Hierbei beschreibt Up die Antwort des Detektors auf eine einfallende Strahlung Ep. Schulz 
(2000) zeigt die Anwendbarkeit zur Temperaturmessung in geschlossenen Kanälen. Die drei 
Parameter r, b, f müssen mittels Stiitzstellen, meistens in Form von Thermoelementen, während 
des Versuches bestimmt werden. Die in der Funktion beschriebene Antwort Up trennt noch 
nicht explizit zwischen Oberflächen- und Offsetstrahlung. 

Zur Kalibrierung müssen folglich im Versuchsaufbau mindestens drei Stützstellen vorhanden 
sein. Die Einbringung dieser ist konstruktiv aufwendig. Des Weiteren hat die Wahl der Stützstel- 
len einen hohen Einfluss auf den Restfehler. Gerade bei hohem Temperaturumfang des Mess- 
objektes kann keine optimale Kalibrierfunktion gefunden werden. 


(2.12) 


Auf Basis des linearen Zusammenhangs in Gleichung 2.11 präsentieren Ochs et al. (2009, 2010) 
eine weiter verbesserte Methodik der Bildverarbeitung und Kalibrierung, die genau diese Pro- 
blematik adressiert. Sie trennen die zu kalibrierende Strahlung in Gleichung 2.12 explizit in 
Objekt- und Offsetstrahlung und erweitern die Funktion zunächst auf vier Unbekannte, die es 
zu bestimmen gilt, wobei / das linearisierte, skalierte Antwortsignal des Detektors beschreibt: 


a ae 
(+) 


Ip — Loge 


Ts = (2.13) 
Von diesen vier Unbekannten können nun 7,b, f, die die optische Konfiguration beschreiben, 
vorher unter vereinfachten Bedingungen bestimmt werden, wenn zwei wichtige Forderungen 
erfüllt werden: 


e Die Temperaturabhängigkeit der Transmission der Komponenten i im optischen Pfad 
(4:/or), ist vernachlässigbar klein. Dies gilt in guter Näherung für Luft bei gegebenen 
Temperaturen. Außerdem ist das optische System der Kamera nahezu isotherm. Einzig 
die Transmission der verwendeten Saphirscheiben zeigt eine starke Temperaturabhän- 
gigkeit. Ochs et al. zeigen nun, dass durch Verwendung eines Tiefpassfilters mit einer 
Grenzwellenlänge von After = 4,1 um diese Temperaturabhängigkeit stark reduziert wer- 
den kann. Dadurch wird die Forderung der Temperaturunabhängigkeit der Transmission 
erfüllt. 
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e Ein lineares Systemverhalten wird vorausgesetzt, die Strahlungsantwort des Detektors 
muss linearisiert werden. Dies wird von Ochs et al. durch eine pixelweise, polynomielle 
Korrektur (Grad 2) erreicht. Die Koeffizientenmatrix der Korrektur kann in einem Vor- 
versuch bestimmt werden. 


Zur Bestimmung von r, b und f genügt dann eine Oberfläche mit bekannter und regelbarer 
Temperatur zur Erzeugung der Stützstellen für Gleichung 2.13. Es müssen die gleiche Ober- 
flächenbeschaffenheit und der gleiche optische Aufbau (Scheiben, Objektiv, Kamera) wie im 
eigentlichen Versuch realisiert werden. In diesem muss nur noch der vierte Parameter, die Off- 
setstrahlung Jff, bestimmt werden. Es genügt eine Stützstelle. Bei mehreren Stützstellen kann 
eine Fehlerminimierung durchgeführt werden. 


Nach Linearisierung der Strahlungsantwort kann das Detektorsignal außerdem ohne Einschrän- 
kung durch die bekannte Belichtungszeit’ geteilt werden. Es entsteht die der Strahlungsleistung 
proportionale Größe /. Die Belichtungszeit kann als zusätzlicher Freiheitsgrad während der Da- 
tenerfassung verwendet werden, um eine optimale Detektorausnutzung und damit minimalen 
Messfehler zu erreichen. Darauf basierend entwickeln Ochs et al. zwei neue Verfahren, die 
maßgeblich zur Einsatzfähigkeit der Infrarotthermographie bei modernen Filmkühlexperimen- 
ten beitragen: 


e Eine belichtungszeitunabhängige Ungleichförmigkeitskorrektur (non-uniformity correc- 
tion, NUC). Durch Einsatz dieser kann über einen weiten Bereich von Oberflächentem- 
peraturen sichergestellt werden, dass eine homogen strahlende Fläche ein homogenes 
Sensorsignal erzeugt. Fehler, die durch Linsenfehler (speziell an den Bildrändern) so- 
wie Abschattungen und Absorption im optischen Strahlengang des Objektivs entstehen, 
werden hierdurch korrigiert. 


e Durch Variation der Belichtungszeit und Aufnahme mehrerer Bilder der gleichen Sze- 
ne kann der Dynamikumfang des Detektors virtuell erweitert werden (sog. high dynamic 
range, HDR Aufnahmen). Nach Linearitäts- und Ungleichförmigkeitskorrektur können 
aus jeder Aufnahme Bildbereiche identifiziert werden, die optimale Belichtungen zeigen. 
Hieraus wird anschließend das fertige Bild zusammengesetzt. Aufgrund der belichtungs- 
zeitunabhängigen Form der Daten kann danach die in-situ Kalibrierung wie zuvor global 
für das komplette Bild bestimmt werden. Dieses Verfahren ist gerade bei filmgekühlten 
Bauteilen mit großer Spreizung der Oberflächentemperaturen unabdingbar. 


Durch diese am Institut für Thermische Strömungsmaschinen entwickelten Verfahren können 
selbst bei komplexen Strahlungssituationen und hohen Umgebungstemperaturen durch Verwen- 
dung der Infrarotthermographie hoch genau und hoch aufgelöst Oberflächentemperaturen be- 
stimmt werden. Eine wichtige Voraussetzung bei der beschriebenen Kalibrierung ist eine hohe 
und diffus-homogene Emissivität € Æ f(A,9,7s) der zu untersuchenden Oberfläche. Ersteres 


7 Die verwendete Kamera ist ohne mechanischen Verschluss aufgebaut. Vom Hersteller wird anstelle der Belich- 
tungszeit der Begriff der Integrationszeit (integration time) verwendet. Aufgrund der Konsistenz zu anderen 
Arbeiten und kamerabasierten Messtechniken wird hier der Begriff der Belichtungszeit verwendet. Beide Be- 
griffe beschreiben die selbe physikalische Zeitspanne. 
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führt zu einem reduzierten Anteil der Offsetstrahlung, letzteres stellt sicher, dass die an weni- 
gen Punkten bestimmte in-situ Kalibrierung für große Bildbereiche (bzw. das komplette Bild) 
gilt. Für die bisher untersuchten Konfigurationen konnte dies durch die Beschichtung der Ober- 
flächen mit Nextel Velvet Coating erreicht werden. Untersuchungen von Lohrengel et al. (1996) 
zeigen die Konstanz der Emissivität dieser Beschichtung im betrachteten Wellenlängen- und 
Temperaturbereich. Neuere Studien von Dury et al. (2007) bestätigen diese Ergebnisse über 
einen breiteren Wellenlängenbereich. Außerdem zeigen sie, dass die geringe Restreflektivität 
praktisch vollkommen diffus ist. Durch thermische Alterung der Beschichtung bei mehreren 
Tagen über ~ 200°C stellen sie einen geringen Anstieg der diffusen Reflektivität (also einen 
Abfall der Emissivität) fest. Da die Einsatztemperatur der Prüflinge im Rahmen dieser Arbeit 
in genanntem Bereich liegt, werden alle Kalibrier- und Prüfoberflächen bei 230°C vor dem 
eigentlichen Einsatz gealtert. Die geringfügig geminderte Emissivität wird zugunsten einer ho- 
mogenen und reproduzierbaren Oberfläche in Kauf genommen. 


Bisher nicht berücksichtigt wurde die Sichtwinkelabhängigkeit der Emissivität. Dies war, wie 
in Unterabschnitt 2.3.2 dargestellt, aufgrund geringer Sichtwinkel nicht notwendig, wurde ver- 
nachlässigt oder Datenbereiche wurden nicht ausgewertet. Bei komplexeren Oberflächen, die 
sich durch eine Erstreckung in alle drei Raumrichtungen auszeichnen, muss diese Winkelab- 
hängigkeit berücksichtigt werden. In Kapitel 4 und Kapitel 5 werden Verfahren vorgestellt, die 
diese Sichtwinkelabhängigkeit berücksichtigen und die thermographische Messung auf beliebig 
geformten Oberflächen mit hoher Genauigkeit ermöglichen. 


Alternative Messtechniken 


In direkter Konkurrenz zur Infrarotthermographie stehen Verfahren, die eine lokal hoch aufge- 
löste Bestimmung der Oberflächentemperatur erlauben. Zwar können durch den Einsatz von 
mehreren Medien und den dabei auftretenden lokalen Konzentrations- oder Partialdruckunter- 
schieden Filmkühl- und Wärmeübergangsgrößen abgeleitet werden, eine genaue Bestimmung 
der Oberflächentemperatur ist aber nicht möglich. 

Als Alternative verbleiben Beschichtungen, die optisch auf ihre Temperatur reagieren. Hier 
muss zwischen Flüssigkristallen unterschieden werden, die einen Farbumschlag beim erreichen 
einer bestimmten Temperatur erzeugen, und thermosensitiven Farben (thermo sensisitve paints, 
TSP), die bei Anregung Licht einer anderen Wellenlänge aussenden. 

Umschlagsfarben für hohe Temperaturen über ~ 1.100K werden bei Maschinentests einge- 
setzt, um Konzepte zu validieren. Die Farbe zeigt nach dem Test nur die maximal erreichte 
Temperatur. Die Auswertung erfolgt qualitativ anhand eines Farbvergleichs. Die erreichte 
quantitative Genauigkeit ist gering, Fehler im Bereich von 100K werden beobachtet (Rolls- 
Royce Deutschland, 2017). 

Gerade im wissenschaftlichen Bereich sind TSP erfolgversprechend. Russin et al. (2009) 
benutzen diese Farben zur Bestimmung von Oberflächentemperaturen. Die Kalibrierung erfolgt 
außerhalb der Messung bei Umgebungsdruck, was einen negativen Einfluss auf die erreichte 
Genauigkeit hat. Jordan et al. (2011) verwenden ein ähnliches Kalibrierverfahren außerhalb 
des eigentlichen Versuchs. Lorenz et al. (2013) wenden erstmals ein in-situ Verfahren zur 
Kalibrierung der TSP an. Sie erreichen so eine hohe Genauigkeit und Auflösung. Die Autoren 
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vergleichen die Ergebnisse mit thermographischen Messungen und können eine vergleichbare 
Qualität erreichen. Ein ähnlicher Vergleich wird von Brack et al. (2016) durchgeführt. Auch sie 
zeigen eine gute Übereinstimmung der Ergebnisse der Thermographie und der TSP. Weiterhin 
weisen sie darauf hin, dass gerade bei instationären Phänomenen die Ansprechzeit der TSP zu 
einer Messabweichung führen kann. 

Pütz (2018) setzt ein auf der Methodik von Lorenz et al. basierendes Verfahren. Er bestimmt 
erfolgreich Temperaturen und leitet aus diesen Wärmeübergangskoeffizienten auf dreidimen- 
sional konturierten Seitenwänden mit glatter sowie rauer Oberfläche ab. 

Weiterhin werden thermosensitive Farben zur Bestimmung von Temperaturen auf internen 
Oberflächen eingesetzt. Durch den Aufbau der zu untersuchenden Kanäle aus transluzentem 
Material und Beschichtung der inneren Flächen mit TSP kann die Reaktion der Farbe von 
außen (also rückseitig) erfasst werden. So bestimmen beispielsweise Domaschke et al. (2012), 
Kunstmann et al. (2013), Waidmann et al. (2013) und Rao et al. (2017) Temperaturen an 
Innenflächen, die dann zur Ableitung von Wärmeübergangskoeffizienten verwendet werden. 
Die Vorteile dieser Methode liegen in der geringen Dicke der Beschichtung: So können die 
von Pütz eingesetzten rauen Oberflächen untersucht werden, die nach Beschichtung mit Nextel 
Velvet Coating für den Einsatz der Thermographie andere Rauheiten zeigen würden. Außerdem 
ist keine Beachtung der Offset-Strahlung notwendig. Aufgrund der Einbettung der aktiven 
Bestandteile in transparente Polymerbinder ist die Einsatztemperatur auf ca. 200°C begrenzt, 
die Infrarotthermographie kann bei deutlich höheren Temperaturen eingesetzt werden. Die hier 
verwendete Beschichtung Nextel Velvet Coating erlaubt den Einsatz bei bis zu 300 °C. Andere 
Beschichtungen werden auch bei weitaus höheren Temperaturen eingesetzt. So untersuchen 
Brandt et al. (2008) erfolgreich hoch emissive Beschichtungen bei bis zu 1.500K. 

Über eine Winkelabhängigkeit bei TSP, wie sie für die Infrarotthermographie Thema dieser 
Arbeit ist, liegen zum jetzigen Zeitpunkt keine Veröffentlichungen vor. 
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2.5 Zielsetzung der Arbeit 


Die präsentierten Studien zeigen die Vor- und Nachteile verschiedener interner Kühlsysteme. 
Da die Grenzen der Leistungsfähigkeit von konventionellen Systemen, zusammengefasst in 
Abschnitt 2.3, erreicht werden, müssen neue Methoden entwickelt werden. Vielversprechend 
ist das Konzept der Drallkühlung, dessen Vor- und Nachteile in Unterabschnitt 2.3.1 erläutert 
wurden. Die bisher durchgeführten Untersuchungen zeigen, dass dieses Konzept in einfachen 
Versuchen vielversprechende Erhöhungen der Nußeltzahl erzielen kann. Eine weitere Unter- 
suchung von modernen Kühlkonzepten unter der Verwendung von Drallkühlung ist daher aus 
technischer Sicht hochaktuell, wie auch aktuelle Patente führender Hersteller zu solchen Kon- 
figurationen zeigen (Rolls-Royce Deutschland et al., 2012). 


Bei der Beurteilung der Konzepte werden bisher immer Kompromisse eingegangen. Die ther- 
mische Leistungsfähigkeit wird meist an einfachen Kanalgeometrien beurteilt, die so nicht in 
realen Schaufelgeometrien zu verwirklichen sind. Einige Untersuchungen versuchen zwar den 
Einfluss der Wechselwirkung verschiedener Parameter zu quantifizieren, es Können aber meist 
nur binäre Zusammenhänge erfasst werden. Auch die Einhaltung von äußeren und inneren Ähn- 
lichkeitskenngrößen, speziell dem Verhältnis von äußerem zu innerem Wärmeübergang hu/nx 
und der Biotzahl, wird nicht gewährleistet. 

Weiterhin wird die Interaktion von Innenströmung und Filmkühlung nur wenig betrachtet. Ein 
Kühlsystem kann nur als Ganzes leistungsfähig sein. Eine komplexe Innenströmung muss auch 
im Hinblick auf die Einströmung in die Filmkühlbohrungen und die daraus resultierende Film- 
kühleffektivität optimiert sein. Auch hier gilt wieder, dass mit immer komplexeren Innenströ- 
mungen, gekennzeichnet durch eine nicht mehr hauptsächlich radiale Strömung in den Kanä- 
len, eine Zunahme der Interaktion erwartet wird. Die Vorhersage der finalen Bauteiltemperatur 
wird weiter erschwert. Numerische Vorhersagen des Wärmeübergangs, die nur an ausgewähl- 
ten Beispielen gezeigt wurden, bleiben ungenau. Aufgrund komplexer Geometrien muss mit 
geometrischen Vereinfachungen sowie lokaler Mittelung und, daraus resultierend, Turbulenz- 
modellen gearbeitet werden. Diese zeigen gerade im für den Wärmeübergang maßgeblichen 
Grenzschichtbereich Unzulänglichkeiten bei der korrekten Prognose der Strömungsgrößen. 
Eine Veränderung der Geometrie oder der Parameter während des Einsatzes, entweder bewusst 
durch nötige Anpassungen zur Implementierung in Schaufelprofilen und Maschinenprozessen 
oder unbewusst durch Abweichungen während der Fertigung oder im Einsatz, kann die vorher- 
gesagten Eigenschaften deutlich verändern. Diese Überlegungen gelten für bewährte Konzepte, 
sind aber umso wichtiger bei Drallkühlverfahren. Die hier erreichte Kühleffektivität hängt maß- 
geblich von der erreichten Drallzahl und der Stabilität der Drallströmung ab. Diese reagieren 
wiederum hoch sensitiv auf geometrische Veränderung. 

Aufgrund dieser Unsicherheiten ist es üblich, neue Kühlkonzepte in kosten- und zeitaufwen- 
digen Tests in der Einsatzumgebung selbst zu validieren. Die erreichte Messgenauigkeit und 
Auflösung der wichtigen Zielgrößen sind aber so gering, dass gezielte Verbesserungen schwer 
auf Basis dieser Daten durchzuführen sind. Eine Möglichkeit, moderne Kühlkonzepte in ihrer 
fertigen und gesamten Funktion vor einem Einsatz in der Maschine mit hoher Auflösung und 
Genauigkeit zu untersuchen und weiter zu optimieren, existiert bisher nicht. Versuche, solche 
Messungen durchzuführen, sind in Unterabschnitt 2.3.2 zusammengefasst. Um die Untersu- 
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chungen zu ermöglichen müssen wieder Kompromisse eingegangen werden. Dies aber macht 
die Übertragung der in den Studien gewonnenen Daten unsicher, was gerade die Entwicklung 
von komplett neuen Konzepten (wie z.B. Drallkühlung) verhindern kann, da Hersteller ein 
großes wirtschaftliches Risiko eingehen müssen. 


Diese Arbeit stellt nun Verfahren zur Verfügung, die genau diese Beurteilung von Kühlkonzep- 
ten bei realistischen Randbedingungen ermöglicht. Dafür wird zunächst in Kapitel 3 ein neuer 
Prüfstand vorgestellt, der Untersuchungen an einer skalierten Schaufelgeometrie ermöglicht. 
Durch weitestgehende Einhaltung aller wichtigen Kennzahlen können anhand der Ergebnisse 
direkt Rückschlüsse auf die Leistungsfähigkeit in der Maschine gezogen werden. 

Die wichtigste Zielgröße, die lokal aufgelöste Oberflächentemperatur, wird mithilfe der 
Infrarotthermographie bestimmt. Zum Erreichen der geforderten weitgehenden Ähnlichkeit 
müssen suboptimale Bedingungen für die Messtechnik in Kauf genommen werden. So müssen 
Oberflächentemperaturen auf stark gekrümmten und gekühlten Oberflächen erfasst werden. 
Die Anforderungen an die lokale sowie physikalische Auflösung sind hoch. Die Prüfkörper 
sind zudem in einem komplex geformten Kanal eingesetzt. Die optische Zugänglichkeit wurde 
für diese Untersuchungen erstmals der geometrischen Ähnlichkeit untergeordnet. Dies führt 
zu hohen Sichtwinkeln auf den Prüfkörper, die sowohl das emissive als auch das reflektive 
Verhalten der Oberfläche beeinflussen. Eine genaue Erfassung der Oberflächentemperatur 
erfordert demnach neue Korrektur- sowie Kalibrierverfahren, die in Kapitel 4 und Kapitel 5 
vorgestellt werden. Diese Verfahren ermöglichen erstmals die hoch genaue Erfassung von 
Oberflächentemperaturen auf nahezu beliebig geformten Bauteilen. 

Die mit diesen neuen Verfahren gewonnenen Ergebnisse werden in Kapitel 6 vorgestellt und 
diskutiert. Neben der Beurteilung der Kühleffektivität einer bewährten Geometrie werden neue 
Kühlkonzepte relativ zu dieser und absolut beurteilt. Durch die hohe Genauigkeit können 
Designempfehlungen formuliert werden, die neue Kühlkonzepte leistungsfähiger machen, 
bevor aufwendige Maschinentests durchgeführt werden. 


3 Entwicklung und Inbetriebnahme eines neuen 
Turbinenprufstandes 


Zur Untersuchung der Kühleffektivität von Rotorschaufeln mit den in Abschnitt 2.5 formulier- 
ten Zielen musste ein neuer Priifstand konzipiert und aufgebaut werden. Die wichtigste Anfor- 
derung an die Auslegung ist die Einhaltung aller fiir den Warmeiibergang wichtigen Ahnlich- 
keitskennzahlen. Wie in Kapitel 2 beschrieben sind dies die Machzahl May, die Reynoldszahl 
Rey, die Nußeltzahl Nuy und die Biotzahl Bin für die Heißgasströmung (Index Ly) sowie 
die Reynoldszahl Rex, die Nußeltzahl Nux und die Biotzahl Bix für die Kühlluft (Index Lx). 
Um die lokale Verteilung dieser Größen im Versuch abzubilden, muss die Abströmung des 
in der Maschine stromauf liegenden Stators nachgebildet werden. Außerdem soll das korrekte 
Temperaturverhältnis von Kühlluft zu HeiBgas I” = 7/7, und letztlich das Verhältnis der Wär- 
meübergangskoeffizienten hu/nx eingehalten werden. Auf die Rotation der Domäne wird im 
experimentellen Aufbau verzichtet. Im Folgenden werden die Skalierung, die Auslegung, und 
darauf aufbauend, die Inbetriebnahme des Prüfstandes beschrieben. 


3.1 Skalierung und Konstruktion 


Die Entwicklung und Inbetriebnahme des Prüfstandes gliedert sich in mehrere Teilschritte: In 
einer Vorauslegung wurde, ausgehend von Betriebsdaten eines Referenztriebwerkes, zunächst 
eine 1D-Skalierung durchgeführt und so ein geeigneter geometrischer Skalierungsfaktor sowie 
daran gekoppelt die thermodynamischen Betriebsparameter festgelegt. Diese Betriebsparame- 
ter und deren Abweichung zu den Maschinendaten sind in Tabelle 3.1 in den Spalten Prüfstand 
und Abweichung 1D + zusammengefasst. Auf Basis dieser Skalierung und der Turbinengeo- 
metrie wurde dann ein Prüfstandsmodell abgeleitet. Im Rahmen einer Inbetriebnahme wurde 
die Machzahlverteilung des Profils an die bekannte Zielverteilung der Maschine angepasst. Die 
finalen, experimentell ermittelten Betriebsparameter und deren Abweichung zur 1D-Skalierung 
sind in den Spalten Experiment, beziehungsweise Abweichung Exp.<+ 1D von Tabelle 3.1 zu- 
sammengefasst. Abschließend wurden zusätzlich qualitative Strömungsvisualisierungen mithil- 
fe von Ölanstrichfarben durchgeführt. 


Bestimmung der Skalierungsfaktoren 


Als primäre Vergleichsgröße der 1D-Skalierung wurde die Abströmmachzahl May > des Heiß- 
gases gewählt, da diese betriebspunktunabhängig und für die hoch kompressible Strömung um 
das Schaufelprofil maßgeblich ist. Zwar wird der Wärmeübergang vor allem von der Reynolds- 
zahl beeinflusst, diese ist aber stark abhängig vom Betriebspunkt und den thermodynamischen 
Randbedingungen. Durch die Einhaltung der Abströmmachzahl sowie der lokalen Machzahl- 
verteilung wird die Kinematik der Strömung und damit die lokale Verteilung der Reynoldszahl 
festgelegt. So muss nur noch überprüft werden, ob die skalare Abströmreynoldszahl Rey > ver- 
gleichbar zu der der Maschine ist. 


36 Entwicklung und Inbetriebnahme eines neuen Turbinenpriifstandes 


Eine Betrachtung der Ahnlichkeit liefert nach Umstellen den geometrischen Skalierungsfaktor 
S in Abhängigkeit der thermodynamischen Verhältnisse xp = *tig/xeng zwischen Prüfstand und 
Maschine (jeweils statisch, an Punkt 2): 


u 105 
SMa = (rie R” TEP PR) . (3.1) 


Bei dieser Skalierung sind für die Prüfstandsgrößen wichtige Grenzen zu beachten. Der Mas- 
senstrom ist aufgrund der Anlagenkapazität begrenzt. Der statische Druck im Prüfstand muss 
aufgrund der komplexen Konstruktion begrenzt sein. Das Temperaturverhältnis zwischen Kühl- 
luft und Heißgas I’ = Tkwot.1/Ti 101 & 0,54 soll eingehalten werden und die Kühlluft soll nicht 
gekühlt werden. Daraus ergibt sich: 


Mrig 7 0,6kg/s, pr = 0,05...0,2, Prig > latm, und Tig > 480K. (3.2) 


Zusätzlich zur Machzahl soll außerdem die Reynoldszahl vergleichbar sein. Diese bestimmt 
den Wärmeübergang durch Einfluss auf die Nußeltzahl (bei näherungsweise konstanter Prandtl- 
zahl) und damit auch die Biotzahl. Analog zu Gleichung 3.1 kann ein Skalierungsgesetz für die 
Reynoldszahl abgeleitet werden: 


SRe = ñR: KR” TR ppl. (3.3) 


Für eine Reynoldszahlähnlichkeit ist also immer ein größerer geometrischer Skalierungsfak- 
tor notwendig, wenn das Druckskalierungsverhältnis pr die oben genannte Grenze nicht über- 
schreiten soll. Bei gegebenem, machzahlbasierten Skalierungsfaktor wird demnach eine im Ver- 
gleich zur Maschine geringere Reynoldszahl im Prüfstand erreicht. 

Auf Basis dieser Überlegungen wurde ein Skalierungsfaktor von S = 3 gewählt. Neben den 
günstigen aerothermischen Betriebsbedingungen führt die so erreichte Vergrößerung der Geo- 
metrie zu einer höheren Auflösung der Ergebnisse bei gleichem optischen System. Die erreichte 
Reynoldszahl liegt dann ~ 30 % unter der erwarteten MTO-Reynoldszahl des Triebwerkes, da- 
mit aber immer noch in einem triebwerkstypischen Bereich. 

Die Nußeltzahlen an verschiedenen Schaufelpositionen ergeben sich für die Maschine aus dem 
bekannten Wärmeübergangskoeffizienten und der temperaturkorrigierten Wärmeleitfähigkeit 
des Heißgases. Über eine Dittus-Bölter-Skalierung (Dittus et al., 1930) bei konstanter Prandtl- 
zahl (Nu œ Re®:8) werden die Nußeltzahlen sowie, mithilfe der Wärmeleitfähigkeit der Luft, die 
Wärmeübergangskoeffizienten des Priifstandes bestimmt. Hieraus ergeben sich letztlich unter 
Verwendung der jeweiligen Materialkoeffizienten und Wärmeübergangskoeffizienten die Biot- 
zahlen. 


Mit korrekter Skalierung der Geometrie und den Betriebsparametern des Heißgases werden 
die Betriebsparameter der Kühlluft implizit festgelegt. Durch Einhaltung des Temperaturver- 
hältnisses I” verbleiben als Freiheitsgrade der Eintrittstotaldruck px tot, der Kühlluft und der 
Kühlluftmassenstromanteil & = k/ry. Diese beiden Größen sind durch die Druckverluste der 
jeweiligen internen Geometrie gekoppelt. 
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Im Rahmen dieser Untersuchungen wurde der Massenstromanteil ¢ als vorgegebene Größe ver- 
wendet!, der sich einstellende Einlasstotaldruck PXK.tot,1 Variiert zwischen verschiedenen Konfi- 
gurationen. 

Anhand der so festgelegten thermodynamischen Parameter kann die Reynoldszahl Rex; und 
die Machzahl Max | am Eintritt bestimmt werden. Aufgrund der niedrigen Stromungsgeschwin- 
digkeiten ist hier vor allem die Reynoldszahl zu beachten. Die Ableitung der Nußelt- sowie der 
Biotzahl erfolgt wieder durch die Dittus-Bölter-Skalierung der Maschinendaten unter Berück- 
sichtigung der temperaturabhängigen Stoffgrößen. 


Tabelle 3.1: Skalierung zwischen Maschine und Prüfstand, Randbedingungen des Experiments 


Maschine" Prüfstand Abweichung 1D +- Experiment’ Abweichung Exp.+> 1D 


Heißgas 
May4.2 0,822 0,822 0,0% 0,83 1,2% 
Rey 2 942.335 709.581 -24,7% 593.754 -16,3% 
NUH,LE 2.717 2.165 -20,3% 1.877 -9,1% 
Bin LE 0,341 0,379 10,0% 0,329 -13,2% 
THtot.i - 559K - 579,4K 3,6% 
PHtot,1 - 150.968 Pa - 125.450 Pa -16,9% 
TH stat,2 z 495K = 497,6 K 0,5 % 
MH - 1,37 kg/s - 1,17 kg/s 14,6% 

Kühlluft 
Max 1 0,274 0,227 -17,2% 0,222 -2,2% 
Rex, 167.613 108.361 -35,4 % 94.754 -12,6 % 
Nuķk,LE 253,8 179,0 -29,4 % 160,8 -10,6% 
Bik LE 0,171 0,168 1,7% 0,151 -10,1% 
Tx tot,1 5 304,3 K = 306,4 K 0,7 % 
PK.tot,1 z 158.466 Pa - 141.039 Pa 11,0% 
¢ - 3,4 - 3,42 0,6% 

Sonstige 
T 0,543 0,548 0,9% 0,529 -3,5% 
hn/hk 1,988 2,256 13,4% 2,18 -3,4 % 


“ Thermodynamische Maschinendaten sind zur Skalierung bekannt, unterliegen aber der Geheim- 
haltung 
' Experimentelle Daten ermittelt mit der Referenzgeometrie der Innenkühlung beim Nennbetriebs- 


punkt 


' Die untersuchten Konzepte sollen in der ersten Laufreihe eines Triebwerkes eingesetzt werden. Hier ist, je nach 
Konfiguration, der Druck der Kühlluft teilweise so hoch, dass dieser verlustbehaftet gedrosselt werden muss, 
um eine sichere Filmkühlung zu gewährleisten (Rolls-Royce Deutschland, 2017). Ein hoher Eintrittstotaldruck 
kann, je nach Anwendungsfall, sogar positiv bewertet werden. Als Vergleichsgröße der Kühlung wird daher 
der Massenstrom verwendet. 
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Neben den Kennzahlen für Heißgas und Kühlluft konnten damit auch die weiteren, in Unterab- 
schnitt 2.3.2 eingeführten Verhältnisse I’ und hu/nx weitestgehend eingehalten werden. 

Die als Ziel gesetzte Reproduktion der betriebspunktunabhängigen Machzahl konnte erreicht 
werden. Die restlichen Kenngrößen der Maschine sind stark betriebspunktabhängig und von 
einem Hochlastbetriebspunkt abgeleitet. Eine Verschiebung zu niedrigeren Reynoldszahlen auf 
bis zu = 25% der ursprünglichen Werte (und damit eine Reduktion der Abweichung zwischen 
Prüfstand und Maschine, speziell der Reynoldszahl) bei Betrachtung von variierenden Betrieb- 
spunkten ist zu erwarten. 
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Prüfstandsaufbau 


Die 1D-Skalierung der Betriebsparameter ist nur eine notwendige Bedingung für die Ähnlich- 
keit des Versuchs. Soll die Schaufelumströmung auch lokal aufgelöst ähnlich sein, muss die 
Geometrie von der Maschine übernommen werden. So wurde auf Basis der Passagengeome- 
trie vom Laufgitter des Triebwerks die Versuchsgeometrie abgeleitet. Hierbei wurden die not- 
wendigen Zugänge für die Messtechnik so implementiert, dass ein Kompromiss aus minimaler 
Strömungsstörung und guter Sichtfeldabdeckung der Schaufel erreicht wurde. Die finale Kon- 
struktion wird in Abbildung 3.1 als Seitenansicht und als planarer Mittelschnitt gezeigt. 

Das Heißgas strömt von rechts nach links durch die Kaskade. Nach einem Beruhigungsgit- 
ter folgt die Messstelle der Anströmbedingungen Myn,ı. Ein folgendes, durch iterative CFD- 
Rechnungen ausgelegtes Strémungsgitter” erzeugt die Anströmung der folgenden Rotorschau- 
fel im Relativsystem. Eingestellt werden der Anströmwinkel ßı = f(r) und die Massenstrom- 
dichteverteilung A = f(r). Dieses Strömungsfeld wird nicht weiter überprüft: Bei korrekter 
Staupunktlage sowie korrekter Machzahlverteilung an der Schaufeloberfläche wird von einer 
korrekten Anströmung ausgegangen. Der Abströmwinkel sowie die saug- und druckseitigen 
Abströmquerschnitte können über verstellbare Seitenwände, sog. Tailboards, angepasst wer- 
den. Hierbei ist das saugseitige Tailboard in drei Freiheitsgraden beweglich. So lässt sich der 
Strömungsquerschnitt an der nabenseitigen und der äußeren Endwand getrennt einstellen. 

Das druckseitige Tailboard ist nur rotatorisch einstellbar. Da die Strömung vor allem durch die 
transsonische Durchströmung des saugseitigen Kanals bestimmt wird, ist diese Einstellmöglich- 
keit hier ausreichend. Die Austrittsbedingungen werden am Punkt My* 2 gemessen und unter 
Verwendung der Isentropenbeziehungen auf die Schaufelhinterkante Ly 2 bezogen. 

Die Kühlluft wird dem Prüfling am Schaufelfuß zugeführt. Die Betriebsparameter werden di- 
rekt am Eintritt in die Schaufel gemessen?. Nach Durchströmen der prüflingsspezifischen In- 
nengeometrie, die detailliert in Abschnitt 3.2 vorgestellt wird, wird ein Teil der Kühlluft an der 
Schaufelspitze entnommen (Messstellen Mx >;). Die für die Filmkühlung verwendete Luftmen- 
ge ergibt sich aus der Massenbilanz. 

Zur Anwendung der Infrarotthermographie wurden drei optische Zugänge vorgesehen. Zur Ver- 
wendung kommt geschliffenes Saphirglas in optischer Qualität. Durch den ersten Zugang (0z1) 
können die Vorderkante und der stromab folgende Teil der Saugseite beobachtet werden. Der 
zweite Zugang (0z2), der wegen des beweglichen Tailboards als Doppelscheibe ausgeführt wur- 
de, ermöglicht die Untersuchung der restlichen Saugseite bis zur Schaufelhinterkante. Diese 
zwei Sichtfelder überlappen theoretisch, die Sichtwinkel sind im Überlappungsbereich aber so 
hoch (> 85°), dass selbst nach Anwendung der später vorgestellten Korrekturverfahren keine 
sichere Datenerfassung in diesem Bereich möglich ist. Der dritte Zugang (0z3) ist in die äußere 
Endwand implementiert und erlaubt die Erfassung der Vorderkante (Überlappung mit (0z1)) 
und eines Teils der stromab gelegenen Druckseite. 


2 Die Auslegung des Strömungsgitters wurde von der Rolls-Royce Deutschland Ltd & Co KG durchgeführt. Das 
Strömungsgitter wurde im Rahmen eines Forschungsprojektes zur Verfügung gestellt. 

3 Die Messung der Eintritts- und Austrittsparameter der Kühlluft K,1; Ok.2 erfolgt in separaten Messadaptern, 
die bei einem Austausch des Prüflings erhalten bleiben. So kann die Messabweichung zwischen den Prüflingen 
verringert werden. 
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Abbildung 3.1: Darstellung der Prüfstandskonstruktion, Außenansicht und Schnitt 
Bilder in ähnlicher Form veröffentlicht in Elfner et al. (2017) 


Herstellung der Prüflinge 


Aufgrund der komplexen internen Geometrie können die Prüflinge nicht mit konventionellen 
Fertigungsverfahren hergestellt werden. Daher werden diese durch das selective laser mel- 
ting, SLM hergestellt. Die allgemeine Formhaltigkeit der Außenflächen wurde mittels opti- 
scher Scanverfahren überprüft und durch iterative Anpassung der Fertigungsparameter auf unter 
0,2 mm Oberflächenabweichung in Richtung der lokalen Normalen reduziert. 

Die Rauheit der Oberflächen ist vergleichbar mit der in Feingussverfahren erreichbaren Rau- 
heiten. Die Außenflächen mussten zur besseren Beschichtung zuerst geglättet werden, bevor 
die für die Infrarotthermographie notwendige hoch emissive Beschichtung aufgetragen wurde’. 
Die Oberflächen der Innengeometrie können nicht nachbearbeitet werden. Sie zeigen nach der 
Fertigung und im Einsatz eine vergleichbare Rauheit wie die Außenflächen. 

Eine Herausforderung für dieses neue Fertigungsverfahren bleiben feine sowie, relativ zur Auf- 
baurichtung gesehen, überhängende Strukturen. An dieser Stelle sind vor allem die Filmkühl- 


* Nach Fertigung liegt die Rauheit bei R, = 30...45um. Durch die Nachbehandlung wird eine Oberflächenrau- 
heit von R, = 1...5um erreicht. Die Rauheit der Beschichtung ist nicht bekannt, ist aber qualitativ mit der 
Fertigungsrauheit vergleichbar. 
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Abbildung 3.2: Interne Geometrien der vorliegenden Untersuchung 


bohrungen zu nennen. Die Genauigkeit des SLM-Prozesses ist hier noch nicht ausreichend, 
die Abweichungen zwischen einzelnen Bohrungen sind zu groß. Dieses Verhalten wurde auch 
von Gritsch et al. (2013) und Kirollos et al. (2016) beobachtet. Um reproduzierbare Ergebnis- 
se zu erzielen, werden diese Bohrungen daher kleiner produziert und darauffolgend auf den 
gewünschten Nenndurchmesser aufgebohrt. So kann die Vergleichbarkeit zwischen verschiede- 
nen Prüflingen sichergestellt werden. 

Aus Gründen der Fertigung und der gewünschten Biotzahl-Ähnlichkeit wurde TiAl6V4 als 
Werkstoff ausgewählt. Das im Bauteil verschmolzene Material zeigt im Vergleich zum Strang- 
material ähnliche physikalische Eigenschaften. Für den Versuchsaufbau nennenswert ist die 
thermische Leitfähigkeit As, die bei As ~ 6,5 W/(mK) liegt (intern überprüft) sowie die Isen- 
tropie derselben. 


3.2 Interne Kühlung 


Die Außengeometrie und der Betriebspunkt der Außenströmung werden in den durchgeführten 
Versuchen nicht variiert. Die untersuchten Konzepte unterscheiden sich nur anhand der 
Innengeometrie und den Betriebspunkten der Kühlluft. Diese Geometrien und die Variation der 
Betriebspunkte werden im Folgenden vorgestellt. 

Im Rahmen dieser Arbeit werden neue, auf Drallkühlung basierende Verfahren untersucht. 
Zur besseren Einordnung der Ergebnisse wird zusätzlich eine konventionelle Innengeometrie 
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Abbildung 3.3: Geometrische Versionen zur Drallerzeugung 


evaluiert. Beide Konfigurationen sind in Abbildung 3.2 dargestellt. Links (a) ist jeweils ein 
repräsentativer Mittelschnitt gezeigt, rechts (b) der schematische Strömungskanal mit Wand- 
elementen zur Wärmeübergangserhöhung sowie den Filmkühlreihen. Ausgehend von der kon- 
ventionellen Geometrie wurde eine Drallgeometrie abgeleitet. Die beiden Vorderkantenkanäle 
mit kleinem Querschnitt wurden zu einem Kanal mit größerem Querschnitt zusammengefasst. 
Die Verrippung der Innenwand wurde entfernt. Am Fuß dieses neuen Kanals wird der Draller- 
zeuger eingesetzt. In den untersuchten Modellen wird nur der Drallerzeuger variiert, mögliche 


Versionen 


zur Drallerzeugung sind in Abbildung 3.3 gezeigt. Verwendet werden axiale 


Leitgitter, dargestellt in Abbildung 3.3(c) - Abbildung 3.3(e). Der Abströmwinkel ist mit ~ 65° 
nahezu konstant zwischen den Konfigurationen. Es kommen folgende Variationen zum Einsatz: 


GVI: 


NSG: 


OSGi: 


Dem generischen Leitgitter nachempfundene Konstruktion. Sechs Leitschaufeln 
mit kurzer Sehnenlänge, daher starke Umlenkung und leichte Falschanströmung 
bei rein axialer Kühlluftströmung. Geringe Versperrung durch die Nabe. 


Gitter aus drei stark verwundenen Leitschaufeln. Große Sehnenlänge, moderate 
Umlenkung. In Strömungsrichtung zunehmende, am Austritt starke Versperrung 
durch die Nabe. 


Aufbau mit zwei überlappenden Leitapparaten mit großer Sehnenlänge. Keine Na- 
be. Dieser Aufbau bietet zwei wichtige Vorteile: Eine Fertigung im Feingussverfah- 
ren ist prinzipiell möglich, dadurch wird eine Serienfertigung ermöglicht. Außer- 
dem ist eine parametrische Anpassung der Geometrie einfach möglich. So wurden 
neun Variationen vorausgelegt und fünf dieser Variationen untersucht. Drei davon 
werden im Detail analysiert (Nummer 0, 3, 8). Hierbei nimmt mit steigender Num- 
mer die Überlappung der beiden Leitapparate und damit die Versperrung des Strö- 
mungskanals zu. So wird eine Steigerung der Drallzahl bei gleichzeitiger Erhöhung 
des Eintrittstotaldrucks px tor, ı erreicht. 
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Abbildung 3.4: Verteilung der isentropen Machzahl bei 50 % Schaufelhöhe 


Die jeweiligen Gitter wurden numerisch (durchgeführt von der Rolls-Royce Deutschland Ltd 
& Co KG) vorausgelegt. Eine weitere, experimentelle Überprüfung der erreichten Drallzahl ist 
nicht vorgesehen. Die genaue Geometrie unterliegt der Geheimhaltung (Rolls-Royce Deutsch- 
land, 2017). 

Neben der Geometrie wird der Betriebspunkt der Kühlluft variiert. Ausgehend vom Nenn- 
betriebspunkt ¢ = 3,4% wird der Kühlluftmassenstromanteil auf &j0 = 120 %¢ sowie Con = 
90 %G eingestellt. 


3.3 Erreichte Ähnlichkeit und Betriebsparameter 


Die Inbetriebnahme und Bestimmung der Betriebsparameter gliedert sich in zwei Teilschritte. 
Zuerst wird die externe Druckverteilung am Prüfling ohne Zuführung von Kühlluft untersucht. 
Anschließend werden am vorher bestimmten Betriebspunkt Ölanstrichfarben eingesetzt, um das 
Strömungsfeld qualitativ zu beurteilen. Hier wird sowohl der Fall ohne Kühlung als auch der 
Fall mit Kühlung analysiert. 

Für die aerodynamischen Untersuchungen wird eine spezielle Druckmessschaufel verwendet. In 
diese Schaufel sind, anstelle der Kühlgeometrie, statische Druckmessbohrungen auf der Ober- 
fläche implementiert, welche die Bestimmung der isentropen Machzahl entlang des Profils er- 
möglichen. Diese Machzahlverteilung um das Profil kann bei 10 %, 50 % und 90 % Schaufel- 
höhe ermittelt werden. Auf jeder Schaufelhöhe sind zwölf Messstellen vorhanden, die jeweils 
Erste liegt auf der geometrischen Staulinie (Erfassung des Totaldrucks, May Lg = 0). Entlang 
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(a) Isentrope Machzahl bei 10 % Schaufelhöhe (b) Isentrope Machzahl bei 90 % Schaufelhöhe 


Abbildung 3.5: Verteilung der isentropen Machzahl, weitere Positionen 


der Saugseite befinden sich jeweils sieben Bohrungen °, auf der Druckseite jeweils weitere vier. 
Die so bestimmten Machzahlverteilungen können durch Verstellung der Tailboards gezielt be- 
einflusst werden. Hierbei werden zwei Ziele verfolgt: Das Erreichen der korrekten Abström- 
machzahl May > und das Erreichen der korrekten Verteilung der isentropen Machzahl entlang 
des Profils auf allen drei Schaufelhöhen. Die resultierende geometrische Verstellung (Winkel 
und Verschiebung) wird nicht eingeschränkt oder überprüft. Als Referenz für die Profilvertei- 
lung dient eine numerische Berechnung der Verteilung für die unskalierte, reale Maschinen- 
geometrie mit Rotation, Kühlluftausblasung und Schaufelspitzenspalt (Daten bereitgestellt von 
Rolls-Royce Deutschland). Die Verteilung auf dem Mittelschnitt ist in Abbildung 3.4 gezeigt, 
die erreichte Übereinstimmung ist sehr gut. In den Messdaten ist über die gesamte Drucksei- 
te und im vorderen Bereich der Saugseite eine leichte Erhöhung der Machzahl zu erkennen. 
Dies muss in Kauf genommen werden, um die korrekte Abströmmachzahl May. zu gewähr- 
leisten. Eine weitere Abweichung ergibt sich im engsten Querschnitt und auch leicht stromauf 
des Selbigen. Dies ist der Bereich höchster Massenstromdichte, minimale geometrische Abwei- 
chungen haben großen Einfluss auf die resultierende Machzahl. Hier befindet sich der Übergang 
von der festen Seitenwand auf das bewegliche Tailboard (siehe Abbildung 3.1), was zu einer 
solchen geometrischen Abweichung zwischen Maschine und Prüfstand führt. 


Die Verteilungen bei 10% Schaufelhöhe und 90 % Schaufelhöhe, gezeigt in Abbildung 3.5, 
konnten ebenso sehr gut an die Referenz angepasst werden. Die Beeinflussung dieser beiden 
Verteilungen ist durch das starre Tailboard gekoppelt, es musste eine mittlere Einstellung ge- 
funden werden. Auf 10 % Schaufelhöhe, nahe des Schaufelfußes, verbleibt eine Abweichung 
der Messdaten zu niedrigeren Machzahlen im Bereich der Vorderkante, speziell entlang der 
Druckseite. Hier befindet sich in der Maschine ein Spalt, aus dem Sperrluft ausgeblasen wird. 
Dieser wird nur in der CFD berücksichtigt und führt lokal zu höheren Massenströmen im Be- 
reich der nabenseitigen Endwand. 

Auf 90 % Schaufelhöhe, nahe der Schaufelspitze, ist die größte Abweichung stromauf des engs- 
ten Querschnittes an der Saugseite zu erkennen. Aus den bereitgestellten Informationen zu den 


5 Die letzte Bohrung auf der Saugseite liegt unmittelbar stromauf der Hinterkante und dient zur Bestimmung der 
Abströmmachzahl. 


3.3 Erreichte Ahnlichkeit und Betriebsparameter 45 


(a) Verteilung ohne Filmkühlung auf glatter Ober- (b) Verteilung mit Filmkühlung 
fläche 


Abbildung 3.6: Ölanstrichvisualisierung der Druckseite 


CFD-Vergleichsdaten und der bekannten Geometrie des Prüfstandes lässt sich hierfür keine Er- 
klärung ableiten. Da die Schaufelspitze der Prüflinge aus fertigungstechnischen Gründen jedoch 
leicht verändert und der Schaufelspitzenspalt entfernt wurde, wird eine geringe Änderung der 
Sekundärströmungsfelder im Bereich der Schaufelspitze vermutet. 

In einem finalen Schritt werden die Drücke der Druckkontrollbohrungen (jeweils sechs pro 
Endwand, siehe auch Abbildung 3.1) am vorher bestimmten Betriebspunkt aufgezeichnet. Die- 
se Druckverteilungen auf den Endwänden werden nach Entfernen der Druckmessschaufel als 
Referenz verwendet. So kann auch bei den thermischen Messungen, bei denen der Einsatz von 
Druckmessbohrungen auf der Schaufeloberfläche nicht möglich ist, der Betriebspunkt kontinu- 
ierlich überwacht werden. Weiterhin ist eine Anpassung des Heißgasmassenstroms möglich, die 
zur Kompensation der durch die Filmkühlung zugemischten Luft nötig ist. 


Neben der quantitativen Bestimmung des Betriebspunktes, der auf der Erfassung von isentropen 
Machzahlen an diskreten Messpunkten basiert, wurden Ölanstrichfarben eingesetzt. Durch die 
Verteilung des hochviskosen Öls auf der Schaufel können zwar nur qualitative Aussagen über 
das Strömungsfeld getroffen werden, dafür wird eine vollflächige Visualisierung möglich. Die- 
se Informationen können bei der Beurteilung der thermischen Daten verwendet werden, um die 
einzelnen physikalischen Effekte und deren Einfluss auf die resultierenden Oberflächentempe- 
raturverteilungen qualitativ zu trennen. Hierfür wurden beim Nennbetriebspunkt zwei Versuche 
durchgeführt: 
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1. Prüfling mit Filmkühlungsgeometrie, Kühlung ebenfalls bei Nennbetriebspunkt. 
Vergleichbar mit der späteren Messsituation. 


2. Keine Kühlung, Prüfling mit glatter Außenfläche ohne Filmkühlungsgeometrie. 


Die Prüflinge werden hierbei an der Vorderkante mit dem Öl benetzt und dann der Strömung 
ausgesetzt. Sobald keine weitere Veränderung der Ölverteilung zu beobachten ist, werden die 
Prüflinge entnommen und umgehend fotografiert. Diese Bilder können mit dem in Kapitel 4 prä- 
sentierten Verfahren auf die Prüflingsgeometrie projiziert und so im dreidimensionalen Raum 
bewertetet werden. Die wichtigsten Erkenntnisse dieser Untersuchungen werden im Folgen- 
den kurz zusammengefasst. Es werden zunächst nur Bilder‘ betrachtet, die genaue Projektion 
auf die Prüflingsgeometrie und der Vergleich mit thermischen Strukturen wird in Abschnitt 6.1 
durchgeführt. 


In Abbildung 3.6 sind die Daten der Druckseite für beide Versuche dargestellt. Der Versuch mit 
glatter Schaufel zeigt eine anliegende Strömung. Eine leichte Akkumulation von Öl ist stromab 
der Vorderkante bis zu einer Höhe von ca. einem Drittel der Schaufelhöhe zu erkennen, die 
auf eine leichte Instabilität in diesem Bereich deutet. Der druckseitige Ast des Hufeisenwirbels 
kann nicht visualisiert werden. Es wird vermutet, dass dieser durch den starken Druckgradienten 
sofort in Richtung der gegenüber liegenden Saugseite transportiert wird. Die Staulinie, erkenn- 
bar als helle Linie und gestrichelt eingezeichnet in Abbildung 3.6(b), ist in guter Deckung mit 
der erwarteten Position. Weiterhin kann eine leichte Ablenkung der Strömung zu niedrigeren 
Radien beobachtet werden. Dies ist auf die fehlenden Zentrifugalkräfte zurückzuführen, die bei 
der ursprünglichen Auslegung des Profils berücksichtigt wurden. Mit Filmkühlung zeigt sich 
ein vergleichbares Ergebnis. Die Umlenkung zu niedrigeren Radien wird durch die Filmkühl- 
strahlen noch deutlicher sichtbar. 

In Abbildung 3.7 wird für den gekühlten Versuch 1 die Verteilung an der Vorderkante, Ab- 
bildung 3.7(a), sowie stromab entlang der Saugseite bis zu den Filmkühlbohrungen, Abbil- 
dung 3.7(b), gezeigt. Wieder lässt sich eine anliegende Strömung erkennen, die stromauf der 
Filmkühlbohrungen auf der Saugseite durch deren Versperrung aufgestaut und nahe der jewei- 
ligen Endwand beschleunigt wird. An der gehäuseseitigen Endwand zeigt sich eine schwache 
Staulinie, die auf einen instabilen Hufeisenwirbel hindeutet. Die Staulinie des nabenseitigen 
Hufeisenwirbels ist deutlich stärker ausgeprägt. Dies ist in Einklang mit den in Abschnitt 2.1 
gezeigten Strömungsstrukturen. Anhand der Bilder der verbleibenden Saugseite stromab der 
Filmkühlbohrungen, dargestellt in Abbildung 3.8, können mehrere Effekte gezeigt werden. Im 
Bereich des engsten Querschnittes kann an der glatten Schaufel, Abbildung 3.8(a), ein Gebiet 
der Strömungsablösung beobachtet werden, das sich aufgrund der Lage sowie Größe mit der 
in Abschnitt 2.1 beschriebenen laminar separation bubble erklären lässt. Sobald der Prüfling 
mit Filmkühlung, Abbildung 3.8(b), untersucht wird, verschwindet das Ablösegebiet auf dem 
Großteil der Schaufel. Die erhöhte Turbulenz der Grenzschicht durch die eingebrachte Störung 
durch Filmkühlung stabilisiert diese. Es verbleiben eine schwache Instabilität (Ölakkumulation 
an der Position der Ablöselinie, rote Linie) und ein lokal begrenztes Ablösegebiet (rot gepunk- 
teter Bereich). An der Hinterkante wurde zusätzlich die nach Sharma et al. (1987) berechnete 


6 Die Ausrichtung der gezeigten Ergebnisse ist zur besseren Orientierung zusätzlich in Abbildung 6.3, Seite 80 
eingezeichnet. 
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(a) Verteilung mit Filmkühlung, Vorderkante (b) Verteilung mit Filmkühlung, vordere Saugseite 


Abbildung 3.7: Ölanstrichfarben im vorderen Schaufelbereich 


Eindringtiefe Zgharma eingetragen. Diese ist in guter Übereinstimmung mit der durch die Öl- 
farben visualisierten Wirbelposition, so dass sich die Grenzlinien im hinteren Schaufelbereich 
dem Hufeisenwirbel zuordnen lassen. Radial innen weiter stromauf ist eine zusätzliche definier- 
te Staulinie (blau) zu erkennen. Diese ist in guter Übereinstimmung mit dem in Abbildung 2.1 
gezeigtem Gebiet der Entstehung des Passagenwirbels durch die Interaktion der beiden Äste 
des Hufeisenwirbels (druckseitiger Ast von der nächsten Schaufel). Radial außen ist die Wirbel- 
staulinie im stromauf liegenden Schaufelbereich nicht klar zu erkennen. Mutmaßlich interagiert 
der Passagenwirbel schwach mit dem saugseitigen Ast des Hufeisenwirbels. Diese stabilere 
Ausprägung des nabenseitigen Hufeisenwirbels im Vergleich zum gehäuseseitigen sowie die 
Bildung eines ausgebildeten Passagenwirbels an der Nabenseite konnten auch numerisch von 
Stichling (2018) gezeigt werden. 


Es konnte damit nachgewiesen werden, dass die Schaufelumströmung wie von Literaturdaten 
prognostiziert ausfällt. Die in Abschnitt 2.1 erläuterten und in Abbildung 2.1 qualitativ gezeig- 
ten Strukturen konnten identifiziert werden. Die positive radiale Beschleunigung durch die ge- 
krümmte, nabenseitige Endwand führt dazu, dass die nabenseitigen Wirbel, wie im rotierenden 
Fall, stabilisiert und in den Kanal getragen werden. Diese Beschleunigung ist deutlich niedriger 
als im eigentlichen rotierenden System, wie sich an den zu niedrigeren Radien umgelenkten 
Filmkühlstrahlen auf der Druckseite zeigt. 
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(a) Verteilung ohne Filmkühlung auf glatter Ober- (b) Verteilung mit Filmkühlung 
fläche 


Abbildung 3.8: Ölanstrichvisualisierung der Saugseite 


Unterschiede zwischen Maschine und Prüfstand 


Die für diese Untersuchungen wichtigen aerodynamischen und thermischen dimensionslosen 
Kennzahlen konnten, wie in Abschnitt 3.1 und Abschnitt 3.2 beschrieben, weitestgehend ein- 
gehalten werden. Es musste auf eine Rotation des Prüfstandes verzichtet werden. Der geringe 
Einfluss der fehlenden Scheinkraftdichten, gerade bei komplexen Innenströmungen, wurde im 
vorherigen Abschnitt beschrieben. Die Abweichungen werden daher im Folgenden vernachläs- 
sigt. 

Durch die Skalierung wurde die Heißgastemperatur abgesenkt. Außerdem müssen die Prüflinge 
für den Einsatz der Infrarotthermographie mit einer Oberflächenbeschichtung versehen werden. 
Diese Effekte wurden in den dimensionslosen Kennzahlen des Wärmeübergangs berücksichtigt. 
Ein weiterer Effekt ist die in der Maschine an die Schaufel übertragene Strahlungsleistung. So 
muss aufgrund der stark nichtlinearen Proportionalität dieses Wärmestroms (Grad X T*) und der 
veränderten radiativen Eigenschaft der Oberfläche (es ~ ag 7) mit einer Steigerung desselben 
gerechnet werden. Dies führt zu einer höheren thermischen Belastung des Schaufelmaterials im 
späteren Einsatz. 


Es wurde ein neuer Prüfstand konzipiert und erfolgreich in Betrieb genommen, der erstmals 
die Einhaltung aller für den Wärmeübergang wichtigen dimensionslosen Kennzahlen garan- 
tiert. Diese Ähnlichkeit konnte in aerodynamischen Messungen nachgewiesen werden. Weiter- 
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hin wurden die verbleibenden Unterschiede zu einer Turbinenschaufel, die in einer Maschine 
eingesetzt wird, erläutert sowie der Resteinfluss dieser Unterschiede qualitativ abgeschätzt. Mit 
dem präsentierten Prüfstand können demnach thermische Messungen durchgeführt werden, die 
die Beurteilung der Oberflächentemperatur verschiedener Prüflinge zum Ziel haben. Aufgrund 
der Ähnlichkeit können diese Daten für finale Auslegungen herangezogen werden. Zur Durch- 
führung der Messungen mit vertretbarem Messfehler müssen, aufgrund der komplexen Geo- 
metrien und ungünstigen Strahlungsverhältnisse, neue Kalibrierverfahren entwickelt werden. 
Diese neuen Verfahren werden in den folgenden zwei Kapiteln präsentiert. 


4 Modellierung und Positionsbestimmung der Kamera 


Wichtige Grundlage für die später in Kapitel 5 vorgestellten Korrekturverfahren ist die genaue 
Modellierung und Positionsbestimmung der Kamera. Die nötigen Annahmen sowie Berech- 
nungsverfahren werden in folgendem Kapitel präsentiert. Zuerst wird auf das Modell der Loch- 
kamera eingegangen, das eine hinreichend genaue Modellierung der verwendeten Infrarotka- 
mera erlaubt. In einem weiteren Schritt wird darauf aufbauend ein neues Verfahren präsentiert, 
das die Bestimmung des Kameramodells anhand eines Bildes und der Geometrie des Prüflings 
erlaubt. Die Anwendbarkeit auf die gegebene Geometrie wird dargelegt. 


4.1 Modell der Lochkamera 


Das Modell der Lochkamera beschreibt die einfachste Umsetzung einer Kamera. Strahlen aus 
dem dreidimensionalen Objektraum fallen durch eine infinitesimal kleine Öffnung auf einen 
hinter dieser Öffnung liegenden Schirm, es entsteht ein zweidimensionales Abbild im Bildraum. 
Dies entspricht einer Zentralprojektion mit der Öffnung als Projektionszentrum und der Skalie- 
rung f, wobei f den Abstand zwischen Öffnung und Schirm in Normalenrichtung definiert. 
Durch geometrische Analyse lässt sich 


(X,Y,Z)' + (fX/Z,fY/Z,f) > (x,y,—) (4.1) 


ableiten. Hierbei (und im Folgenden) beschreibt der Vektor X= (X,Y, Ze die Koordinaten im 
Objektraum und der Vektor ¥ = (x,y)? die Koordinaten im Bildraum. Aufgrund der Skalierung 
ist die Beschreibung der Transformation durch eine einzelne lineare Abbildung im euklidi- 
schen Raum nicht möglich. Für die mathematische Beschreibung dieser Projektion ist daher 
die Einführung des projektiven Raumes und der homogenen Koordinaten zweckmäßig. Für 
die Berechnung werden die Koordinatenvektoren um den Skalierungsfaktor W bzw. w ergänzt. 
Dies bedeutet, dass Punkte, die im euklidischen Raum im Unendlichen liegen, nicht geson- 
dert behandelt werden müssen. Das Konzept der Parallelität (Schnittpunkt im Unendlichen) 
wird fallengelassen. Euklidisch parallele Linien schneiden sich im Unendlichen, erfasst durch 
W = w =O. Außerdem lassen sich komplexe euklidische Transformationen, die Skalierungen 
und Translationen beinhalten, mittels einer einzelnen projektiven Transformationsmatrix P ana- 
log zu linearen Abbildungen im euklidischen Raum darstellen. 

Der Wechsel von euklidischen in homogene Koordinaten ergibt sich durch das Ergänzen des 
Koordinatenvektors mit dem Skalierungsfaktor W = w = 1, die euklidischen Koordinaten erge- 
ben sich aus den Homogenen zu (X/W,Y/W,Z/W) bzw. (x/w,y/w). Punkte im Unendlichen 
(W = w = 0) können so weiterhin identifiziert werden. Die Projektion der Lochkamera nach 
Gleichung 4.1 kann dann als lineare Abbildung dargestellt werden: 


x f 0 : 5 
y| = f (0) =P (4.2) 
w 1 0 2 2 

W W 
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Für die Erweiterung auf FPA-Kameras mit Sensoren mit diskreten Elementen (Pixeln) werden 
wichtige Begrifflichkeiten eingeführt: 


e Das Kamerazentrum Č beschreibt das Projektionszentrum der Zentralprojektion. Wird 
im Folgenden von der Position der Kamera gesprochen, ist die Lage des Kamerazentrums 
gemeint. 


e Der Vektor durch das Kamerazentrum und senkrecht auf die Bildebene wird als optische 
Achse bezeichnet. Wird im Folgenden von der Ausrichtung der Kamera gesprochen, ist 
die Richtung der optischen Achse gemeint. 


Der Schnittpunkt der optischen Achse mit der Bildebene wird als Bildmitte ¢ bezeichnet. 
Diese Bildmitte liegt nicht zwingend in der Mitte des Sensors, der später in der Bildebene 
liegt. 


Der Bildabstand b und der Objektabstand g bezeichnen die Abstände von der Bildebene 
zum Kamerazentrum bzw. vom Objektpunkt zum Kamerazentrum, jeweils projiziert auf 
die optische Achse. 


e Die Brennweite f beschreibt den Abstand des Projektionszentrums zur Bildebene b. 
Wenn ein Objektiv verwendet wird, wird hierunter häufig die physikalische, durch das 
Linsensystem festgesetzte Brennweite verwendet. Diese ist aber nur bei Fokussierung 
auf © gleich dem Bildabstand b. 


So ergibt sich für beliebige FPA-Kameras die Kameramatrix P zu: 


xX fs cx O Y 
y| = fc OO], (4.3) 
w 1 0 

W 


Hierbei beschreibt s die Scherung eines jeden Pixels und damit des gesamten Detektors (im 
Allgemeinen gilt s = 0 aufgrund rechtwinkliger Detektoren). Die intrinsische Kameramatrix k 
kann hier direkt aus der linken 3 x 3-Teilmatrix abgelesen werden. Soll nun noch berücksichtigt 
werden, dass die Kamera in einem beliebig gewählten Koordinatensystem nicht im Ursprung 
liegt, muss die Transformation um eine Rotation R und eine Translation 7 bzw. die extrinsische 
Kameramatrix r = |[R|r]ergänzt werden: 


P=K{Rif]. (4.4) 


Diese Zusammenhänge sind graphisch in Abbildung 4.1 gezeigt. Hierbei wird eine exemplari- 
sche Strahlungsverteilung von einer komplexen Oberfläche mittels der gezeigten Transforma- 
tionen auf ein diskretes Gitter (den Bildsensor) abgebildet. Die Bildebene wird im Bildabstand 
b vor das Projektionszentrum C gelegt. Dies ist geometrisch äquivalent zu dem bei Kameras 
üblichen Aufbau mit dem Bildsensor hinter dem Projektionszentrum. Die Rotation des Bildes 
um 180° entfällt". 


! Die nötige Rückrotation bei Lage der Bildebene hinter dem optischen Zentrum C wird von den meisten Kameras 
intern vor Bereitstellung der Rohdaten durchgeführt. 
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Abbildung 4.1: Zentralprojektion mittels Lochkameramodell. Koordinatensystem O’ zur bes- 
seren Anschauung verschoben, es gilt O’ = C 


Es ist zu beachten, dass die Kameramatrix (als Kombination aus intrinsischer und extrinsischer 
Teilmatrix) des Lochkameramodells bei Abbildung vom dreidimensionalen Objektraum in den 
zweidimensionalen Bildraum nicht invertierbar ist. Diese mathematische Tatsache lässt sich 
leicht geometrisch erklären: Zwar kann einem Punkt im Raum exakt ein Punkt im Bild zuge- 
ordnet werden, ein Punkt im Bild definiert jedoch einen kontinuierlichen Strahl, der durch das 
Kamerazentrum sowie den Punkt selbst verläuft. Die Tiefeninformation im Objektraum geht 
demnach verloren. Bei Verwendung von nur einer Kamera kann die Rekonstruktion der Szene 
aus Bilddaten nur erfolgen, wenn Annahmen über die Objekte in der Szene getroffen werden. 
Für die Korrekturverfahren, die im Rahmen dieser Arbeit vorgestellt werden, ist die Annahme 
mit geringstem Fehler eine Diskretisierung der zu untersuchenden Oberfläche. Diese Diskre- 
tisierung wird so fein gewählt, dass die lokale Auflösung unter der der Kamera liegt. Anhand 
dieser Diskretisierung lassen sich Bilddaten in den Objektraum transformieren. Das genaue 
Verfahren wird in Abschnitt A.1 beschrieben. 


4.2 Einschränkungen des Lochkameramodells 


Das vorgestellte Modell der Lochkamera ist für die Fragestellungen dieser Arbeit ausreichend 
umfangreich. Für die weitere Anwendung müssen Einschränkungen beachtet werden, um eine 
genaue Modellierung der verwendeten Kamera zu erreichen. 

So wird von einer idealen Abbildung ausgegangen. Eventuelle radiale oder tangentiale Verzer- 
rungen, die durch Objektivfehler entstehen, können nicht berücksichtigt werden. Eine Möglich- 
keit zur Verwendung der vorgestellten Modelle bei vorhandener Verzerrung ist die Korrektur im 
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Bildraum. Dies ist bei dem verwendeten Kamerasystem nicht notwendig, da diese Verzerrungen 
sehr gering sind (Willmann, 2014). 

Ein weiterer Unterschied ergibt sich in der Definition der Brennweite. Diese ist beim Modell 
der Lochkamera immer der Abstand b zwischen dem optischem Zentrum C und der Bildebene, 
eingetragen in Abbildung 4.1. Bei Verwendung einer realen Kamera mit Sammellinse, bezie- 
hungsweise eines komplexen Objektivs, beschreibt die Brennweite tiblicherweise die physika- 
lische Eigenschaft der Sammellinse oder des Objektivs und damit den Abstand des optischen 
Zentrums C zur Bildebene bei Fokussierung auf Unendlich. Dies ist in den meisten Fallen eine 
gute Näherung: Der Auszug yp, der bei realen optischen Systemen mit fixer Bildebene benötigt 
wird, um Punkte näher als Unendlich in der Bildebene abzubilden, ist wesentlich kleiner als die 
Brennweite. Zusätzlich ist der Objektabstand meist wesentlicher größer als der Bildabstand 


Wf ER. (4.5) 


Der Winkelfehler durch Verlagerung des optischen Zentrums in Richtung der optischen Achse 
ist daher auch gering. Bei ansteigender Brennweite steigt zwar der nötige Auszug, aber auch 
Bildabstand b und Objektabstand g. 

Ist ein großer Abbildungsmaßstab B gefordert, werden diese Bedingungen verletzt. Durch Ver- 
wendung von Zwischenringen? steigt der Auszug um das Maß der verwendeten Verlängerung 
Xzr an, bei Auszug um die Brennweite xzr = f wird genau der Abbildungsmaßstab B = 1 
erreicht. Daraus folgt also: 


e Wird das Lochkameramodell benutzt, um eine reale Kamera im Nahbereich zu berechnen, 
muss die Brennweite der intrinsischen Matrix auf das Maß der physikalischen Brennweite 
plus Auszugsverlängerung gesetzt werden: 


fe =f +x2R. (4.6) 


e Der Fehler durch den nicht genau bekannten Auszug yp ist weiterhin gering. Glei- 
chung 4.5 gilt auch nach Ersetzten von f durch fx. Eine genaue Bestimmung von 
xr müsste für jede Aufnahme durchgeführt werden und ist praktisch nur mit großem 
Aufwand durchzuführen. Da die Kamera üblicherweise weder an der Nah- noch an der 
Ferneinstellgrenze fokussiert ist, wird der Fehler durch Setzen von 


XF,max 
2 


fk =f+4zR+ (4.7) 


weiter reduziert. 
e Eine Berechnung der objektseitigen Schärfeebene ist durch die Linsengleichung möglich, 


wobei die Bildweite b durch die geschätzte Brennweite der Lochkamera fx ersetzt werden 
muss: 


1 1 | 
—+-=-, (4.8) 
fk & f 

2 Im Rahmen dieser Untersuchungen werden zwei verschiedene Zwischenringe mit den Längen xzr ı = 0,5” 


sowie xzr = |” eingesetzt. 
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Nach der Schätzung der Kamera durch die später vorgestellten Methoden kann die Brenn- 
weite der Lochkamera von der Vorgabe abweichen und der Auszug xr angepasst worden 
sein (siehe Gleichung 4.11c). Da aber gerade im Nahbereich kleine Auszugsänderungen 
Axr einen großen Einfluss auf die Lage der objektseitigen Schärfeebene haben, bleibt die 
Berechnung ungenau. 


4.3 Bestimmung der Kameraparameter 


Zur vollständigen Beschreibung der Kamera ist die Kameramatrix P ausreichend. Diese Matrix 
hat im Allgemeinen elf Freiheitsgrade: drei für die Rotation, drei für die Translation (extrin- 
sisch), zwei Brennweiten fx und fy, die Bildmitte č und letztlich die Schiefe s. Der zwölfte 
Eintrag verbleibt als abhängiger Skalierungsfaktor. Die Bestimmung erfolgt nach dem DLT?- 
Algorithmus (Hartley et al., 2003). So soll die Kamera Gleichung 4.2 für alle Punkte erfüllen. 
Im projektiven Raum kann dies durch Kolinearität erfüllt werden, was zu einem linearen Pro- 
blem führt: 


xxPX=0 bzw. Ap=0 mit p=| : (4.9) 


Durch Bestimmung von korrespondierenden Punkten in Bild und Objektraum {xi Xi} lässt sich 
die Matrix A bilden. Hier werden aus jedem Punkt zwei Zeilen abgeleitet, also zwei Freiheits- 
grade bestimmt. Zur Lösung sind mindestens 5,5 Bildpunkte notwendig. Die Lösung des Glei- 
chungssystems erfolgt numerisch vorteilhaft durch Singulärwertzerlegung von A. Nach erfolg- 
ter Bestimmung von P lassen sich die gesuchten Kameraparameter ableiten. Zusätzlich werden 
die Reihen von P als p;, die Zeilen von P als p! und die Matrix M3x3 mit P = [M|p*] eingeführt. 


e Das Kamerazentrum Č ergibt sich aus PC = Ü durch Singulärwertzerlegung von P. 


e Die Richtung der optischen Achse ergibt sich zu OA = det(M)mz mit mz als dritte Zeile 
von M. Hierbei bestimmt 173 wegen P3X = 0 die Richtung und det(M) die Orientierung 
der Achse. 


e Die intrinsische Kameramatrix k folgt aus Gleichung 4.4 durch RO-Zerlegung von M. 
Durch Multiplikation mit dem Pixelabstand ¢, = &, = ¢ lässt sich die Brennweite in mm 
bestimmen. 


Diese Bestimmung der Kameramatrix setzt voraus, dass die gewählten Punktkorrespondenzen 
entweder fehlerfrei sind oder dass eine starke Überbestimmung vorliegt, also wesentlich mehr 
Korrespondenzen als nötig vorhanden sind. Beides kann nicht gewährleistet werden. Eine feh- 
lerfreie Auswahl ist bei gegebener Auflösung nicht möglich, eine Überbestimmung, wie sie in 
anderen Anwendungsbereichen praktisch umgesetzt wird, setzt oft mehrere tausend Korrespon- 


3Direkte Lineare Transformation 
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denzen voraus. Auch ein nicht systematisches Probieren ist bei gegebener Problemgröße nicht 
zielführend.* Es müssen daher neue Möglichkeiten gefunden werden, die Kameramatrix mit 
wenig Korrespondenzen robust zu bestimmen. 


4.3.1 Genetischer Algorithmus 


Wenn der komplexe funktionale Zusammenhang eines Problems nicht bekannt ist, eignen sich 
genetische Algorithmen zum Lösen. Diese optimieren durch geschicktes Probieren ohne genaue 
Kenntnis der Übertragungsfunktion und können so schnell zu guten Lösungen finden. Hierbei 
wird eine mögliche Lösung als Individuum bezeichnet. Jedes dieser Individuen besitzt ein ein- 
zigartiges Genom, welches die nötigen Eingangsgrößen des Problems enthält. Die Lösung des 
Problems wird dann für jedes Individuum bestimmt und die Qualität der Lösung diesem als 
Eignung zugeordnet. Dies wird für eine Menge an Individuen (eine Population) durchgeführt. 
In einem weiteren Schritt werden nur die bestgeeigneten Individuen der Population genutzt, um 
neue Permutationen dieser Individuen, eine neue Population, durch zufällige Veränderungen 
abzuleiten. Mit dieser neuen Population (einer neuen Generation) beginnt der Prozess erneut 
(Mitchell, 1998). 

Nachteile dieser Algorithmen sind eine schwierige Beurteilung der Konvergenz (sowohl der 
Verlauf als auch das finale Ergebnis) und eine eventuell erhöhte Rechenzeit durch unsystemati- 
sche Optimierung. 

Diese Art der Bestimmung der Kameramatrix wird erstmals von Kneer (2017) eingesetzt. Er 
definiert die Objektkoordinaten X; der Korrespondenzen als unveränderlich und führt die Bild- 
koordinaten x; als Genom eines jeden Individuums ein. Auf Basis einer ersten manuellen Aus- 
wahl der Bildkoordinaten wird durch zufällige Veränderung dieser eine Menge an j Individuen 
erzeugt. So lässt sich für jedes Individuum eine Kameramatrix P; bestimmen. Als einzigen Op- 
timierungsparameter definiert Kneer nach Hartley et al. den Reprojektionsfehler 


(4.10) 


I 


Der von Kneer (2017) eingesetzte genetische Algorithmus wird von einem proprietären 
MATLAB-Paket ausgeführt. Weitere Informationen über die Methoden zur Reproduktion 
und Selektion des Algorithmus können der Dokumentation (Stand R2018A) entnommen 
werden: Reproduktion kann durch elitäre Weitergabe, Mutation oder Kreuzung an einzelnen 
Teilen des Genoms (sog. one-point crossover) erfolgen. Selektion erfolgt gewichtet durch 
den Qualitätswert Arp. Zwar lassen sich die einzelnen Teilschritte konfigurieren, der Aufbau 
eines erweiterten genetischen Algorithmus mit gegebenenfalls mehreren Qualitätskriterien pro 
Individuum ist aber nicht möglich (Matlab, 2018). 


4 Mit beschriebenen Freiheitsgraden und dem verwendeten, verhältnismäßig niedrig auflösenden Kamerasys- 

tem ergibt sich die Anzahl möglicher Lösungen, ähnlich einem beliebten Glücksspiels, zu ($): Hier gilt 
für n im schlechtesten Fall n = no = npixeı © 350.000 und selbst bei sehr genauer Wahl der Punkte immer 
noch n = nı = nro & 10.000 unter Annahme von einem 4mm x 4mm großen Toleranzfenster. Die zuge- 
hörigen Lösungsmengen müssen komplett ausgewertet werden und ergeben sich zu N(no) = @(10°°) sowie 


N(n;) = @(107!). 
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Bei komplexen Geometrien und wenigen Korrespondenzen findet dieser Algorithmus keine zu- 
friedenstellende Lösung und muss daher erweitert werden. Als Grundlage dient der (u +A)- 
Algorithmus, bei dem in jedem Iterationsschritt aus u bestgeeigneten Elternindividuen A Neue 
erzeugt werden. Aus diesen u + A Individuen werden für die nächste Generation wieder nur die 
u bestgeeigneten weiterverwendet. Beyer et al. (2002) geben einen Überblick über die verschie- 
denen Methoden und fassen die Vorteile der (u + A )-Strategie zusammen. Bei dieser Version 
ist Divergenz nicht möglich, maximal eine Stagnation bei lokalen Minima. Die Evaluation von 
vielen Individuen pro Generation ist gut parallelisierbar und damit numerisch schnell. Der so 
umgesetzte hohe Selektionsdruck 4/1 kann die benötigten Generationen bis zum Erreichen der 
gewünschten Lösung verringern. Der Ablauf ist im Folgenden beschrieben, ein Ablaufplan ist 
im Anhang in Abbildung A.3 zu finden. 


0. Als Eingangsgröße dient weiterhin eine erste Auswahl an Bildpunkten x;, wobei die Ob- 
jektkoordinaten X; festgehalten werden. Die Erzeugung der Initialpopulation, bestehend 
aus u Individuen, erfolgt durch elementweise Variation des Genoms x;. Mit 50 % Wahr- 
scheinlichkeit wird ein Element variiert. Falls variiert wird, erfolgt eine Verschiebung um 
Apx = U (—t, +t), wobei aus der Gleichverteilung mit Toleranz t gezogen wird. Die Tole- 
ranz (Einheit px) sollte abhängig von der Auflösung des Kamerasystems im Objektraum 
gewählt werden. Geeignet sind Werte, die ca. 1 mm entsprechen. Bei größerer Unsicher- 
heit in der Wahl der Punkte kann dieser Wert auch erhöht werden. 


1. Die Population der j Individuen wird evaluiert: Es werden zuerst alle Kameramatrizen 
P; nach Gleichung 4.9 berechnet. Eine reine Optimierung des Reprojektionsfehlers nach 
Kneer kann nicht sicherstellen, dass die restlichen Kameraparameter, also die Position 
und die intrinsische Matrix, korrekt bestimmt werden. Daher werden weitere Zielgrößen 
eingeführt, die die Optimierung dieser Größen sicherstellen: 


Ac =|le-g;|| = ]|E- (Ms); Gta) 
Ag j = arctan((12/k22) j) (4.11b) 
Aaj = fx —||kj-(64)]]- (4.1 1c) 


Die erste Größe nach Gleichung 4.1 1a beschreibt die Abweichung der Bildmitte von der 
erwarteten Bildmitte. Bei einer perfekten Kamera sollte diese Bildmitte in der Mitte des 
Sensors liegen. Die Zielgröße wird minimiert. Die zweite Größe nach Gleichung 4.11b 
beschreibt den Scherungswinkel des Detektors. Dieser wird minimiert, da die verwende- 
ten FPA-Detektoren rechtwinklig sind. Die letzte Größe nach Gleichung 4.1 1c beschreibt 
die Abweichung des Abstandes von Detektor zu Kamerazentrum vom erwarteten Ab- 
stand nach Abschnitt 4.2. Auch diese Zielgröße wird minimiert. Zusammen mit dem vor- 
her definierten Reprojektionsfehler ergeben sich vier Zielgrößen. Abschließend wird die 
Konvergenz beurteilt. Sollten die definierten Kriterien nicht erfüllt sein, wird mit Schritt 
2 fortgefahren. 


2. Anhand der bestimmten Kameramatrizen P; und der daraus abgeleiteten Zielgrößen 
nach Gleichung 4.10 sowie Gleichung 4.11 müssen aus (A + u) Individuen die u 
besten ausgewählt werden. Hierfür wird der von Deb et al. (2002) vorgeschlagene 
NSGA-II-Algorithmus mit einem verbesserten Sortierverfahren (Fortin et al., 2013) 
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verwendet. Dieser sortiert alle Individuen auf Paretofronten im vierdimensionalen 
Lösungsraum. Falls auf den definierten Fronten nicht genügend Individuen vorhanden 
sind, werden die Übrigen anhand ihres normierten Abstandes zu den umliegenden 
Individuen sortiert und die mit großem Abstand gewählt. Diese Methode hat mehrere 
Vorteile: Es muss kein skalarer Parameter, abgeleitet aus den Zielgrößen, definiert 
werden: Eine echte Pareto-Optimierung wird durchgeführt. Elitäre Individuen werden in 
die nächste Population übernommen. Durch die Verbesserung der Sortierungsverfahren 
ist die Methode, im Vergleich zu anderen Verfahren, schnell. Eine Implementierung des 
NSGA-II-Algorithmus wird von Fortin et al. (2012) zur Verfügung gestellt. 


. Aus den ausgewählten u Individuen muss eine neue Generation von A Indivi- 


duen erzeugt werden. Dies geschieht durch Kreuzung, Mutation oder Weiter- 
gabe. Jedes der A neuen Individuen kann nur durch einen der drei Modi ent- 
stehen, jedes Elternindividuum kann beliebig oft eingesetzt werden. Für die 
Mutationswahrscheinlichkeit P,, = P(Mutation) = 0,4 und die Kreuzungswahr- 
scheinlichkeit Pex = P(Kreuzung) = 0,45 gilt damit Pm + Pe < 1. Die ver- 
wendete Methodik für Kreuzung und Mutation wird in Abbildung 4.2 gezeigt. 
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Abbildung 4.2: Mutation und Kreuzung am Beispiel eines Genoms aus sieben Punkten, 


* 


zweidimensional. Mutierte Größen mit 


Bei einer Kreuzung werden sowohl Ein- als auch Aussprungspunkt zufällig gewählt. 
Nach Übernahme des Genoms von Elternindividuum 1 wird dann der durch die 
Grenzen bestimmte Teil von Elternindividuum 2 ersetzt (sog. two-point crossover). Es 
wird immer der gesamte Punkt (beide Koordinaten) übernommen. Bei einer Mutati- 
on werden, ausgehend vom Elternindividuum 1, eine zufällige Anzahl von Punkten 
(Pini = P(Mutation, Indiv.) = 0,35) verändert. Dies kann auch nur einzelne Koordinaten 
betreffen. Der Betrag der Mutation wird aus einer Gaußschen Verteilung N(v,2) mit 
v=U(-v,,v,), v = 1,5px und £ = 3 gezogen. Bei Weitergabe wird das Elternindivi- 
duum übernommen. Weitergabe kann nur bei Pm + Pex < 1 stattfinden. Dies sichert in 
Verbindung mit dem (u + A)-Schema sowie dem NSGA-II-Algorithmus den Elitismus 
der Evolution. 


4. Der Algorithmus wird bei Schritt 1 fortgesetzt. 
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Nach Beendigung der Evolution verbleiben  paretooptimale Individuen, die jeweils eine Ka- 
meramatrix P definieren. Im vorliegenden Fall kann das Individuum mit dem geringsten Repro- 
jektionsfehler ausgewählt werden. Dies kann bei Betrachtung der anderen Zielgrößen erklärt 
werden: Die Bildmitte weicht aus fertigungstechnischen Gründen und Objektivfehlern von der 
theoretischen Detektormitte ab (Willson et al., 1993). Die Optimierung der Schiefe gelingt sehr 
gut, der Restfehler liegt bei @(10~*). Der Bildabstand des Sensors b bzw. f ist a priori nicht 
genau bekannt (siehe Abschnitt 4.2), eine gewisse Abweichung muss daher toleriert werden. 
Die durch das Genom x; des Individuums mit geringstem Reprojektionsfehler aus den u pa- 
retooptimalen Individuen und den fixen Objektpunkte X; definierte Kameramatrix ist dann die 
optimale Kameramatrix. 


Konvergenzverhalten 


Die Beurteilung der Konvergenz stellt bei genetischen Algorithmen eine Herausforderung dar, 
da die Optimierung nicht systematisch verläuft. So können je nach Güte der ersten Auswahl an 
Punkten schon in den ersten Iterationen optimale Lösungen gefunden werden. Außerdem kann 
der Optimierer, je nach Fall und Mutationsstärke, auch über mehrere Iterationen (im schlech- 
testen Fall auch dauerhaft) bei lokalen Minima stehen bleiben. Dies kann aufgrund des kom- 
plexen funktionalen Zusammenhangs des Problems und der daher unbekannten Lösungsfläche 
nicht verhindert werden. Im vorliegenden Fall wird daher der minimale Reprojektionsfehler 
und dessen arithmetisches Mittel der momentanen Population zur Beurteilung der Konvergenz 
verwendet. Als Abbruchkriterien dienen sowohl ein Schwellwert, der für diese Größen unter- 
schritten werden muss, als auch die Änderungsrate des minimalen Reprojektionsfehlers. Durch 
diese Kombination und die Verwendung der in Unterabschnitt 4.3.1 beschriebenen Methoden 
wird sichergestellt, dass eine verwertbare Lösung gefunden wird. Der Elitismus des Algorith- 
mus sichert weiterhin, dass diese Lösung uneingeschränkt in jeder Generation auftreten kann. 
Als maßgebliche Einflussfaktoren auf den Konvergenzverlauf werden demnach die Mutations- 
stärke v, und die Vielfalt einer Population A untersucht. Hierbei werden als Ausgangswerte 
die als optimal ermittelten Werte v, = 1,5px und u,A = 500,7.000 verwendet. Die Mutation 
wird dann bei konstanter Vielfalt reduziert (Pmi = 0) und verstärkt (v, = 10px). Die Ergebnis- 
se werden in Abbildung 4.3 gezeigt. Der Mittelwert zeigt starke Schwankungen, da dieser von 
wenigen Individuen mit sehr hohem Reprojektionsfehler stark beeinflusst wird. Sobald Mittel- 
wert und Minimalwert für mehrere Generationen einen vergleichbaren Wert annehmen, ist von 
einem erreichten Minimum auszugehen. 

Ohne Mutation, also Pmi = 0, zeigt sich ein frühes Erreichen eines lokalen Minimums mit ver- 
gleichbar hohem Reprojektionsfehler. Mittelwert sowie Minimalwert erreichen denselben Wert, 
der Algorithmus kann das erreichte lokale Minimum ohne Mutation nicht verlassen. Nach ca. 
120 Generationen findet der Abbruch statt. Bei Verwendung einer hohen Mutationsstärke mit 
vy = 10px zeigt sich zunächst ein steilerer Verlauf zu niedrigen Reprojektionsfehlern. Es wer- 
den nach einigen Generationen immer wieder lokale Minima erreicht, da die starke Mutation 
zu großen Sprüngen auf der Lösungsfläche führt und die mutierten Individuen schlechtere Eig- 
nung zeigen. Die Konvergenz ist langsam. Der Abbruch erfolgt nach ca. 230 Generationen. Mit 
einer moderaten Mutationsstärke kann eine kontinuierliche Reduktion des Reprojektionsfehlers 
erreicht werden. Eine kurzzeitige Stagnation an lokalen Minima kann nicht vermieden werden 
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Arp, min: Pmi = 0 
© Arp, min: vy = 10 px 
Arp, 0: vy = 1,5 px 
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Abbildung 4.3: Konvergenzverhalten des genetischen Algorithmus bei Variation der Mutati- 
onsstärke 


(siehe Generation 390 bis 410). Die Lösung bei Abbruch wegen anhaltender Stagnation (ca. 
Generation 540) erreicht einen um eine Größenordnung reduzierten Reprojektionsfehler. 
Weiterhin wird die Vielfalt bei konstanter Mutationsstärke reduziert (A = 700) und erhöht 
(A = 21.000). Hierbei wird der Selektionsdruck #/r konstant gehalten. Die Ergebnisse dieser 
Variation werden in Abbildung 4.4 gezeigt. Bei zu geringer Vielfalt A = 700 ist eine Stagna- 
tion bei hohem Reprojektionsfehler zu beobachten. Bei zu kleiner Population ist die Anzahl 
an Lösungskandidaten zu gering, um effektiv lokale Minima zu überwinden. Die Vielfalt muss 
daher erhöht werden, bis eine hinreichende Qualität der Lösung erreicht wird. Bei Steigerung 
der Vielfalt steigt die Rechenzeit aufgrund des Selektionsverfahrens quadratisch (Deb et al., 
2002). Daher muss ein optimaler Wert gefunden werden, bei dem hinreichend schnell eine aus- 
reichende Konvergenz erreicht wird. Bei einer Steigerung der Vielfalt kann zwar eine schnellere 
Konvergenz in Bezug auf die Anzahl der Iterationen erreicht werden, bezogen auf die benötigte 
Rechenzeit fallen Optimierungen mit zu großen Populationen aber langsamer aus. Dies wird in 
Abbildung 4.4 exemplarisch an der Population für A = 21.000 gezeigt. Die Berechnung wurde 
nach ca. 150 Generationen abgebrochen. Die Rechenzeit lag schon eine Größenordnung höher 
als bei den anderen präsentierten Fällen. 
Für den vorliegenden Fall wurde v, = 1,5px und u, = 500, 7.000 als optimal ermittelt. Mit 
diesen Parametern wird der gewählte Schwellwert für den Reprojektionsfehler von Ap < 0,4px 
über eine Stagnationszeit von 30 Generationen nach vertretbarer Rechenzeit zuverlässig unter- 
schritten. 
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Abbildung 4.4: Konvergenzverhalten des genetischen Algorithmus bei Variation der Vielfalt 


4.3.2 Plausibilitätsprüfung 


Die final berechnete Kameramatrix kann nicht experimentell validiert werden. Das berechnete 

Kamerazentrum ist das Projektionszentrum und liegt an einem unbekannten Ort im optischen 
Pfad des Objektivs. Die tatsächliche Bildmitte ist nicht bekannt. Der variierende Auszug be- 
stimmt Perspektive und Schärfentiefe. Letztere wird im Modell der Lochkamera nicht berück- 
sichtigt. Der Objektabstand g und die Perspektive reagieren insensitiv auf kleine Veränderungen 
des Auszugs. Alle bestimmten Auszüge waren in technisch vertretbarem Rahmen. So kann eine 
berechnete Kamera nur auf ihre Plausibilität überprüft werden. Dies wurde für die berechneten 
Kameras durchgeführt: Die berechnete Position wurde mit einer im Raum gemessenen Position 
unter Annahme eines mittig im Objektiv liegenden Projektionszentrums verglichen. Die Ab- 
weichungen wurden auf IA z | < 4mm bestimmt. Diese Grenzabweichung wird auch später für 
die Fehlerabschätzung eingesetzt. 

Weiterhin kann anhand der aufgenommenen Daten und deren Überlappung mit der Geome- 
trie überprüft werden, ob die Kameramatrix die reale Kamera repräsentiert. Die diskretisierte 
Geometrie wird mittels der Kameramatrix in den Bildraum transformiert. Anhand von mar- 
kanten geometrischen Konturen, die auch in den Strahlungsdaten identifizierbar sind, lässt sich 
die Überlappung beurteilen. Diese Überlappungsbilder sind in Abbildung 4.5 für verschiede- 
ne Kamerapositionen dargestellt. Hier wird im Hintergrund das Strahlungssignal eines Bildes 
der Kamera dargestellt. Dunkle Bereiche, also geringe Oberflächenstrahlungen, lassen sich den 
Filmkühlbohrungen zuordnen. In weiß sind die mit dem mathematischen Modell der Kamera in 
den Bildraum transformierten Kanten der Filmkühlbohrungen überlagert. Es zeigt sich eine her- 
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Abbildung 4.5: Übereinstimmung der projizierten Geometrie mit den Bilddaten 


vorragende Überlappung und Qualität der berechneten Kamerapositionen. Die Restabweichun- 
gen liegen im Bereich von wenigen px. Selbst die in Abbildung 4.5(b) gezeigte Aufnahme, die 
unter unvorteilhaftem Sichtwinkel durchgeführt wurde und aufgrund der geringen Schärfentiefe 
aus mehreren Aufnahmen zusammengesetzt wird (siehe auch Abschnitt 5.1), zeigt, bis auf den 
extremen Randbereich links unten, keine nennenswerten Abweichungen. Weiterhin muss neben 
der Qualität der Kameraprognose auch berücksichtigt werden, dass die Geometrie selbst durch 
das in Abschnitt 3.1 beschriebene Fertigungsverfahren SLM Abweichungen in der Größenord- 
nung des Restfehlers zeigt, die Geometriedaten aber aus dem idealen CAD-Modell stammen. 


Die entwickelte Methode erlaubt auch bei schwierigen geometrischen Bedingungen und weni- 
gen Korrespondenzen eine Schätzung der Kameramatrix. Die erreichte Genauigkeit ist auch 
bei unvorteilhaften Kamerapositionen und Sichtwinkeln sehr hoch. Kompositbilder, die mit 
leicht unterschiedlichen Kamerapositionen aufgenommen wurden, wie es für das später ge- 
zeigte Fokus-Stapeln notwendig ist, können durch eine gemittelte Kameramatrix repräsentiert 
werden. Die so bestimmten Kameramatrizen können demnach für die im Folgenden Kapitel 5 
eingeführten Korrekturverfahren verwendet werden. 


5 Die Auflösung der gezeigten Bilder liegt bei ca. 10 px/mm. Eine zylindrische Bohrung (elliptische Austrittsflä- 
che) hat entlang der kurzen Achse einen Halbdurchmesser von 1,2 mm. 


5 Infrarotthermographie auf dreidimensionalen Oberflachen 


Ein Ziel dieser Arbeit ist die Bereitstellung von Methoden, die die Erfassung von Oberflä- 
chentemperaturen mittels Infrarotthermographie auf beliebig geformten Oberflächen ermög- 
licht. Nachdem im vorangegangenen Kapitel die notwendigen Grundlagen für die Modellierung 
einer beliebigen Kamera thematisiert wurden, werden im Folgenden die Korrekturverfahren 
vorgestellt. Hierbei ist vor allem der in Abschnitt 5.2 vorgestellte zweistufige Korrekturprozess 
maßgeblich für die erreichte Genauigkeit der Ergebnisse. 


5.1 Erweiterung des Schärfebereichs im Objektraum 


Der Bereich der noch akzeptierbaren Schärfe im Objektraum ist für das gegebene Kamerasys- 
tem gering. Dies ist auf die relativ große Eintrittspupille des Objektivs von f/2.3 sowie den ho- 
hen Abbildungsmaßstab ß durch Verwendung von Zwischenringen zurückzuführen. Ist es nun 
aufgrund der Krümmung der abzubildenden Oberflächen oder der Lage der optischen Zugänge 
nicht möglich, die abzubildende Oberfläche parallel zur Bildebene zu positionieren, kommt es 
zur unscharfen Abbildung von Teilen der Oberfläche. Bereiche mit hohen Strahlungsgradienten 
können dann nicht mehr sicher aufgelöst werden. 

Durch den Einsatz einer Aperturblende und der damit verbundenen Reduktion der Eintrittspu- 
pille Kann der Bereich der scharfen Abbildung erhöht werden. Je kleiner der Durchmesser der 
Aperturblende dap gewählt wird, desto größer wird der scharf abgebildete Bereich im Objekt- 
raum. Hier gilt es zwei Grenzen zu berücksichtigen. Durch Verkleinerung des Durchmessers 
dag wird das Strahlenbündel, welches in das Objektiv einfällt, radial beschnitten. Dies führt 
zu einer Reduktion der Intensität in der Bildebene und muss durch längere Integrationszeiten 
ausgeglichen werden. Außerdem haben bei kleinen Durchmessern dag die durch die Blenden- 
öffnung enstehenden Beugungsstrukturen negativen Einfluss auf die Bildschärfe. Durch die hö- 
here Wellenlänge der Infrarotstrahlung ist dies schon bei verhältnismäßig großen Durchmessern 
gegeben. Der Mindestdurchmesser der Aperturblende dag min wurde für den vorliegenden Fall 
experimentell ermittelt und eine Blende mit dag = 6mm eingesetzt. Die resultierende Erhöhung 
der Schärfentiefe wird in Abbildung 5.1(a) gezeigt. 


Sollte der Bereich der Schärfentiefe trotz Blende nicht ausreichen, müssen weitere Verfahren 
eingesetzt werden, die den scharfen Objektbereich erhöhen. Ein gängiges Verfahren aus der 
Makro-Fotografie ist die Fokus-Stapelung (engl. focus-stacking). Hierzu wird eine Bildreihe 
eines Objektes aufgenommen. Zwischen den Bildern wird die Schärfeebene in Richtung der 
optischen Achse durch Verschieben der Kamera oder Ändern des Auszuges verschoben. Wenn 
die Methode der Auszugsveränderung gewählt wird (diese ist praktisch einfacher umzusetzen), 
muss später beim Ausrichten der Bilder eine perspektivische Transformation erlaubt werden, 
da sich Öffnungswinkel und Perspektive der Kamera ändern. Bei Verschiebung des gesamten 
Kamerasystems ist eine Ähnlichkeitstransformation ausreichend, was die Rechenzeit reduziert. 
Nachfolgend werden die Bilder aneinander ausgerichtet, auf Schärfe analysiert und zusammen- 
gesetzt, wobei aus jedem Teilbild nur der scharf abgebildete Bereich verwendet wird. 
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Für die Anwendung auf Strahlungsaufnahmen wurde ein daran orientiertes Verfahren entwi- 
ckelt. Im Vergleich zu schon verfiigbaren Verfahren wird in dem hier vorgestellten keine Filte- 
rung der Messdaten durchgefiihrt. Das neue Verfahren wird im Folgenden vorgestellt. 


1. Bei einem stationären Betriebspunkt wird eine Infrarot-Aufnahmereihe des Objektes er- 
stellt. Die Einzelbilder können ohne Einschränkung der Allgemeinheit durch Verwendung 
der Verfahren von Ochs et al. (2010) bei verschiedenen Integrationszeiten t; aufgenom- 
men werden oder HDR-Bilder sein. Dabei wird die Schärfeebene in diskreten Teilschrit- 
ten (ca. 12) durch den gesamten interessierenden Objektbereich verschoben. 


2. Die Bilder werden aneinander ausgerichtet. Zur Vermeidung von Ausrichtungsfehlern 
durch Strömungsstrukturen werden die Kanten der Geometrie mithilfe eines gradienten- 
basierten Kantenfilters extrahiert. Bei gegebenem Dynamikumfang einer Aufnahme I; der 
Reihe wurde ein richtungsunabhängiger Sobelfilter (Sobel, 2014) mit den Kernel Sy, Sy 
1. Ordnung zur Berechnung der Gradienten G; als optimal ermittelt: 


Gyi = Sy * Íi, Gyi = Sy * Íi, (5.la) 


Gi= \/ Git O; (5.1b) 


Die eigentliche Ausrichtung der Bilder erfolgt dann mittels der ECC-Methode von Evan- 
gelidis et al. (2008). Zur weiteren Stabilisierung wird jeweils die Transformation zum 
nächsten Bild der Aufnahmereihe berechnet, da hier nur kleine Veränderungen vorliegen. 
Wird das mittlere Bild als (willkürliche) Referenz mit Pmit = 1 definiert, ergibt sich für 
das i-te Bild die Transformationsmatrix P; kom der Form 3 x 3 zu 


> 


A t 
V vo 


Pi kom = P; 2 P;-ı ee! Pit (5.2b) 


P= ; (5.2a) 


im Falle der perspektivischen Transformation mit den affinen Parametern A = [a1 ...aa], 
den perspektivischen Parametern V = [v, v2], dem Translationsvektor f = [t1 ,f2]™ und dem 
abhängigen Skalierungsfaktor vo. 


3. Jedes ausgerichtete Bild wird mittels eines Sobelfilters 2. Ordnung auf seine Scharfe ana- 
lysiert. Eine Schärfekarte mit gleichen Dimensionen wie das Bild wird erstellt. Diese 
enthält an jedem Pixel die Information, welches Bild der Aufnahmereihe die höchsten 
Gradienten an diesem Pixel aufweist. Diese Karte ist stark verrauscht. 


4. Die Schärfekarte wird extensiv gefiltert. Durch subsequente Anwendung eines Medianfil- 
ters! (OD, 3), eines Bilateralfilters (0,31, I = Xp = 138), eines Medianfilters (0,27) und 
schlussendlich eines Gaussschen Filters (LJ, 9) Kann eine glatte Karte erzeugt werden. Die 
eigentlichen Messdaten werden hierbei nicht verändert. 


5. Das Bild wird nach den in der Schärfekarte hinterlegten Informationen aus den Teilbildern 
zusammengesetzt. Eine weitere Filterung der Daten ist nicht notwendig. 


!Kernelform und Größe in Klammern 
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(a) Scharfentiefe einer Aufnahme mit Blende (b) Fokus-Stacking 


Abbildung 5.1: Vergleich der Schärfentiefe bei Anwendung von Fokus-Stacking 


Durch dieses Verfahren lassen sich beliebige Objekte scharf abbilden, solange deren Erstre- 
ckung im Fokusbereich der verwendeten Kamera liegt. Durch die extensive Filterung der Schär- 
fekarte wird ohne Filterung (und damit ohne Veränderung) der Messdaten ein glattes Bild er- 
zeugt. Ein beispielhaftes Ergebnis ist in Abbildung 5.1 dargestellt. Hierbei beträgt die Schär- 
fentiefe ca. 10mm, die Erstreckung des Objektes in Richtung der optischen Achse ca. 110mm. 
Das Ergebnis wurde aus 14 Bildern zusammengesetzt. 


5.2 Korrektur des Sichtwinkeleinflusses 


In Abschnitt 2.3 und Abschnitt 2.4 wurden die bisher eingesetzten Methoden zur Untersuchung 
verschiedener Geometrien mittels Infrarotthermographie beschrieben. Die bisher immer getrof- 
fene Annahme einer konstanten Emissivität soll im Folgenden genauer untersucht werden. Da 
die Wellenlängen- und Temperaturabhängiskeit der Emissivität der verwendeten Beschichtung 
im spektralen Band der verwendeten Kamera gezeigt wurde (Lohrengel et al., 1996; Dury et 
al., 2007), verbleibt die bisher vernachlässigte Sichtwinkelabhängigkeit. Dafür wird zuerst der 
Sichtwinkel relativ zur Kamera genauer definiert. Darauffolgend wird ein funktionaler Zusam- 
menhang zwischen Emissivität und Sichtwinkel vorgestellt. Da der Einfluss auf die Messge- 
nauigkeit gerade bei höheren Sichtwinkeln enorm ist, wird ein zweistufiges Kalibrierverfahren 
vorgestellt, das eine Reduktion des Messfehlers erlaubt. 


Wenn eine beliebige Oberfläche mit einer Kamera nach dem Lochkameramodell beobachtet 
wird, lassen sich zwei verschiedene Winkel definieren. Diese sind in Abbildung 5.2 dargestellt. 
Der Winkel zwischen dem Vektor vom Kamerazentrum zum lokalen Oberflächenpunkt und der 
lokalen Oberflächennormale wird im Folgenden als lokaler Oberflächenwinkel © bezeichnet. 
Der, vom Betrag kleinere, Winkel zwischen optischer Achse und lokaler Oberflächennormale 
wird als globaler Oberflachenwinkel 6 bezeichnet. Bei konsistenter Durchführung der Korrek- 
tur können theoretisch beide Winkel verwendet werden, praktisch ist die Wahl des globalen 
Winkels @ besser geeignet: Bei Kalibrierung des Kamerasystems (NUC) auf einer ebenen, iso- 
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gekrümmte Oberfläche 


Bildebene 


Abbildung 5.2: Definition der Oberflächensichtwinkel 


thermen Oberfläche wird garantiert, dass eine solche ebene Fläche mit 0 = 0° aber © £ 0° ein 
isometrisches Strahlungssignal produziert. Der Winkeleinfluss von © — 0 wird implizit korri- 
giert. Werden gekrümmte Oberflächen betrachtet, muss lediglich der hinzukommende Einfluss 
des globalen Winkels 6 korrigiert werden, diese Korrekturvorschrift kann unter Verwendung 
vorkalibrierter Strahlungsdaten außerhalb des Versuchsstandes erstellt werden. 


Zur Bestimmung dieser Vorschrift wird ein vereinfachter Aufbau verwendet. Dafür wird eine 
isotherme Oberfläche mit der im Versuch verwendeten Beschichtung bei verschiedenen Tempe- 
raturen mit dem Infrarotkamerasystem betrachtet. Diese Oberfläche beinhaltet zusätzlich einen 
schwarzen Strahler als Referenz. Der Versuchsaufbau wird in Abbildung 5.3(a) gezeigt. Bei 
verschiedenen Temperaturniveaus unter Variation des Sichtwinkels wird das Strahlungssignal 
der Oberfläche registriert. Bei Draufsicht @ = 0° wird zusätzlich die Schwarzkörperstrahlung 
und daraus abgeleitet die normale Emissivität €} = Es/E,„ bestimmt. So lassen sich die in Ab- 
bildung 5.3(b) gezeigten Datenpunkte ermitteln: Hier wird die mit der normalen Emissivität 
e€; normierte Emissivität der Oberfläche e = f(@) als Funktion des Sichtwinkels aufgetragen. 
Es zeigt sich ein stark nichtlinearer Abfall der Emissivität mit zunehmendem Sichtwinkel. 
Bis 0 ~ 50° ist keine nennenswerte Veränderung festzustellen, dies erklärt die Validität vie- 
ler der in Abschnitt 2.3 gezeigten Untersuchungen auf größtenteils ebenen Geometrien. Bei 
höheren Sichtwinkeln führt die Veränderung der Emissivität zu einer lokalen Veränderung des 
Strahlungssignals trotz homogener Oberflächentemperatur und muss zur Sicherstellung eines 
vertretbaren Messfehlers korrigiert werden. Zur Modellierung des gezeigten Zusammenhangs 
muss eine geeignete Funktion gewählt werden. Bei Oberflächenrauheiten von vergleichbarer 
Größenordnung wie die Wellenlänge schlägt Wen et al. (2006) eine bidirektionale Reflektanz- 
verteilungsfunktion BRDF für die Reflektivität r vor. Mit Gültigkeit des Kirchhoffschen Strah- 
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Abbildung 5.3: Betimmung der Winkelabhängigkeit der Emissivität 


lungsgesetzes € = 1 — r wird dann 


OSN Q ch, +6) (5.3) 


verwendet. Die modellierten Kurven können so die Messdaten hervorragend abbilden. Da die 
Funktion nicht auf Werte en < 1 beschränkt ist, kann im Bereich niedriger Sichtwinkel eine 
zu hohe Emissivität berechnet werden. Die Abweichungen sind gering, die betroffenen Punkte 
werden daher nachträglich auf eine maximale normierte Emissivität von en = 1 begrenzt. Eine 
weitere wichtige Feststellung ist der nicht vorhandene Einfluss der Oberflächentemperatur. Wie 
in Abschnitt 2.4 beschrieben, setzt sich die aufgenommene Strahlung aus Objektstrahlung und 
reflektierter Umgebungsstrahlung zusammen. Das Verhältnis dieser beiden Signale wird, neben 
dem Verhältnis von Oberflächentemperatur und der (Konstanten) Temperatur der Umgebung, 
von der Emissivität bzw. der Reflektivität beeinflusst. Wenn nun die Daten für verschiedene 
Oberflächentemperaturen zusammenfallen, kann der Betrag der Umgebungsstrahlung keinen 
relevanten Einfluss auf die ermittelten Korrekturfunktionen haben. Die Kalibrierung kann da- 
her unabhängig vom Betrag der Umgebungsstrahlung angewendet werden. 

Vergleichbare Daten wurden für die verwendete Beschichtung auch schon von Lohrengel et al. 
(1996) ermittelt. Während der dort gezeigte Verlauf ähnliche Form besitzt, ist die Absenkung 
der Emissivität bei höheren Winkeln moderater. Es ist zu beachten, dass neben der eigentli- 
chen Zusammensetzung der Beschichtung auch die Oberflächenbeschaffenheit einen wichtigen 
Einfluss hat. So wurden die Messungen von Lohrengel et al. bei niedriger Temperatur und ver- 
mutlich nicht gealtertem (siehe Abschnitt 2.4) Lack durchgeführt. Das Beschichtungsverfahren 
ist ähnlich, die angegebene Schichtdicke aber verhältnismäßig hoch. Die genaue Zusammenset- 
zung der Beschichtung, und ob eine Änderung derselben in gegebenem Zeitraum stattgefunden 
hat, ist außerdem nicht bekannt. Neuere Daten von laniro et al. (2010), die durch ein neues Ver- 
fahren ermittelt wurden, zeigen quantitativ vergleichbare Daten, geben jedoch keine konkreten 


68 Infrarotthermographie auf dreidimensionalen Oberflächen 


Informationen zur eingesetzten Beschichtung (matte black paint). Aufgrund dieser Unsicher- 
heiten werden die hier ermittelten Korrekturwerte verwendet. AuBerdem wird empfohlen, die 
Korrekturen nach jeder Neubeschichtung neu zu ermitteln und dabei die Korrekturoberfläche in 
der gleichen Charge wie die eigentlichen Priiflinge zu verarbeiten. Nur so kann ein Fehler durch 
abweichende Korrekturwerte, der durch die starke Nichtlinearität bei hohen Winkeln verstärkt 
wird, vermieden werden. 


Die Anwendung des Korrektors kann dann entweder in einem Schritt direkt bei der Temperatur- 
kalibrierung oder in zwei Schritten angewendet werden. Bei der Anwendung in zwei Schritten 
wird zuerst das Strahlungssignal korrigiert. Darauf folgend wird die Temperaturkalibrierung 
unter Berücksichtigung der veränderten Reflektivität durchgeführt. Dieses aufgeteilte Verfah- 
ren hat sich im Rahmen der durchgeführten Messungen als vorteilhaft erwiesen. So wird zum 
einen erreicht, dass eine isotherme Oberfläche im halbunendlichen Raum ein unabhängig vom 
Sichtwinkel homogenes Strahlungssignal produziert. Zum anderen ist die nötige Strahlungskor- 
rektur beim Durchführen der in-situ-Kalibrierung kleiner, was das Verfahren numerisch stabiler 
macht. Die Korrektur des Strahlungssignals ergibt sich direkt aus dem gemessenen Strahlungs- 
feld ID obj (X, 6) zu 


=> 


Tp obj,c(X) = Ip onj(X, 0) -en(0) " (5.4) 


mit der vorher bestimmten Korrekturfunktion ¢(@) und dem Sichtwinkel @(X), der sich aus der 
Richtung der optischen Achse OA (siehe Abschnitt 4.3) sowie dem lokalen Oberflachennorma- 
lenvektor zu 


(5.5) 


ergibt. 


5.3 Temperaturkalibrierung unter Berücksichtigung des Sichtwinkels 


Die Bestimmung der Temperatur-Strahlungs-Relation kann nach Korrektur des Sichtwinkelein- 
flusses ähnlich der von Ochs et al., gemäß Gleichung 2.13 präsentierten Methodik erfolgen. 
Nach Vorbestimmung der drei Parameter r,b, f verbleibt die vom Objekt reflektierte Strahlung 
loff als freier Parameter, der im eigentlichen Versuch bestimmt werden muss. Zur Bestimmung 
wird eine in-situ Kalibrierung anhand von diskreten Stützstellen durchgeführt. Hierbei muss 
pro Sichtbereich mindestens eine Stützstelle vorhanden sein. Zur Fehlerminimierung im Sin- 
ne eines LSO-Schätzers sind mehrere Stützstellen sinnvoll. Zusätzlich zu Gleichung 2.13 muss 
berücksichtigt werden, dass die reflektierte Strahlung sichtwinkelabhängig ist. Diese Abhän- 
gigkeit ergibt sich aus dem Kirchhoffschen Gesetz mit r = 1 — e, die reflektierte Strahlung œ r 
nimmt bei steigendem Sichtwinkel und sinkender Emissivität e zu. Bei gekühlten Bauteilen in 
heißen Umgebungen kann dies, entgegen der Intuition, zu einem Anstieg des unkorrigierten 
Strahlungssignals trotz fallender Emissivität führen. 

Daher muss das Kalibrierverfahren erweitert werden. Hierzu wird zuerst an jedem Stützpunkt i 
Gleichung 2.13 gelöst und so der Anteil der reflektierten Strahlung /orr,; bestimmt. Dieser kann 
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mit dem bekannten Korrekturfaktor en(0;) zu dem virtuellen Wert bei 6 = 0 skaliert werden: 
Ioft.i,c = lott, i 2 €n(6;) : (5.6) 


Die so bestimmten Wertepaare {T;,/orr,,.} werden dann mittels LSQ- Verfahren zur Bestimmung 
von Joff| go nach Gleichung 2.13 verwendet. Die Berechnung eines Temperaturfeldes aus ei- 
nem Strahlungsfeld erfolgt unter Beachtung der Identität an jeder Kalibrierstelle global durch 


> b 
T,(X,0) = : , (5.7) 
ern 


Eine weitere Herausforderung ergibt sich aus der Biotzahlähnlichkeit der Prüflinge. Bedingt 
durch die hohen Biotzahlen ergeben sich starke Temperaturgradienten über die zu untersuchen- 
de Wand. Unter Annahme stationärer ID-Wärmeleitung durch eine Wand der Dicke t ergibt 
sich dieser Gradient zu 


_ Biy- ATy 
an 


dy ($F) =hu-ATy bzw. (42) (5.8) 
W Ww 
Bei gegebenen Randbedingungen wird selbst in moderat gekühlten Bereichen ein Gradient von 
x~ 10K/mm erreicht. In stark gekühlten Bereichen ist ein weitaus höherer Gradient zu erwarten. 
Selbst bei geometrisch kleinen Thermoelementen ergibt sich eine hohe Messunsicherheit. Eine 
Reduktion dieser Messunsicherheit kann durch eine Reduktion der Biotzahl erreicht werden. 
Dies wird durch die Fertigung eines Kalibrierprüflings aus Kupfer erreicht. Bei gleicher Au- 
Bengeometrie (ohne Filmkühlung) und vergleichbaren aerodynamischen Bedingungen kann ei- 
ne Reduktion der Biotzahl und damit des Temperaturgradienten in der Wand erreicht werden. 
Mit Bin,riau/Bip,cu © 50 wird der Temperaturgradient deutlich reduziert. Der Kupferprüfling kann 
zusätzlich innen gekühlt werden und erlaubt so das Einstellen von verschiedenen Oberflächen- 
temperaturen. In diesen Prüfling werden an der Außenfläche Nuten eingebracht, in die Thermo- 
elemente geklebt werden (zehn Thermoelemente Typ K mit je 0,25 mm Durchmesser), die als 
Stützstellen dienen. Letztlich wird der Prüfling mit Nextel Velvet Coating beschichtet. 
Die Kalibrierung wird dann anhand der Temperatur-Strahlungs-Wertepaare auf der Kupfer- 
schaufel nach obigem Verfahren für jede Kameraposition bestimmt und später auf die Strah- 
lungsmessung des eigentlichen Prüflings angewandt. Hierbei wird die Annahme getroffen, dass 
die so bestimmte reflektierte Strahlung Zoff unabhängig von der eigentlichen Temperatur des 
Prüflings ist (kein thermodynamischer Strahlungsleistungsfluss zwischen Kanalwänden und 
Prüfling für den am Prüfling reflektierten Strahlungsanteil). So muss während der Kalibrierung 
nur sichergestellt werden, dass die den reflektierten Strahlungsanteil maßgeblich beeinflussende 
Wandtemperatur des Strömungskanals konstant über alle Versuche ist. Dies wird durch Einhal- 
tung des Betriebspunktes realisiert. Zusätzlich können durch Variation der internen Kühlung 
und damit der Wandtemperatur des Kupferprüflings aus einem Thermoelement mehrere Wer- 
tepaare generiert werden, was die Qualität der Kalibrierfunktion erhöht und einen Extrapolati- 
onsfehler bei der späteren Messung verringert. 
Die finalen Kalibrierkurven für verschiedene Positionen sind in Abbildung 5.4 gezeigt. Hierbei 
wurden exemplarisch zwei Sichtbereiche ausgewählt, von denen einer im Mittel hohe Sichtwin- 
kel auf den Prüfling hat ((0z3)) und der andere im Mittel geringe Sichtwinkel aufweist ((0z2)). 
Bei beiden Sichtbereichen zeigt sich eine nur geringe Restabweichung zwischen der verscho- 
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(a) Finale Kalibrierkurve, hohe Sichtwinkel, (0z3) (b) Finale Kalibrierkurve, niedrige Sichtwinkel, 
(0z2) 


Abbildung 5.4: Finale Kalibrierkurven fiir verschiedene Sichtpositionen. Vorkalibrierung un- 
terschiedlich durch Doppelscheibe bei (022) 


benen Kalibrierfunktion und den Stiitzstellen. Das Verfahren der in-situ Kalibrierung kann also 
auch bei komplexen Geometrien angewendet werden. Die Stiitzpunkte fiir die Sichtposition in 
Abbildung 5.4(b) verschieben sich durch den neu eingeführten Reflektivitätskorrektor praktisch 
nicht (Sichtwinkel 0 im Mittel < 50°). Der RMS Fehler ist niedrig und verändert sich nicht wei- 
ter durch die Korrektur. Im Gegensatz dazu steht der in Abbildung 5.4(a) gezeigte Sichtbereich 
mit hohem Sichtwinkel (9 bis zu 85°). Auch hier ist eine gute Übereinstimmung festzustellen. 
Der RMS Fehler kann durch die Korrektur weiter reduziert werden. Der verbleibende Fehler 
ist im Vergleich zu geometrisch einfachen Sichtpositionen trotz Korrektur höher. Dies ist auf 
die Restunsicherheit in der Schätzung der Kameraposition sowie der Güte der Korrekturfunk- 
tion nach Gleichung 5.5 zurückzuführen. Die starke Nichtlinearität der Funktion führt schon 
bei kleinen Abweichungen zu betragsmäßig großen Fehlern, wenn diese Abweichungen bei ho- 
hen Sichtwinkeln auftreten. Für die später durchgeführte Fehlerabschätzung wird daher eine 
Unsicherheit bei hohen Sichtwinkeln betrachtet. 


5.4 Fehleranalyse 


In den vorangegangenen Abschnitten wurden die nötigen neuen Verfahren beschrieben, die die 
Bestimmung von Oberflächentemperaturen auf komplexen Oberflächen erlauben. Abschließend 
soll der Fehler abgeschätzt werden, der nach Anwendung der neuen Verfahren verbleibt. Tech- 
nisch relevant ist die Abweichung der gemessenen Oberflächentemperatur zur wahren Tem- 
peratur der Oberfläche, formuliert als Unsicherheit der Totalkühleffektivität. Nur wenn diese 
Unsicherheit hinreichend gering ist, Können die Daten zur Prognose der Effektivität der Kühl- 
systeme herangezogen werden. Zur Quantifizierung dieser Abweichung wird eine konservative 
Abschätzung durchgeführt. An einem Punkt, der unter hohem Sichtwinkel von 0 = 75° betrach- 
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tet wird, kann zunächst der Fehler in der Emissivität Ag bestimmt werden: 


CEN 
00 
Ohne Korrektur Unc ist dies genau die Differenz zur korrigierten Emissivität. Mit Korrektur 
c wird das Differential bestimmt. Die Winkelabweichung d@ ergibt sich aus der in Unter- 
abschnitt 4.3.2 genannten Grenzabweichung von Az | < 4mm bezogen auf einen minimalen 
Objektabstand von 200 mm zu 1,15°. Angenommen werden daher 1,5°. Der Fehler im Tempe- 
raturbereich berechnet sich dann zu 

oT, oT, oT, 


= -Ah -AI — -A 5.10 


Aennc=1-e(O)n|,_ 7. bzw. Aenc = dal eens (5.9) 


AT; 


mit der oben bestimmten Abweichung der normierten Emissivität A ey, der Abweichung der in- 
situ bestimmten reflektierten Strahlung A off = 6% (aus Kalibriermessungen bestimmt) und der 
Unsicherheit der Strahlungsmessung Alp. Diese ergibt sich zu einer konstanten Restunsicher- 
heit des Detektors von 0,5 % nach Ochs et al. (2009). Zu beachten ist, dass die Detektorstrahlung 
Ip nach Gleichung 5.4 auch von der normierten Emissivität ey abhängt. 
Letztlich bleibt noch ein konstanter Fehler durch die Messung mit Thermoelementen. Dieser 
wird mit A Trc = 0,75K abgeschätzt. Dieser Fehler ist sowohl bei den Infrarotmessungen als 
auch bei der Bestimmung der Randbedingungen gegeben. Final kann die fehlerbehaftete Total- 
kühleffektivität berechnet werden: 

TCE(X) = WERU mit 7;(X) = R(X) +(AT, + ATzc). (5.11) 

TH — Tk 

Das numerische Ergebnis dieser Abschätzung bei gegebenen Parametern der Kalibrierfunktio- 
nen wird in Abbildung 5.5 gezeigt. Auf Basis der für einen typischen Temperaturbereich berech- 
neten virtuellen Einstrahlungen wird Gleichung 5.10 gelöst. Mit gegebenen Randbedingungen 
nach Tabelle 3.1 kann durch Gleichung 5.11 der Fehler bestimmt werden. Auf der linken Achse 
des Diagramms ist die Totalkühleffektivität abgetragen, auf der rechten Achse deren prozen- 
tualer Fehler. Die Darstellung erfolgt über einen prüfstandstypischen Bereich von Oberflächen- 
temperaturen. Der Zusammenhang der nominellen Totalkühleffektivität ist per Definition linear. 
Die im Betrieb zu erwartenden Grenzwerte der Totalkühleffektivität sind als senkrechte, unter- 
brochene Linien eingetragen. Der berechnete Fehler ist als Band um die nominelle Totalkühlef- 
fektivität absolut (linke y-Achse) und relativ (rechte y-Achse) auf den jeweiligen Nominalwert 
bezogen abgebildet. Zusätzlich ist der mittlere Fehler im eingezeichneten, maschinentypischen 
Bereich der Kühleffektivität Øeng angegeben. Durch Anwendung der Korrektur kann der Fehler 
maßgeblich reduziert werden. Der mittlere Fehler sinkt auf Ong © 2,3 %, obwohl die Restun- 
sicherheiten verhältnismäßig groß angenommen wurden. Diese Restunsicherheiten wirken sich 
vor allem auf die korrigierten Daten aus. Von einer weiteren Reduktion des Restfehlers, spezi- 
ell bei Anwendung der Korrekturverfahren, ist auszugehen. Dies kann auch bei den in Abbil- 
dung 5.4 gezeigten Kalibrierkurven festgestellt werden: Selbst bei hohem Sichtwinkel liegt der 
verbleibende RMS-Fehler bei 4,266 K - dies entspricht nur etwa 1 % der gemittelten Oberflä- 
chentemperatur. 
Weiterhin soll quantifiziert werden, ab wann die Anwendung der Korrekturverfahren sinnvoll 
ist. Dies hängt allein von den zu erwartenden Sichtwinkeln ab. Die Abschätzung des Fehlers 
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Abbildung 5.5: Fehlerabschatzung der Kalibriermethodik bei gegebenen Randbedingungen fiir 
einen hohen Sichtwinkel 0 = 75° 


wird bei gleichen Annahmen für einen Sichtwinkelbereich von 0 < 0 < 85° durchgeführt, es 
wird jeweils der mittlere Fehler Øeng bestimmt. Das Ergebnis dieser Untersuchung wird in Ab- 
bildung 5.6 dargestellt. Hier zeigt sich der deutlich spätere Anstieg des Messfehlers bei Ver- 
wendung der Korrektur, der Messbereich kann erweitert werden. Wird beispielsweise davon 
ausgegangen, dass schon eine Änderung der Metalltemperatur um 20 K einen deutlichen Ein- 
fluss auf die Lebensdauer hat (vgl. Kapitel 1), Kann so eine nötige minimale Unsicherheit der 
Kühleffektivität berechnet werden. Hierbei muss auf Basis der linearen Abschätzung mittels 
Differentation der Fehler auf die maximale Temperaturdifferenz Ty — Tx skaliert werden. Dar- 
aus ergibt sich unter Annahme der Randbedingungen der Maschine eine tolerierbare Unsicher- 
heit von 2,5 % bei einer approximierten Oberflächentemperatur von 900 K. Dieser tolerierbare 
Fehler ist in Abbildung 5.6 markiert. Ohne Korrektur ergibt sich ein maximaler geometrischer 
Messbereich von 6 = 49°, mit Korrektur eine Erweiterung des Messbereiches bis zu einem 
Sichtwinkel von 6 = 76°. Aufgrund der deutlich konservativen Abschätzung werden daher für 
die Praxis die vorher beschriebenen Grenzen von 0 = 50° sowie 6 = 80° empfohlen. 


Restfehler bei Anwendung der in-situ Kalibrierung 


Ein bisher noch nicht beachteter Fehler ergibt sich bei der in-situ Kalibrierung mit in den se- 
paraten Kalibrierprüfling eingelassenen Thermoelementen. Der Aufbau des Kalibrierkörpers ist 
maßstabsgetreu in Abbildung 5.7 gezeigt. Dieser muss ebenfalls mit Nextel Velvet Coating be- 
schichtet werden. Aufgrund sehr geringer Wärmeleitfähigkeit der Beschichtung ergeben sich 
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Abbildung 5.6: Fehlerabschätzung der Kalibriermethodik bei Veränderung des Sichtwinkels 


trotz geringer Schichtdicke hohe Temperaturgradienten normal zur Oberfläche (Daten zur Wär- 
meleitfähigkeit nach Lohrengel et al. (1996), Schichtdicke an Proben bestimmt zu 0,45 um). 
Daher erfassen die unter der Lackschicht liegenden Thermoelemente nicht die eigentliche Ober- 
flächentemperatur. Dies wurde durch die Annahme eines ID-Wärmestroms und dem daraus re- 
sultierenden Temperaturabfall ausgeglichen. Dieser kann im Versuch durch die Wärmebilanz 
der Kalibrierschaufel flächengemittelt bestimmt werden. Aus dem dadurch geschätzten spezi- 
fischen Wärmestrom lässt sich der Temperaturabfall über die Lackschicht ermitteln (Verfahren 
nach Ochs, 2011). 

Eine weitere Herausforderung ergibt sich aus dem hohen lateralen Wärmewiderstand der Be- 
schichtung. So ist davon auszugehen, dass hohe thermische Gradienten auf der Oberfläche von 
den Thermoelementen, eingebettet im gut leitenden Substrat, nicht abgebildet werden können. 
Aufgrund der diffusiven Natur der Wärmeleitungsgleichung ist eine Rekonstruktion nicht mög- 
lich. An dieser Stelle muss in Zukunft überprüft werden, ob durch eine Optimierung der Wärme- 
widerstände von Substrat und Beschichtung, auch unter Berücksichtigung der Temperaturgradi- 
enten in Bezug auf den Durchmesser der Thermoelemente, eine Verringerung dieser Unsicher- 
heit möglich ist. Dieser Kalibrierfehler verschiebt das mittlere Niveau des Temperaturfeldes. 
Die auf der Schicht abgebildeten Gradienten können trotzdem auf den Einsatz in der Maschine 
übertragen werden. Zum einen ist die Biotzahl des Prüflings (und die der Schaufel in Maschi- 
nenumgebung) deutlich höher als die des Kupferkörpers, zum anderen werden die Schaufeln 
vor dem Einsatz meist mit den in Abschnitt 2.2 beschriebenen thermischen Dämmschichten be- 
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Abbildung 5.7: Temperaturdämpfungsverhalten von Lackschicht und stark konduktivem Sub- 
strat 


schichtet. Diese haben ihrer Aufgabe geschuldet sehr geringe Wärmeleitfähigkeiten. Zusätzlich 
werden diese Schichten oft mit wandnormalen Rissen aufgebracht, was eine laterale Wärmelei- 
tung weiter reduziert. 


Durch die in diesem Kapitel vorgestellten Korrekturverfahren kann nun erstmals die Tempera- 
turverteilung auf beliebig geformten Oberflächen mittels Infrarotthermographie bestimmt wer- 
den. In einem letzten Schritt müssen dann die aus mehreren diskreten Sichtbereichen aufgenom- 
menen Daten auf der Prüflingsgeometrie zusammengesetzt werden. Dafür werden die Tempe- 
raturfelder zunächst durch Berechnung der Totalkühleffektivität nach Unterabschnitt 2.3.2 mit 
den Totaltemperaturen von Heißgas und Kühlluft entdimensioniert. Die verwendeten Tempera- 
turen werden zeitaufgelöst gemessen und zum jeweiligen Zeitpunkt der Bildaufnahme ausge- 
wertet. Dies erlaubt einerseits den Vergleich mit anderen Daten, andererseits werden eventuelle 
Schwankungen der Randbedingungen und deren Effekt auf die Oberflächentemperatur redu- 
ziert. 

Falls Datenbereiche in Objektkoordinaten überlappen, wie beispielsweise zwischen (ozl) und 
(0z3), müssen diese Bereiche überblendet werden. Dafür wurde ein neues Verfahren entwickelt, 
welches in Abschnitt A.3 erläutert wird. Das finale Ergebnis ist eine in 3D-Objektkoordinaten 
vorliegende Verteilung der Totalkühleffektivität. 


6 Ergebnisse der experimentellen Untersuchungen 


Durch die in den vorangegangenen Kapiteln eingeführten neuen Verfahren können Oberflä- 
chentemperaturen mit hoher Genauigkeit auf den stark gekrümmten und gekühlten Prüflingen 
bestimmt werden. Die Ergebnisse dieser Untersuchungen werden im Folgenden vorgestellt. 
Das Ziel der Untersuchungen ist die Beurteilung der neuen auf Drallkühlung basierenden Kon- 
zepte zur Erhöhung des inneren Wärmeübergangs und damit zur Absenkung der Metalltempe- 
ratur. Beurteilt werden die Kühlkonzepte anhand der aerothermischen Randbedingungen sowie 
der Oberflächentemperatur auf der Heißgasseite. Zur besseren Vergleichbarkeit der absoluten 
Werte wird zusätzlich eine konventionelle Innenkühlung mit den gleichen Verfahren bewertet. 
So kann durch eine zusätzliche relative Analyse die Restunsicherheit weiter reduziert werden. In 
einem ersten Schritt wird die qualitative Verteilung der Totalkühleffektivität auf der Außenflä- 
che der Turbinenrotorschaufel anhand der konventionellen Geometrie untersucht. Diese Vertei- 
lung ist für alle Geometrien ähnlich. Anschließend wird die konventionelle Referenzgeometrie 
quantitativ bewertet. Auf Basis dieser Daten können die neuen internen drallbasierten Geome- 
trien untersucht werden. Bei den neuen Geometrien wird nur der Einsatz zur Drallerzeugung 
variiert. Es kommen die in Abschnitt 3.2 gezeigten Einsätze zur Verwendung. Hieraus ergeben 
sich sechs Konfigurationen: Die konventionelle Referenzgeometrie (CV) und fünf Drallgeome- 
trien (GVI), (NSG), (OSG0), (OSG3), (OSG8). 

Hierbei dient das Eintrittstotaldruckverhältnis der Kühlluft II (Gleichung 6.1) bei ansonsten 
konstanten Größen zur Klassifizierung des Aufwandes analog zu fk/fxo. Die aus der Oberflä- 
chentemperatur abgeleitete Totalkühleffektivität TCE soll bei konstanten Heißgasbedingungen 
und damit konstantem Heißgaswärmeübergang hy als Analogie zu Nux/Nugo optimiert werden. 
Neben dem Ziel einer möglichst hohen Totalkühleffektivität muss auch die Verteilung dieser 
berücksichtigt werden. So sollen lokale Abweichungen mit starken Gradienten (hot spots oder 
cold spots) vermieden werden. 

Der Betriebspunkt für die Heißgasströmung wird konstant gehalten. Für die Kühlluftströmun- 
gen werden drei Betriebspunkte analysiert. Der Nennbetriebspunkt wird nach Abschnitt 3.1 
bei ioo = 100% - & festgelegt. Davon ausgehend wird ein Betriebspunkt mit hohem Kühl- 
luftmassentrom von &120 = 120% - & und ein Betriebspunkt mit niedrigem Massenstrom von 
Son = 90% - & untersucht. Durch Variation des Massenstroms wird der Eintrittstotaldruck der 
Kühlluft variiert. Hierbei wird angenommen, dass ein Großteil des für die auf Drallkühlung 
basierenden Konzepte notwendigen Überdrucks im Drallerzeuger abgebaut wird, die Dichte 
im Kühlkanal selbst ist näherungsweise gleich! , ein Effekt auf den Wärmeübergang durch 
eine Änderung der statischen Dichte aufgrund des erhöhten Vordruckes wird vernachlässigt. 
Die Eintrittstotaltemperatur und der split, also die Massenaufteilung auf Filmkühlung und die 
beiden Abnahmen an der Schaufelspitze purge I (¢)1) sowie purge 2 (Cp2), werden konstant 
gehalten. So wird sicher gestellt, dass der gleiche Massenstrom zur Filmkühlung zur Verfügung 
steht. Lokale Unterschiede in der Filmwirkung sind dann auf lokal veränderte Anströmwinkel 
B der internen Strömung zurückzuführen. 


! Aufgrund von unterschiedlicher Geometrie der Kanäle und Verrippung im konventionellen Fall wird der Druck- 
abbau gradueller verlaufen. Am Ende des bzw. der Vorderkantenkanäle ist der Druck dann, bedingt durch die 
konstanten purge-Massenstromverhältnisse angeglichen. Im hinteren Schaufelbereich sind beide internen Geo- 
metrien gleich. 
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Tabelle 6.1: Thermodynamische Randbedingungen der Untersuchungen 


T II Cpt / Spe 
Coo C100 C120 Coo So 120 Lo ioo S120 
Konventionell (CV) 0,530 0,529 1,09 1,16 1,28 17/10 16/10 


Drall (GVI) 0,530 0,529 1,07 1,12 1,25 17/10 15/9 


Drall (NSG) © 0,529 0,528 0,527 1,13 1,21 1,36 16/10 15/10 


Drall(OSGO) __. 0,529 0,528 1,12 1,19 1,34 17/10 16/10 
Drall (OSG3) 0,530 0,529 0,528 1,15 1,23 1,39 16/10 
Drall(OSG8) ~~ 0,528 0,527 1,17 1,26 1,44 16/10 15/10 


Die Randbedingungen fiir alle durchgeführten Versuche sind in Tabelle 6.1 gegeben. Hier- 
bei beschreiben I’, II das Temperatur- sowie das Druckverhältnis und ¢p1, ¢p2 die purge- 
Massenstromverhältnisse nach 

MK.22 


1) = Gern mit PH,LE = 0,97 - PH,tot,1 (6.1) 


Tk,1 PK.tot,1 MK21 
= I= Cpl = > 
Ta,ı PHLE &- mH 


mit den direkt stromab der Schaufelspitze (Mx 2;, Abbildung 3.1) gemessenen purge- 
Massenströmen mg/zi. Die Abweichungen der Randbedingungen während dem Testzeitraum 
einer bestimmten Konfiguration an einem bestimmten Betriebspunkt lagen unter 0,5 %rel 
für Temperatur- und Druckverhältnis. Die Schwankung der Massenstromaufteilung konnte 
unter 1,5 %reı gehalten werden. Bei pı = 16% Nennwert und Nennmassenstrom der Kühlluft 
entspricht diese Abweichung ~% 0,05g/s. Diese Abweichungen werden daher nicht explizit 
angegeben. Weiterhin wird die Massenstromaufteilung auf eine Ganzzahl gerundet. Der 
maximale Rundungsfehler beträgt 0,1 g/s. 


6.1 Analyse der qualitativen Temperaturverteilung 


Zuerst wird die qualitative Verteilung der Totalkühleffektivität anhand der konventionellen Geo- 
metrie untersucht. Die gezeigten Strukturen sind maßgeblich von der Heißgasströmung aufge- 
prägt und für alle Konfigurationen ähnlich. Durch eine Zuordnung zu charakteristischen Strö- 
mungsstrukturen aus den Ölvisualisierungen können die Effekte beurteilt werden. Im Folgenden 
werden die finalen Kühleffektivitäten der konventionellen Geometrie am Nennbetriebspunkt 
C100 in Falschfarbendarstellung mit den Konturlinien der Ölfarben auf der 3D Geometrie über- 
lagert. Die Konturwerte der Öllinien haben keine physikalische Bedeutung und werden daher 
nicht explizit erwähnt. Sprünge an den Bildgrenzen der Ölvisualisierung können nicht vermie- 
den werden, da leichte Belichtungsunterschiede zwischen den einzelnen Bildern vorliegen. Zur 
besseren Orientierung sind die drei wichtigsten Ansichten in Abbildung 6.3 eingetragen. 

Anhand der niedrigen Kühleffektivität entlang der Vorderkante, dargestellt in Abbildung 6.1(a) 
rechts bzw. im Detail in Abbildung 6.1(b), lässt sich die Lage der Staulinie bestimmen. Auf- 
grund der minimalen Grenzschichtdicke und dem daraus resultierenden hohen Wärmeüber- 
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(a) Verteilung der Kühleffektivität und Konturen (b) Verteilung der Kühleffektivität, Detailansicht 
der Ölvisualisierung. Ansicht Front Staupunkt 


Abbildung 6.1: Verteilung der Kühleffektivität und Konturen der Ölvisualisierung. Vorderer 
Schaufelbereich. Konventionelle Geometrie, Betriebspunkt {00 


gangskoeffizienten wird hier das globale Minimum der Kühleffektivität erreicht. Die Staulinie 
liegt entlang der gesamten Schaufelhöhe auf der geometrisch definierten Vorderkante. Die kor- 
rekte Anströmung des Prüflings kann so nochmals bestätigt werden. 

Stromab der Vorderkante, entlang der Saugseite, ist ein in Strömungsrichtung schneller Anstieg 
der Kühleffektivität zu beobachten. Durch Ausbildung der Grenzschicht nimmt der heißgasseiti- 
ge Wärmeübergangskoeffizient ab, die Wandstärke wird reduziert und die Wärmeabfuhr an den 
Bohrungsinnenflächen der folgenden, stark angestellten (siehe auch Abbildung 3.2(a)), saugsei- 
tigen Filmkühlreihe nimmt zu. Anhand der Ölvisualisierung ist eine Aufstauung der Strömung 
an den saugseitigen Filmkühlbohrungen zu erkennen, gefolgt von einer Beschleunigung nahe 
der jeweiligen Endwände. 

Der Bereich stromab der Filmkühlungsreihe wird in Abbildung 6.2(a) gezeigt. Der mittlere 
Schaufelbereich von = 15% bis ~ 85% Schaufelhöhe ist durch die stromauf liegende Film- 
kühlreihe effektiv gekühlt. Die verwendete Diffusorgeometrie begünstigt die hohe laterale Ab- 
deckung bei gleichmäßiger Effektivitätsverteilung sowie die in Strömungsrichtung langsame, 
graduelle Ausmischung des Kühlfilms. Die erreichte Totalkühleffektivität wird hier maßgeblich 
durch die Filmkühleffektivität beeinflusst. Ein deutlicher Unterschied der Totalkühleffektivi- 
tät ergibt sich stromab der Bereiche ohne Filmkühlbohrungen nahe der Endwände. Die in der 
Ölvisualisierung sichtbare Staulinie des nabenseitigen Passagenwirbels kann mit einem Be- 
reich niedriger Totalkühleffektivität verbunden werden. Der Passagenwirbel mischt vermutlich 
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(a) Verteilung der Kühleffektivität und Konturen (b) Verteilung der Kühleffektivität und Konturen 
der Ölvisualisierung. Ansicht Saugseite der Ölvisualisierung. Ansicht Druckseite 


Abbildung 6.2: Verteilung der Kühleffektivität und Konturen der Ölvisualisierung. Konventio- 
nelle Geometrie, Betriebspunkt C199 


Heißgas in die Grenzschicht. Dieser Bereich kann bis zur Hinterkante verfolgt werden. An der 
Schaufelspitze ist dagegen ein Bereich hoher Künhleffektivität zu erkennen. Dieser lässt sich mit 
dem in Abschnitt 3.3 beschriebenen Ablösegebiet erklären, das durch die vorher beschriebene 
Aufstauung verstärkt wird. Dieses wirkt als Pufferzone mit geringer Mischung in Wandnorma- 
lenrichtung, der Heißgaswärmeübergang wird reduziert. Der kleinere, endwandseitige Hufei- 
senwirbel wandert erst stromab dieses Gebietes an der Schaufel in Richtung der Nabe. Während 
diese Grenzlinie durch die Ölvisualisierung deutlich zu erkennen ist, kann dieser in den Daten 
der Totalkühleffektivität keine klar zu erkennende Struktur zugeordnet werden. 

Kurz stromauf der Hinterkante sind mehrere Sprünge der Kühleffektivität zu erkennen. Diese 
sind auf eine inhomogene Strahlung der Scheibenrahmen zurückzuführen und damit an keine 
Änderung der Oberflächentemperatur bzw. der Kühleffektivität gekoppelt. 

Die Darstellung stromab der Vorderkante entlang der Druckseite wird in Abbildung 6.2(b) ge- 
zeigt. Durch die Anstellung der Filmkühlbohrungen am Staupunkt ergibt sich eine stark in- 
homogene Verteilung der Totalkühleffektivität stromab der Vorderkante. Im Mittel steigt diese 
aber wie auch entlang der Saugseite an. Weiterhin ist die schon in Abschnitt 3.3 beschriebene 
Umlenkung der Strömung zu niedrigeren Radien zu beobachten. Die Strukturen der Kühleffek- 
tivität bestätigen die Prognose der Ölvisualisierung. 

Durch die hohe Auflösung des Messsystems können nun auch Details wie die Struktur der ein- 
zelnen Filmkühlstrahlen dargestellt werden. Abbildung 6.1(b) zeigt diese Strukturen im Bereich 
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der nabennahen Vorderkante in starker Vergrößerung. Durch die starke Anstellung der Bohrung- 
en sowohl bezüglich der Heißgasströmung als auch der Kühlluftströmung aus fertigungstech- 
nischen Gründen weicht die Kühleffektivität stark von den Prognosen aus bisherigen Untersu- 
chungen ab (zusammengefasst in Unterabschnitt 2.2.2). Direkt stromab des Bohrungsaustrittes 
ist eine niedrige Kühleffektivität zu beobachten. Das Heißgas (Hauptströmungsrichtung in ge- 
gebener Darstellung näherungsweise horizontal) trifft auf den austretenden Kühlluftstrahl und 
bildet ein Wirbelsystem, das den Wärmeübergangskoeffizienten erhöht. Der eigentliche Kühl- 
film bildet sich oberhalb der Bohrung und überdeckt nur einen schmalen, scharf abgegrenzten 
Bereich (in Abbildung 6.1(b) schraffiert). Darüber hinaus lässt sich eine Aufspaltung des Kühl- 
films erkennen, die sich vor allem bei den nabennahen Bohrungen der Filmkühlreihen an der 
Vorderkante ausbildet. Der hier höhere Totaldruck der Kühlluft führt zu einer Erhöhung des 
Massenstroms durch diese Bohrungen und damit zu einer stärkeren Ausprägung des Nieren- 
wirbelpaares. Bei gegebenem Betriebspunkt der konventionellen Geometrie führt die Ausbil- 
dung des hoch intensiven Wirbelpaars zur Aufspaltung des Kühlluftstrahls beim Austritt aus 
der Bohrung. Das Heißgas Kann dann zwischen den Teilstrahlen eindringen. Die Wirbelstruktur 
und deren Veränderung bei verschiedenen Betriebspunkten konnte von Stichling (2018) gezeigt 
werden, die Aufspaltung des Strahles wurde in den durchgeführten numerischen Untersuchun- 
gen jedoch nicht beobachtet. 


Die im vorangegangenen Abschnitt gezeigten Strukturen der Totalkühleffektivität konnten den 
wichtigsten Strukturen der Sekundärströmung im Heißgaskanal zugeordnet werden. Diese sind, 
zumindest qualitativ, bei allen untersuchten Geometrien zu finden. Ein quantitativer Vergleich 
verschiedener Konzepte kann daher anhand von lokal gemittelten Daten durchgeführt werden. 


6.2 Untersuchung der Referenzgeometrie 


Zusätzlich zur qualitativen Untersuchung des Temperaturfeldes werden die Daten im Folgenden 
auch quantitativ beurteilt. Aufgrund der erreichten Ähnlichkeit zwischen Versuch und Maschi- 
ne können die so bestimmten Kühleffektivitäten direkt zur Bestimmung von absoluten Ober- 
flächentemperaturen von Bauteilen in Maschinenumgebung verwendet werden. Die ermittel- 
ten Kühleffektivitäten werden an mehreren diskreten Positionen über die Schaufelhöhe unter 
Berücksichtigung der Fläche (bei gleicher Erstreckung des Einzugsbereiches in Richtung der 
Laufkoordinate also der Höhe des Bereiches) gemittelt und über die gemittelte Position auf 
der Prüflingsoberfläche aufgetragen. Die Positionen der Mittelung sind in Abbildung 6.3 als 
Diamanten dargestellt. Diese Positionen werden anhand einer Parametrisierung der Lauflänge 
entlang der Saugseite (positives Vorzeichen, rechts des Staupunktes) und entlang der Druckseite 
(negatives Vorzeichen, links des Staupunktes) entdimensioniert und über diese Laufkoordinate 
aufgetragen. Der Staupunkt ist als senkrechte, gepunktete Linie markiert. Die lineare Interpola- 
tion dient nur zur Visualisierung des Trends. 


Die über die Schaufelhöhe gemittelten Kühleffektivitäten für die konventionelle Geometrie bei 
drei Betriebspunkten sind in Abbildung 6.4 dargestellt. Die minimale Totalkühleffektivität TCE 
wird, wie schon aus den dreidimensionalen Daten in Abschnitt 6.1 ersichtlich, am Staupunkt er- 
reicht. Eine Steigerung des Kühlluftmassenstroms führt hier zu einer annähernd linearen Stei- 
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Abbildung 6.3: Mittelschnitt Seitenverhältnis 1 : 1 mit Punkten der Mittelung über die Schau- 
felhöhe und Positionen der Filmkühlbohrung und Blickrichtungen fiir die 3D 
Daten. Druckseite schwarz, Saugseite rot. Kühlluftreihen *, Auswertepositio- 
nen 


gerung der Kühleffektivität. Im Mittel wird eine Kühleffektivität von TCE ~ 42,5 % erreicht, 
was dem in Unterabschnitt 2.3.2 genannten Minimum entspricht. Da diese Daten von der schon 
im Betrieb erprobten Geometrie stammen, können so die Messmethodik und die Genauigkeit 
der quantitativen Ergebnisse nochmals bestätigt werden. 

In Hauptströmungsrichtung entlang der Druckseite kann zuerst ein schneller Anstieg der To- 
talkühleffektivität beobachtet werden. Die leicht stromab vom Staupunkt liegenden Kühlluft- 
bohrungen senken die effektive Temperatur für den Wärmeübergang auf der Heißgasseite. Die 
Grenzschicht auf der Heißgasseite dickt auf und reduziert so den Wärmeübergangskoeffizienten 
hx. Außerdem wird hier der Bereich der inneren Verrippung erreicht. Letztlich sinkt die Wand- 
stärke. Weiter stromab bleibt die Kühleffektivität annähernd konstant, in diesem Bereich über- 
lagern sich mehrere Effekte. Zwar wird durch eine weitere Reihe Filmkühlbohrungen zwischen 
P3 und P4 ein neuer Kühlfilm erzeugt und so die Filmkühleffektivität wieder erhöht, jedoch 
steigt die Temperatur der Kühlluft am Übergang zwischen Vorderkantenkanal und erstem rück- 
läufigen Multipass an. 

Entlang der Saugseite stromab der Vorderkante steigt die Kühleffektivität ebenfalls schnell an. 
Dies ist auf vergleichbare Effekte wie auf der Druckseite und den weitestgehend symmetri- 
schen Aufbau der vorderen Kühlkanäle zurückzuführen. Die maximale Kühleffektivität wird 
im Bereich der saugseitigen fanshaped-Bohrungen erreicht (S2 ca. zwei Diffusordurchmesser 
stromauf, S3 ca. zwei Diffusordurchmesser stromab der jeweiligen Diffusorkanten). 


6.2 Untersuchung der Referenzgeometrie 81 


0,60 


P4 P3 P2 Pi SP S1 S2 S3 S4 S586 S7 S8 S9 
Position, param. Lauflänge [-] 


Abbildung 6.4: Verlauf der gemittelten Totalkühleffektivität der Referenzgeometrie, drei Be- 
triebspunkte 


Aufgrund der hohen Anstellwinkel und damit großen inneren Oberflächen der Bohrungen kann 
die Kühlluft schon vor dem Austritt aus dem Diffusor Wärme aufnehmen. 

Am Punkt S3 direkt nach den Bohrungen kann von der TCE direkt auf die Filmkühleffektivität 
geschlossen werden, da diese an dieser Position der maßgebliche Effekt der Kühlung ist. Hier 
wird ein vom Betriebspunkt unabhängiger Maximalwert erreicht. Eine Erhöhung des Kühl- 
luftmassenstromes (und damit des Druckverhältnisses IT bzw. des Innendrucks px jot) führt zu 
keiner weiteren Erhöhung der Filmkühleffektivität. Dies ist auf ein zu hohes Impulsverhältnis 
und damit ein Ablösen der Kühlluftstrahlen, vor allem nahe des Schaufelfußes, zurückzuführen. 
Dieser Effekt konnte auch von Stichling (2018) numerisch gezeigt werden. 

Im stromab liegenden Schaufelbereich nimmt die Totalkühleffektivität graduell ab. Zwar dickt 
die Grenzschicht auf und der heißgasseitige Wärmeübergangskoeffizient hx nimmt ab, dies wird 
aber von einer abnehmenden Kühlwirkung überlagert: Wie in Abbildung 6.2(a) gezeigt, wan- 
dern die saugseitigen Äste der naben- und gehäuseseitigen Hufeisenwirbel an der Schaufelober- 
fläche nach innen und mischen dabei Heißgas ein. Dies führt im Mittel über die Schaufelhöhe 
zu einer reduzierten Kühleffektivität. Weiterhin nimmt die Filmkühleffektivität im mittleren 
Schaufelbereich durch Ausmischung in Strömungsrichtung ab. Letztlich ist die interne Küh- 
lung der Hinterkante aufgrund der kleinen Kanäle und damit ungünstigen Flächenverhältnisse 
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zwischen Außenfläche und Kühlkanalfläche? und der im hinteren Bereich des Kühlsystems 
deutlich ansteigenden Temperaturen der Kühlluft ineffektiv. Die Kühleffektivität sinkt an der 
Hinterkante nahezu auf den Minimalwert der Vorderkante. 


6.3 Untersuchung der Drallgeometrien 


Die Untersuchung der Kühleffektivität der Drallgeometrien wird analog zur Untersuchung der 
konventionellen Geometrie durchgeführt. Neben den absoluten Daten werden zusätzlich relati- 
ve Effektivitäten untersucht, als Referenz dient der konventionelle Fall am jeweiligen Betriebs- 
punkt. In den folgenden Darstellungen sind links (durchgezogene Linien) die absoluten Wer- 
te der Totalkühleffektivität analog zu Abbildung 6.4 dargestellt, rechts daneben (gestrichelte 
Linien) werden die relativen Daten gezeigt. Da für die konventionelle Geometrie bereits Ma- 
schinendaten vorhanden sind, ergibt sich so eine zweite Möglichkeit der Skalierung zurück zur 
Maschine. Diese relative Skalierung ist nochmals robuster, da so auch Auslegungsfehler des 
Prüfstandes ausgeglichen werden. 


6.3.1 Künhleffektivität der GVI-Geometrie 


Die (GVD-Geometrie ist als axialer Drallerzeuger mit Nabe und mehreren Schaufeln mit kurz- 
en Sehnen aufgebaut (siehe auch Abbildung 3.3). Bei gegebenen Massenströmen benötigt die- 
ser Drallerzeuger das geringste Eintrittstotaldruck verhältnis II aller Konfigurationen inklusive 
der konventionellen Geometrie. Dies führt zu einer niedrigen Drallzahl und damit zu einem 
vergleichsweise niedrigen internen Wärmeübergang. Dies muss speziell im Bereich der Vor- 
derkante sichtbar sein, da die Drallgeometrien hier auf zusätzliche Turbulatoren an den Innen- 
wänden verzichten. Die absolute und die relative Performance der (GVI)-Geometrie wird in 
Abbildung 6.5 gezeigt. 

Der qualitative Verlauf der absoluten Kühleffektivität, gezeigt in Abbildung 6.5(a), ist mit dem 
der konventionellen Geometrie vergleichbar. Dieser ergibt sich maßgeblich aus den in Ab- 
schnitt 6.1 beschriebenen, vom Heißgas aufgeprägten Strukturen. Die absolute Totalkühleffek- 
tivität am Staupunkt ist mit TCE ~ 32% deutlich geringer und unter dem für einen sicheren 
Betrieb notwendigen Wert. Die Gradienten stromab des Staupunktes, insbesondere entlang der 
Saugseite, sind größer, was zu einer unerwünschten Erhöhung der thermischen Spannungen 
führen Kann. Die Betriebspunktabhängigkeit der Totalkühleffektivität ist wieder annähernd line- 
ar. Das bei der konventionellen Geometrie zu beobachtende Zusammenfallen der Kurven nach 
den saugseitigen fanshaped-Bohrungen bei S3 entfällt. Die Rotationsrichtung des Dralls ist ge- 
genläufig zur Richtung der Filmkühlbohrungen. Dies führt nach Unterabschnitt 2.3.1 zu einer 
Reduktion des Durchflusskoeffizienten und damit auch des Impulsverhältnisses. Durch diese 
Reduktion hebt der Kühlluftstrahl heißgasseitig nicht von der Oberfläche ab, ein geschlossener 


2 Die Hinterkante der verwendeten Rotorschaufelgeometrie wird durch stark angestellte, diskrete Bohrungen 
auf der Druckseite gekühlt. Der momentane Stand der Technik sieht weitaus komplexere Kühlsysteme für die 
Hinterkante vor. Hier werden inzwischen großflächige Hinterkantenschlitze, eine Entfernung der druckseitigen 
Fläche nahe der Hinterkante und eine Vielzahl interner Turbulatoren eingesetzt, die alle zu einer Steigerung 
der Kühleffektivität im Hinterkantenbereich führen. 
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TCE rel. [-] 


Position, param. Lauflange [-] Position, param. Lauflange [-] 


(a) Totalkühleffektivität (GVI) absolut (b) Totalkühleffektivität (GVD relativ 


Abbildung 6.5: Verlauf der gemittelten Totalkühleffektivität, (GVI), drei Betriebspunkte 


Kühlfilm wird gebildet. Eine Veränderung zu Betriebspunkten mit erhöhtem Massenstrom führt 
dann zu einer Erhöhung des internen Totaldruckes px to, und damit zu einer Erhöhung der Aus- 
blaserate bei geschlossenem Kühlfilm, was die Filmkühleffektivität erhöht. 

Der relative Verlauf der Kühleffektivität wird in Abbildung 6.5(b) gezeigt. Auch im Vergleich 
zeigt sich die stark reduzierte Kühlwirkung im Bereich der Vorderkante. Das relative Verhalten 
ist, bis auf die Abweichung an S3, im Rahmen der Messgenauigkeit betriebspunktunabhängig. 
Im Mittel über die gesamte Schaufel stellt sich eine verringerte Kühleffektivität ein. Nur im 
Bereich der saugseitigen Filmkühlbohrungen, vor allem unmittelbar stromab dieser, kann eine 
Steigerung der Totalkühleffektivität erreicht werden. Da hier die Filmkühleffektivität den größ- 
ten Einfluss auf die Totalkühleffektivität hat, lässt sich dieser Sachverhalt durch zwei Effekte 
erklären: Durch die oben beschriebenen Veränderungen in der Kinematik der Durchströmung 
dieser Bohrungen wird die Filmkühleffektivität erhöht. Weiterhin ist die Kühlluft durch die 
verringerte Wärmeaufnahme im Vorderkantenbereich im Mittel kälter, was insbesondere im 
Bereich starker Filmkühlung zu einer Reduktion des treibenden Temperaturgradienten auf der 
Heißgasseite führt. 

Aufgrund der im Mittel stark gesunkenen Totalkühleffektivität TCE, der zu niedrigen Total- 
kühleffektivität im Vorderkantenbereich und der hohen Gradienten kann dieser Drallerzeuger 
nicht für die Verwendung empfohlen werden. Des Weiteren ist die Fertigung mit konventio- 
nellen Verfahren (z.B. Feinguss) nicht sicherzustellen. Daher müssen weitere Drallerzeuger 
untersucht werden. 


6.3.2 Künhleffektivität der NSG-Geometrie 


Die (NSG)-Geometrie ist weiterhin als axialer Drallerzeuger mit zentraler Nabe aufgebaut. Die 
einzelnen Leitschaufeln mit kurzer Sehne der (GVI)-Geometrie werden durch drei stark ver- 
wundene Leitschaufeln mit hoher Sehnenlänge ersetzt (siehe auch Abbildung 3.3). Das Ein- 
trittstotaldruckverhältnis II steigt und liegt für die untersuchten Betriebspunkte über dem der 
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Abbildung 6.6: Verlauf der gemittelten Totalkühleffektivität, (NSG), drei Betriebspunkte 


konventionellen Geometrie. Die gemittelten Ergebnisse der experimentellen Untersuchung wer- 
den in Abbildung 6.6 gezeigt. Die absoluten Daten, Abbildung 6.6(a), zeigen qualitativ einen 
der (GVI)-Geometrie vergleichbaren Verlauf. Das mittlere Niveau der Totalkühleffektivität so- 
wie die Totalkühleffektivität am Staupunkt sind höher im Vergleich zur (GVI)-Geometrie. Dies 
ist auf die mutmaßlich (numerische Voruntersuchungen von Rolls-Royce Deutschland, 2017) 
erhöhte Drallzahl Y, abgeleitet vom erhöhten Eintrittstotaldruckverhältnis II, zurückzuführen. 
Stromab der Vorderkante in Richtung der Saugseite sind wieder starke Gradienten zu beob- 
achten. Die maximale Totalkühleffektivität wird an S3 erreicht. Die Betriebspunktabhängig- 
keit nimmt ab. Mit steigender Drallzahl steigt die tangentiale Geschwindigkeitskomponente 
und damit der relative Einströmwinkel in die Bohrungen und begrenzt die maximale Ausbla- 
serate schon bei Betriebspunkten mit niedrigem Massenstrom. Der Vergleich zur konventionel- 
len Geometrie, dargestellt in Abbildung 6.6(b), zeigt weiterhin eine schlechtere Performance. 
Die Kühleffektivität am Staupunkt und unmittelbar stromab desselben sowohl in Richtung der 
Druckseite als auch der Saugseite ist geringer als die der konventionellen Geometrie. Die interne 
Kühlung und die Versorgung der zylindrischen Filmkühlbohrungen am Staupunkt kann durch 
den internen Drall nicht gewährleistet werden. Die Druckseite stromab zeigt eine mit dem kon- 
ventionellen Fall vergleichbare Totalkühleffektivität. Die Saugseite zeigt speziell im Bereich der 
Diffusorbohrungen eine breiter verteilte, leicht verbesserte Kühleffektivität. Dies kann wieder 
durch die kältere Kühlluft erklärt werden. Im Bereich der Hinterkante fällt die Kühleffektivität 
dann stärker ab. Da die interne Geometrie an dieser Position bei allen Prüflingen gleich ist, kann 
eine mögliche Erklärung das im Mittel höhere Temperaturniveau der Schaufel und die damit ge- 
ringere Wärmeleitung von der dünnen Hinterkante zu stärker gekühlten Bereichen sein. Durch 
eine Steigerung der Kühleffektivität am Staupunkt konnten die Gradienten reduziert werden. 
Im Mittel liegt die Totalkühleffektivität immer noch unter der der konventionellen Geometrie. 
Speziell an der Vorderkante und direkt stromab dieser ist die Totalkühleffektivität geringer. Wei- 
terhin ist die Fertigung der (NSG)-Geometrie mit Klassischen Gussverfahren schwer möglich. 
Der Einsatz dieser Geometrie kann also für eine Verwendung in der Maschine ebenfalls nicht 
empfohlen werden. 


6.3 Untersuchung der Drallgeometrien 85 


TCE rel. [-] 


P4 SP S9 
Position, param. Lauflänge [-] Position, param. Lauflänge [-] 
(a) Totalkühleffektivität (OSG0) absolut (b) Totalkühleffektivität (OSG0) relativ 


Abbildung 6.7: Verlauf der gemittelten Totalkühleffektivität, (OSG0), drei Betriebspunkte 


6.3.3 Künhleffektivität der OSG-Geometrien 


Die zuvor präsentierten Ergebnisse der beiden Drallgeometrien (GVD und (NSG) konnten 
noch kein positives Urteil im Sinne einer technischen Nutzung dieser neuen Konzepte zulas- 
sen. Aus diesem Grund wurde eine dritte geometrische Variation von Drallerzeugern gefertigt. 
Die (OSGi)-Geometrien (siehe Abbildung 3.3) werden ohne zentrale Nabe gefertigt. Zwei Leit- 
gitter werden am äußeren Umfang angebracht und erzeugen so die Drallströmung. Im Bezug 
auf die vorher aufgezeigten Mängel ergeben sich so neue Potentiale: Durch den Aufbau ohne 
Nabe ist die (Serien-) Fertigung im Feinguss- Verfahren der eigentlichen Rotorschaufel prinzipi- 
ell möglich. Durch eine parametrische Auslegung lässt sich die Versperrung und damit sowohl 
das Eintrittstotaldruck verhältnis II als auch die Drallzahl Y variieren (numerische Voruntersu- 
chungen durchgeführt von Rolls-Royce Deutschland, 2017). 

Im Folgenden werden die Ergebnisse der experimentellen Untersuchungen von drei Variatio- 
nen der (OSG)-Geometrie vorgestellt: (OSG0), (OSG3) und (OSG8). Hierbei beschreibt der 
Laufindex im Namen die geometrische Versperrung (steigend). Aus Tabelle 6.1 können die 
Randbedingungen der Versuche entnommen werden. Wie aus den Ergebnissen der Vorausle- 
gung zu erwarten war, steigt das Eintrittstotaldruck verhältnis II mit steigendem Laufindex an, 
die Drallzahl Y nimmt mutmaßlich zu. Während (OSG0) und (OSG3) vergleichbare Druckver- 
hältnisse wie die (NSG)-Geometrie zeigen, ist die (OSG8)-Geometrie auf ein deutlich höheres 
Eintrittstotaldruckverhältnis ausgelegt. 


Die Ergebnisse der Untersuchung der ersten Geometrie (OSG0) werden in Abbildung 6.7 darge- 
stellt. Die absoluten Daten, Abbildung 6.7(a), zeigen wieder den erwarteten qualitativen Verlauf 
der vorher untersuchten Geometrien. Am Staupunkt konnte das mittlere Niveau der Kühleffek- 
tivität trotz niedrigerem Eintrittstotaldruck verhältnis gegenüber der (NSG)-Geometrie gehalten 
werden. Auch analog zu dieser Geometrie stellt sich der Verlauf der Totalkühleffektivität ent- 
lang der Saugseite ein. Das Maximum der Kühleffektivität verschiebt sich leicht stromab in 
Richtung S4, der folgende Gradient in Strömungsrichtung ist geringer. Dies führt zu einer wei- 
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Abbildung 6.8: Verlauf der gemittelten Totalkühleffektivität, (OSG3), drei Betriebspunkte 


teren Überdeckung des gut gekühlten Bereichs und damit einer im Mittel sinkenden Prüflings- 
temperatur. Dies ist vor allem im Bereich der Hinterkante im Vergleich zu dem konventionellen 
Konzept (Abbildung 6.7(b)) zu beobachten. Hier wird erstmals durchgehend eine erhöhte Kühl- 
effektivität beobachtet. Weiterhin lässt sich anhand der relativen Kühleffektivitäten die immer 
noch präsente, reduzierte Kühlwirkung des Staupunktes und des Bereiches unmittelbar stromab 
des Staupunktes auf Druck- sowie Saugseite erkennen. Auf der Druckseite ergibt sich eine bis- 
her nicht beobachtete Anomalie, bei der der niedrigste Betriebspunkt eine im Vergleich deut- 
lich erhöhte Kühleffektivität aufzeigt. Dies ist auch bei den folgenden (OSG)-Geometrien in 
unterschiedlicher Ausprägung zu erkennen. Es wird vermutet, dass komplexe, betriebspunkt- 
abhängige Sekundärströmungen im ersten radial nach innen durchströmten Kanal (siehe auch 
Abbildung 3.2(b)) die Einströmung in die Filmkühlbohrungen beeinflussen. 

Zusammenfassend erreicht die (OSG0)-Geometrie eine der (NSG)-Geometrie vergleichbare 
Kühlwirkung bei geringerem Eintrittstotaldruckverhältnis II. Gleichzeitig ist diese Gestaltung 
des Drallerzeugers in den Fertigungsprozess integrierbar. 


Im Folgenden werden die in Abbildung 6.8 dargestellten Verläufe der gemittelten Totalkühl- 
effektivitäten der (OSG3)-Geometrie untersucht. Diese Geometrie zeigt eine erhöhte Versper- 
rung und damit ein erhöhtes Eintrittstotaldruckverhältnis II. Schon in den in Abbildung 6.8(a) 
dargestellten Verläufen der absoluten Totalkühleffektivität lässt sich eine deutliche Steigerung 
derselben erkennen. Im Bereich der fanshaped-Bohrungen auf der Saugseite an S3 werden mit 
TCE ~ 65% die bisher höchsten Effektivitäten beobachtet. Die Abhängigkeit vom Betriebs- 
punkt, gerade nach den Filmkühlbohrungen, ist weiterhin niedrig. Der Gradient in Strömungs- 
richtung nach der Position S3 ist ebenso niedrig, eine breite Überdeckung der gut gekühlten 
Regionen kann erreicht werden. Am Staupunkt wird mit TCE ~ 40% erstmals eine der kon- 
ventionellen Geometrie ähnliche Kühlwirkung erreicht. 

Die relativen Daten, dargestellt in Abbildung 6.8(b), können diese positive Prognose bestätigen. 
Bis auf den kritischen Bereich am Staupunkt und direkt stromab auf der Saugseite Konnte die 
Totalkühleffektivität erhöht werden. Auf der Druckseite kann eine Steigerung von = 10% er- 
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Abbildung 6.9: Verlauf der gemittelten Totalkühleffektivität, (OSG8), drei Betriebspunkte 


reicht werden, die auf die verbesserte Wärmeabfuhr an der Kanalinnenseite zurückzuführen ist. 
Auch auf der Saugseite konnte die Totalkühleffektivität erhöht werden. Diese Erhöhung liegt 
im besten Fall bei bis zu ~ 22 %. Hierzu trägt im vorderen Bereich der Saugseite die verbesser- 
te Innenkühlung bei, an den Diffusorbohrungen (S2 und S3) dominiert mutmaßlich das besser 
angepasste Impulsverhältnis. Der stromab liegende Bereich bis zur Hinterkante profitiert von 
einem generell gesunkenem Temperaturniveau des Prüflings. Durch Verwendung der (OSG3)- 
Geometrie kann eine global verbesserte Kühlung erreicht werden. Der Bereich am Staupunkt 
und der direkt stromab liegenden Saugseite verbleibt mit leicht verschlechterter Kühlwirkung. 


Zuletzt werden die Kühleffektivitäten der (OSG8)-Geometrie diskutiert, die in Abbildung 6.9 
dargestellt werden. Bei dieser Geometrie liegt das Eintrittstotaldruckverhältnis mit II = 1,26 
am Nennbetriebspunkt über dem der anderen untersuchten Drallerzeuger. Trotzdem kann keine 
weiter verbesserte Kühlwirkung erreicht werden. Am Staupunkt ist die erreichte Totalkühlef- 
fektivität ähnlich der der (OSG3)-Geometrie. Der Verlauf stromab entlang der Saugseite ist 
ebenfalls vergleichbar, die erreichte maximale Totalkühleffektivität, speziell nach S3, bleibt un- 
ter der der (OSG3)-Geometrie. Die weiter erhöhte Drallzahl führt mutmaßlich zu stark gestörten 
Einströmbedingungen für die Diffusorbohrungen und damit zu einer gesenkten Filmkühleffek- 
tivität. 

Auch die relativen Daten, gezeigt in Abbildung 6.9(b), können diesen Sachverhalt abbilden. 
Zwar wird wieder eine, bis auf den Bereich am Staupunkt, durchgehend erhöhte Totalkühlef- 
fektivität erreicht, das mittlere Niveau sowie das Maximum auf der Saugseite ist aber geringer 
als beim vorher untersuchten Fall. Die bei der (OSG0)-Geometrie beschriebene Anomalie auf 
der Druckseite ist wieder zu beobachten. So wird trotz einer weiteren Erhöhung des Eintrittsto- 
taldruckverhältnisses II keine maßgebliche Verbesserung der Totalkühleffektivität erreicht. 
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Abbildung 6.10: Verlauf der gemittelten Totalkühleffektivität, Vergleich der Geometrien bei 
Nennbetriebspunkt 


6.4 Vergleich und abschließende Bewertung der Konzepte 


Abschließend werden im Folgenden alle untersuchten Geometrien nochmals verglichen. Es 
werden wieder die absoluten und relativen, über die Schaufelhöhe gemittelten, Totalkühlef- 
fektivitäten, dargestellt in Abbildung 6.10, verwendet. Aufgrund der geringen, meist annähernd 
linearen Abhängigkeit des Betriebspunktes werden für jedes Konzept nur die Daten am Nenn- 
betriebspunkt Cjo9 gezeigt. Die absoluten Daten sind in Abbildung 6.10(a) zusammengefasst, 
die relativen Daten in Abbildung 6.10(b). Beide Abbildungen sind außerdem in vergrößerter 
Darstellung im Anhang in Abschnitt A.4 und Abschnitt A.5 zu finden. Die gezeigte Spanne von 
Kühleffektivitäten von TCE = 30% bis TCE ~ 70% können so auch in der Literatur, Unterab- 
schnitt 2.3.2, gefunden werden und sind in guter Übereinstimmung mit den erwarteten Werten 
aus der Betriebserfahrung von Rolls-Royce Deutschland. Der Prüfstand in Kombination mit 
dem Messverfahren kann daher zur quantitativen Bestimmung solcher Totalkühleffektivitäten 
verwendet werden. 


Am Staupunkt und direkt stromab in Richtung der Saugseite kann keines der neuen Konzep- 
te die Totalkühleffektivität der konventionellen Geometrie erreichen. Stichling (2018) konn- 
te anhand von auf die experimentellen Ergebnisse angepassten CFD-Berechnungen numerisch 
zeigen, dass der durch die zwei Reihen der zylindrischen Kühlluftbohrungen ausgeblasene Mas- 
senstrom global trotz veränderter Innengeometrie nur geringfügig variiert. Aufgrund stark ge- 
störter Einströmbedingungen bei Verwendung der Drallkühlung führt ein vergleichbarer Ge- 
samtmassenstrom zu hohen Abweichungen zwischen den Bohrungen und damit zu hohen lo- 
kalen Geschwindigkeiten in einzelnen Bohrungen. Dies führt bei Verwendung der kurzen, zy- 
lindrischen Bohrungen zu Ablöseeffekten und einer reduzierten Filmkühleffektivität stromab 
der heißgasseitigen Öffnung. Die sich dadurch einstellende Erhöhung des externen Wärme- 
stroms kann nicht von der verbesserten Innenkühlung ausgeglichen werden. Des Weiteren ist 
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Abbildung 6.11: Verlauf der gemittelten Totalkühleffektivität über die normierte Schaufelhöhe 
an drei Positionen nahe der Vorderkante 


die Wandstärke am Staupunkt durch die Einbringung eines einzelnen Vorderkantenkanals mit 
möglichst kreisförmigem Querschnitt höher. Dies bedingt einen erhöhten thermischen Wider- 
stand, der nicht durch den erhöhten inneren Wärmeübergang ausgeglichen werden kann. Sobald 
die Wandstärken stromab der Vorderkante vergleichbar sind, steigen die von den neuen Kon- 
zepten erreichten Totalkühleffektivitäten. 

Stromab des Staupunktes entlang der Druckseite (Pi) verschiebt sich das mittlere Niveau der To- 
talkühleffektivität mit der mittleren Schaufeltemperatur, der Verlauf bleibt qualitativ vergleich- 
bar. Da die druckseitigen fanshaped-Bohrungen aus dem zweiten, radial nach innen laufenden 
Kanal gespeist werden, ist hier keine maßgebliche Veränderung bei Variation der Geometrie des 
Vorderkantenkanals zu erwarten. 

Auch stromab des Staupunktes entlang der Saugseite (Si) lässt sich eine Parallelverschiebung 
der Totalkühleffektivitäten feststellen. Diese ist im vorderen Schaufelbereich auf die verbes- 
serte Innenkühlung mit steigendem Eintrittstotaldruck verhältnis IT und daraus mutmaßlich re- 
sultierender erhöhten Drallzahl Y zurückzuführen. Deutliche Unterschiede ergeben sich kurz 
stromauf sowie stromab der saugseitigen Diffusorbohrungen durch die komplexe Wechselwir- 
kung zwischen interner Strömung und Durchströmung der Filmkühlbohrungen. Hier kann für 
den Drallerzeuger (OSG3) ein Optimum festgestellt werden. 


Abschließend soll der radiale Verlauf der Totalkühleffektivität (über die Schaufelhöhe) beurteilt 
werden. Die relative Totalkühleffektivität TCE (also das lokale Verhältnis der Totalkühleffek- 
tivität der Drallgeometrie zur konventionellen Referenzgeometrie) ist für zwei verschiedene 
Drallgeometrien in Abbildung 6.11 über die Schaufelhöhe aufgetragen dargestellt. Lokal unter- 
schiedliche Heißgaswärmeübergänge, die für die drei geometrischen Variationen vergleichbar 
sind, können so vernachlässigt werden. Es wurden vier Positionen gewählt: Der Staupunkt (SP), 
die Position unmittelbar stromab in Richtung der Druckseite (Pl) und zwei Positionen unmit- 
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telbar stromab in Richtung der Saugseite (S1, S2). Das Verhalten ist für beide Drallkonzepte 
vergleichbar, lediglich das mittlere Niveau der Kühleffektivität ist bei der in Abbildung 6.11(b) 
dargestellten (OSG3)-Geometrie höher als bei der in Abbildung 6.11(a) dargestellten (NSG)- 
Geometrie. Dies deckt sich mit den im vorherigen Abschnitt diskutierten Ergebnissen. Bei der 
Verteilung der relativen Kühleffektivität lässt sich, bis auf die Position S1 unmittelbar nach 
den Filmkühlbohrungen, eine Reduktion mit zunehmender Schaufelhöhe feststellen. So erreicht 
die (OSG3)-Geometrie am Staupunkt am Schaufelfuß eine der konventionellen Geometrie ver- 
gleichbare Kühleffektivität. Mit steigender Schaufelhöhe fällt im Innenkanal die Drallzahl auf- 
grund von Dissipation, der interne Wärmeübergangskoeffizient sinkt. Dies führt zu einer sin- 
kenden Kühleffektivität mit steigender Schaufelhöhe. Beim konventionellen Fall wird durch die 
über die gesamte Schaufelhöhe eingebrachten Rippen erneut Turbulenz generiert, der Wärme- 
übergangskoeffizient an der Innenseite wird kontinuierlich hoch gehalten. Weiter stromab vom 
Staupunkt (P1, S2) verbleibt dieser Effekt. Aufgrund sinkender Wandstärken ist die relative 
Erhöhung der Kühleffektivität über die Schaufelhöhe durch Verbesserung des internen Wärme- 
übergangskoeffizienten ausgeprägter als am Staupunkt. 

Die Reduktion des internen Wärmeübergangskoeffizienten in Richtung der Schaufelhöhe kann 
durch die Korrelation von Chang et al. (1994) mit der Drallzahl verknüpft werden. Hierfür wird 
die Korrelation, die aus Ergebnissen für Drallerzeugung durch tangentiale Einblasung abgeleitet 
wurde, anhand eines Korrekturfaktors an die erwarteten Drallzahlen angepasst. Zur Darstellung 
werden die berechneten Verläufe entlang der y-Achse verschoben. Somit kann nur der Verlauf, 
nicht aber der absolute Wert interpretiert werden. Diese Verläufe sind in guter Übereinstimmung 
mit den Verläufen der relativen Totalkühleffektivität an den Position Pl und S2 mit geringer 
Wandstärke. 

An der Position S1, die auch im Vergleich des vorherigen Abschnittes bei jeder Drallkonfigura- 
tion die minimale Totalkühleffektivität gezeigt hat, ist der Trend genau umgekehrt. Hier ergibt 
sich die starke Abhängigkeit der Filmkühleffektivität von den internen Anströmbedingungen 
und dem externen compound-Winkel. Erst bei weitestgehend abgebautem Drall im Innenka- 
nal (steigende Schaufelhöhe) steigt die Kühleffektivität aufgrund steigender Filmkühleffekti- 
vität wieder an. Auf der Druckseite ist dieser Effekt nicht zu beobachten. Diese Unterschiede 
sind vermutlich auf die Drallrichtung der internen Strömung (abgewandt von den Einström- 
querschnitten der saugseitigen Bohrungen, zugewandt für die druckseitigen Bohrungen) sowie 
die unterschiedlichen Heißgasbedingungen zurückzuführen (Beschleunigung auf der Saugseite, 
Verzögerung entlang der Druckseite). 


Die in diesem Kapitel präsentierten Ergebnisse der experimentellen Untersuchungen erlauben 
erstmals die Beurteilung neuer, auf Drallkühlung basierender Kühlkonzepte unter Beachtung 
aller Randbedingungen und der Einhaltung der dimensionslosen Kennzahlen. Die in vorange- 
gangenen Untersuchungen als positiv bewerteten Konzepte konnten so auf ihre Eignung bei 
Implementation in realistische Schaufelgeometrien überprüft werden. 

Es konnte gezeigt werden, dass erst ab einem erhöhten Eintrittstotaldruckverhältnis von 
II ~ 1.2 und mit geometrisch komplexen axialen Leitgittern Verbesserungen erzielt werden 
konnten. So konnte die Kühleffektivität der Prüflinge um bis zu ~ 20% erhöht werden. 
Große Bereiche der Druckseite sowie der Saugseite zeigen bei dem besten Drallkühlkonzept 
(OSG3) deutlich erhöhte Totalkühleffektivitäten. Im thermisch kritischen Schaufelbereich 
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um die Vorderkante konnte keines der untersuchten Drallkonzepte eine im Vergleich zu dem 
konventionellen Konzept verbesserte Kühlwirkung erzielen. Die Gründe hierfür sind die 
notwendigen geometrischen Änderungen des Kühlluftkanals und die starke Beeinflussung der 
Filmkühlbohrungen, die aus diesem gespeist werden. Hier kann bei zylindrischen Bohrungen 
eine sinkende Filmkühleffektivität beobachtet werden, bei Diffusorbohrungen kann im vorlie- 
genden Fall eine Steigerung beobachtet werden. 

Die Verwendung von Drallkühlung als direkter Austausch eines bewährten, konventionellen 
Kühlkonzeptes kann anhand der vorliegenden Untersuchungen noch nicht uneingeschränkt 
empfohlen werden. Eine genaue Abstimmung des System ist notwendig. So muss zuerst 
die Geometrie der internen Kühlkanäle im Hinblick auf homogene Wandstärken angepasst 
werden. Nach Auswahl eines geeigneten Drallerzeugers muss dann die Filmkühlung genau 
auf die interne Strömung optimiert werden. Hier ergibt sich vor allem an der Vorderkante 
bei Verwendung zylindrischer Bohrungen ein erhöhter Kühlbedarf, der z.B. durch größere 
Bohrungsdurchmesser gewährleistet werden kann. Weiterhin muss die Reduktion der Drallzahl 
und damit des Wärmeübergangskoeffizienten berücksichtigt werden. Bei konventionellen 
Konzepten wird diese Größe durch Rippen oder ähnliche Strukturen erneut erhöht. Diese kon- 
tinuierliche Erneuerung entfällt bei den hier untersuchten Konzepten und führt zu reduzierten 
Künhleffektivitäten im Schaufelspitzenbereich. 


Aus den hoch aufgelösten, dreidimensionalen Temperaturfeldern lassen sich, auch unabhän- 
gig der verwendeten internen Kühlungsgeometrien, weitere Erkenntnisse ableiten. So muss zur 
Vermeidung von lokalen Temperaturgradienten das Sekundärströmungsfeld des Heißgases kon- 
trolliert und optimiert werden. Hierbei ist der sich an der Vorderkante bildende, nabenseitige 
Hufeisenwirbel (dessen Intensität und Verlauf in die Kanalmitte) maßgeblich. Dies wird nach 
aktuellem Stand der Technik mit konturierten Seitenwänden erreicht. Sollte dies nicht ausrei- 
chen, ist eine lokal angepasste Filmkühlung in Betracht zu ziehen. Außerdem muss die Pro- 
gnose der Kühleffektivität von stark angestellten zylindrischen Bohrungen, wie sie im Bereich 
der Vorderkante eingesetzt werden, in Zukunft weiter untersucht werden. Zwar existieren sol- 
che Untersuchungen in Bezug auf die Veränderung der Anströmwinkel auf der Kühlluftseite, 
Untersuchungen zur heißgasseitigen Anströmung unter hohen compound-Winkeln bei maschi- 
nenähnlichen Randbedingungen existieren aber bisher nicht. 


7 Zusammenfassung und Ausblick 


Die Kühlung von Komponenten in Gasturbinen, speziell den rotierenden Laufschaufeln der 
ersten Turbinenstufe von Flugzeugtriebwerken, verbleibt weiterhin als große technische Her- 
ausforderung. Neue Entwicklungstrends führen zu weiter steigenden Turbineneintrittstempe- 
raturen, steigenden Druckverhältnissen und damit steigenden Kühllufttemperaturen. Die Ein- 
satztemperaturen verfügbarer Materialien verbessern sich dagegen nur langsam. Eine effektive 
Kühlung der von Heißgas umströmten Bauteile kann höhere Wirkungsgrade für neue Maschi- 
nen ermöglichen. Während eine breite Datenbasis zum Verhalten der einzelnen Teilkomponen- 
ten eines Schaufelkühlsystems existiert, gibt es nur wenige Untersuchungen, die die kombinier- 
te Kühlwirkung von Innen- und Filmkühlung beschreiben. Keine dieser Untersuchungen konnte 
bisher eine volle aerothermische Ähnlichkeit und gleichzeitig eine hoch aufgelöste Messung der 
Zielgrößen sicherstellen. 

Das Ziel der vorliegenden Arbeit war daher die Untersuchung von innovativen Kühlkonzep- 
ten für Turbinenrotorschaufeln hinsichtlich ihres aerodynamischen Verhaltens und ihrer thermi- 
schen Leistungsfähigkeit. Diese wurde durch die Bestimmung der Totalkühleffektivität, also der 
Temperatur des Schaufelmaterials relativ zu den Freistromtemperaturen von Heißgas und Kühl- 
luft, dargestellt. Damit die im Versuch bestimmten Zielgrößen für die Auslegung von Bauteilen 
herangezogen werden können, muss eine weitestgehende aerothermische Ähnlichkeit eingehal- 
ten und trotzdem ein geringer Restmessfehler garantiert werden. 

Um diese Ähnlichkeit zu erreichen, wurde ein neuer Prüfstand konzipiert und aufgebaut. In 
umfangreichen Voruntersuchungen Konnte gezeigt werden, dass alle für die Aerodynamik und 
den konjugierten Wärmeübergang wichtigen dimensionslosen Kennzahlen eingehalten werden 
können. Lediglich auf den Effekt der Rotation muss verzichtet werden. So können unter Ver- 
wendung des neuen Prüfstandes komplette Kühlkonzepte hinsichtlich ihrer Effektivität beurteilt 
und die Ergebnisse für die finale Bauteildimensionierung verwendet werden. 

Diese Analyse unter voller Ähnlichkeit stellt hohe Anforderungen an die in dieser Arbeit als 
Messtechnik eingesetzte Infrarotthermographie. Die Oberflächentemperatur muss großflächig, 
aber trotzdem hoch aufgelöst bestimmt werden. Hier müssen mehrere Größenordnungen an 
Längenmaßen berücksichtigt werden. Gleichzeitig ergibt sich bei der Untersuchung der stark 
gekühlten Bauteile in geschlossenen Kanälen eine hohe Strahlungsdynamik sowie ein unvorteil- 
haftes Verhältnis zwischen Umgebungs- und Objektstrahlung. Diese hochkomplexen Messun- 
gen werden nun erstmals zusätzlich an uneingeschränkt dreidimensionalen Geometrien durch- 
geführt. Diese Dreidimensionalität führt zu lokalen Veränderungen der radiativen Eigenschaften 
der Oberfläche abhängig vom Betrachtungswinkel. Des Weiteren wird die Zuordnung der Da- 
ten zwischen Bild und Objekt erschwert. Trotz dieser extremen Anforderungen sind die expe- 
rimentellen Untersuchungen nur sinnvoll, wenn der verbleibende Restfehler auf ein Minimum 
reduziert wird: Die finalen Bauteile reagieren hoch sensitiv auf eine Veränderung der Einsatz- 
temperatur, der Restfehler muss daher bekannt und mit ca. 2,5 % niedrig sein. 

Um dies zu gewährleisten wurden neue Verfahren entwickelt, die eine Bestimmung der Ober- 
flächentemperatur mittels Infrarotthermographie auf beliebig komplexen Oberflächen auch bei 
unvorteilhaften Strahlungssituationen ermöglichen. Die Grundlage dieser Verfahren ist die ge- 
naue Kenntnis der Kameraposition, genauer dem optischen Zentrum. Da dieses messtechnisch 
schwierig zu bestimmen ist, wurde ein neuer genetischer Algorithmus präsentiert, der die Be- 
stimmung des Kamerazentrums erlaubt. 
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Die erreichte Genauigkeit bezogen auf die resultierenden Sichtwinkel liegt bei < 1,5°. Ausge- 
hend von dem so bestimmten optischem Zentrum können die neuen Kalibrierverfahren ange- 
wendet werden. Diese korrigieren den Einfluss der Veränderung der radiativen Eigenschaften 
der Oberfläche mit dem Betrachtungswinkel. Zusätzlich wurde gezeigt, dass durch eine Diskre- 
tisierung der Messoberfläche eine Übertragung der Ergebnisse der thermischen Messungen aus 
dem Bildraum in den Objektraum ermöglicht wird. Die Datenfelder werden dadurch entzerrt 
und der physikalischen Geometrie zugeordnet. In einer finalen Fehlerabschätzung konnte ab- 
schließend gezeigt werden, dass der Fehler in einem weitem Bereich bis zu einem Sichtwinkel 
von 80° unter dem geforderten minimalen Restfehler liegt. 

Anhand der so kalibrierten Messdaten konnten darauffolgend verschiedene neue Konzepte zur 
Erhöhung des internen Wärmeübergangs von Turbinenrotorschaufeln untersucht werden. Hier- 
bei wurde nur die interne Geometrie variiert, die Außengeometrie und die Geometrie der Film- 
kühlung wurden nicht verändert. Der Fokus der Untersuchungen lag auf der Analyse der To- 
talkühleffektivität von Drallkühlkonzepten im Vorderkantenkanal der Schaufel. Verschiedene 
axiale Drallerzeuger wurden im Schaufelfuß eingesetzt. Zur besseren Einordnung der Daten 
und einem korrekten Vergleich zu bisher verwendeten Konzepten wurde zusätzlich eine kon- 
ventionelle Geometrie untersucht. 

Die Ergebnisse der experimentellen Untersuchungen zeigen, dass die Verwendung von Drall- 
strömungen ein vielversprechendes Konzept zur Erhöhung des internen Wärmeübergangs ist, 
wenn ein erhöhter Eintrittstotaldruck der Kühlluft mit dem restlichen Maschinenkonzept zu 
vereinbaren ist. Auf der Druckseite konnten Verbesserungen der Totalkühleffektivität von bis 
zu 10 % erreicht werden, auf der Saugseite sogar von bis zu 20 %. Die Verbesserung auf der 
Saugseite ist zusätzlich stark abhängig von der Filmkühlung, die wiederum durch die komple- 
xe Innenströmung beeinflusst wird. Eine Trennung der Effekte ist bisher nicht möglich. Am 
thermisch hoch belasteten Staupunkt konnte keines der untersuchten Drallkühlkonzepte eine 
verbesserte Kühlwirkung erzielen. Durch erhöhte Wandstärken und eine reduzierte Filmkühlef- 
fektivität verbleiben trotz erhöhtem internen Wärmeübergang keine Netto-Vorteile. 

Die Verwendung von Drallkühlung kann demnach zu einer deutlichen Reduktion der Materi- 
altemperatur und damit zu einer erhöhten Lebensdauer von Turbinenschaufeln beitragen. Eine 
isolierte Anpassung von Teilen der internen Geometrien zur Implementierung eines Drallerzeu- 
gers sind jedoch nicht ausreichend. Aufgrund der von starker Querströmung geprägten Innen- 
strömung, deren Drallintensität und damit deren Wärmeübergang sich stark mit der Lauflänge 
ändert, müssen gänzlich neue Kühlungsgeometrien entworfen werden. Aus den präsentierten 
Ergebnissen können zwei grundlegende Empfehlungen abgeleitet werden: Die Innengeometrie 
und der Drallerzeuger müssen zuerst auf einen hohen internen Wärmeübergang ausgelegt wer- 
den. Auf Basis der resultierenden internen Strömung, speziell der Richtung und dem Betrag der 
Querströmung in den Tangentialebenen, muss dann eine abgestimmte Filmkühlung entwickelt 
werden. Hierbei ist die additive Fertigung als wichtige Technologie zu sehen, die eine freie Ge- 
staltung der Innengeometrie und damit die Abstimmung komplexer Kühlsysteme ermöglicht. 
Über die Notwendigkeit und die Möglichkeiten einer Erneuerung der Drallintensität müssen 
weitere Untersuchungen durchgeführt werden. 


Auch in Zukunft sollte das Potential des vorgestellten Prüfstandes und der vorgestellten Metho- 
den ausgeschöpft werden. Neue Konzepte können durch Verwendung der additiven Fertigungs- 
verfahren schnell realisiert und hoch genau evaluiert werden. Hier ist die Weiterverfolgung der 
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Drallkühlung nur eine Option. Eine Untersuchung verbesserter konventioneller Verfahren oder 
von Prallkühlverfahren, die zunehmend für den Einsatz in Rotorschaufeln erforscht werden, 
ist ebenso möglich. Einzig die Außengeometrie des Prüflings ist vorerst durch die Konturen 
des externen Kanals beschränkt. Soll diese auch verändert werden, müssen neue Seitenwände 
gefertigt werden. Eine neue Fertigung der Seitenwände ermöglicht weiterhin auch die Imple- 
mentierung weiterer optischer Zugänge. So ist gerade die Sicht auf die stromab liegende Druck- 
seite und druckseitige Hinterkante mit dem jetzigen Aufbau nicht möglich. Diese Bereiche sind 
aber bei modernen Schaufeln thermisch hoch belastet, die Kühlung der Hinterkante selbst kann 
als eigener Teilbereich der Fachliteratur gesehen werden. Untersuchungen bei der gegebenen 
aerothermischen Ähnlichkeit können wichtige Informationen liefern. Die vorgestellten Kali- 
brierverfahren können ohne Einschränkung verwendet werden. 

Trotz der möglichen Verwendung der vorgestellten Kalibrierverfahren muss die in-situ Kali- 
brierung weiter verbessert werden. Es wurde gezeigt, dass die eingebrachten Thermoelemente 
aufgrund der unterschiedlichen Biotzahlen von Substrat und Beschichtung nicht zwingend die 
Oberflächentemperatur abbilden. Eine Anpassung der Biotzahlen von Substrat und Beschich- 
tung auf dem Kalibrierkörper führt dagegen zu hohen Temperaturgradienten. Hier ist die Ent- 
wicklung von neuen Korrekturverfahren nötig. Alternativ zeigen Ianiro et al. (2010) sowie Hi- 
jazı et al. (2011) Verfahren, die durch Kombination von Messungen bei verschiedenen Wellen- 
längen die Anforderung an die Oberflächen reduzieren können. 

Bei komplexen geometrischen Aufbauten muss die Annahme einer konstanten Offsetstrahlung 
Eorr genauer untersucht werden. Zwar wird diese durch den Sichtwinkel lokal korrigiert, ihr 
Betrag bleibt aber konstant. Trotz technisch vollkommen diffuser Reflexion der Oberfläche ist 
davon auszugehen, dass sich geometrische Strukturen unterschiedlicher Körperstrahlung tief- 
passgefiltert auf der Messoberfläche abbilden. Dies ergibt sich zum Beispiel bei Fenstern (ge- 
ringe Eigenstrahlung), die nahe an den zu untersuchenden Oberflächen liegen, den Sichtbereich 
aber nicht vollständig überdecken. Hier sollte in Zukunft durch Verwendung von raytracing 
bestimmt werden, welche Anteile der Umgebung zur reflektierten Strahlung an jedem Objekt- 
punkt beitragen. Eine Korrektur kann ähnlich der Sichtwinkelkorrektur mittels eines weiteren 
skalaren, ortsabhängigem Faktors durchgeführt werden. 

Weiterhin müssen die zur Totalkühleffektivität beitragenden Effekte besser getrennt und ein- 
zeln quantifiziert werden. Als wichtigste Einflussgröße ist die adiabate Filmkühleffektivität zu 
nennen. Durch deren Bestimmung könnten die Effekte der Filmkühlung von denen der internen 
Kühlung getrennt werden. Da die Informationen an den Innenflächen der eingesetzten Geome- 
trien nicht bekannt sind, kann das Superpositionsprinzip nicht angewendet werden. Eine Be- 
stimmung dieser internen Größen ist nur unter größtem Aufwand möglich. Egger et al. (2013) 
zeigt ein transientes Verfahren, dass die Bestimmung der internen Größen durch eine Messung 
an den Außenflächen und Korrektur von lateraler Wärmeleitung im Material ermöglicht. Eine 
wichtige Bedingung hierbei ist jedoch eine niedrige Biotzahl, was bei Einhaltung der Maschi- 
nenähnlichkeit nicht gegeben ist. Eine weitere Korrektur könnte aufgrund der diffusiven Natur 
der Wärmeleitungsgleichung nur näherungsweise erfolgen. 

Zur Bestimmung der isolierten Filmkühlgrößen können aber andere Methoden verwendet wer- 
den. Stimpson et al. (2018) zeigen ein Verfahren, das eine Näherung der Größen durch Untersu- 
chung eines weiteren Prüflings ohne Filmkühlung ermöglicht. Eine weitere Möglichkeit bietet 
die Verwendung von druck- oder temperatursensitiven Farben. Aufgrund der Temperaturbe- 
schränkungen könnten diese Untersuchungen nur reynolds-, aber nicht machähnlich durchge- 
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fiihrt werden. Letztlich zeigen Die Methoden zur Transformation der Bilddaten in den Objekt- 
raum könnten bei den genannten Verfahren uneingeschränkt verwendet werden. 

Außerdem gilt es, die verwendeten Berechnungsverfahren (Schätzung der Kameraposition, In- 
terpolation der Daten) hinsichtlich der Geschwindigkeit zu optimieren. Durch konsequente Um- 
setzung in maschinennaher Programmiersprache ließe sich die Dauer der Auswertung reduzie- 
ren. Eine erste Betrachtung und grobe Korrektur der Daten im Objektraum könnte so schon 
während des Versuchs ermöglicht werden. 
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A Anhang 


A.1 Interpolation der Bilddaten in 3D 


Die Abbildung mittels der Kameramatrix ist, wie in Abschnitt 4.1 beschrieben, nicht invertier- 
bar. Zur Entzerrung und besseren Interpretation der Daten ist eine Darstellung derselben auf 
der Priiflingsgeometrie wiinschenswert. Dafiir werden die Objektpunkte auf die Bildebene pro- 
jiziert. Dies kann vorwärts mit der bestimmten Kameramatrix ausgeführt werden: 


Xoi = P Xoi: (A.1) 


Punkte, die außerhalb des Bildsensor liegen, werden maskiert. Die verbleibenden Punkte ha- 
ben als Dimension px in Bildkoordinaten. Dieser Raum ist noch kontinuierlich, während die 
Bilddaten nur an diskreten (also ganzzahligen) Koordinaten vorliegen. Die jeweiligen Daten an 
den kontinuierlichen Punkten werden mittels lokaler linearer Interpolation ermittelt. So kann je- 
dem Objektpunkt ein Bilddatenwert zugeordnet werden. Falls verschiedene Bilder der gleichen 
Kameraposition interpoliert werden, setzt dies voraus, dass die Kameraposition vollkommen 
unverändert bleibt. Da dies technisch nicht zu realisieren ist, werden die Bilder vor der Daten- 
interpolation an einem Referenzbild mittels dem ECC-Algorithmus (Evangelidis et al., 2008) 
ausgerichtet. So können kleine Änderungen der Kameraposition ausgeglichen werden, die auf- 
wändige Berechnung der Position muss pro Sichtbereich nur einmal durchgeführt werden. 

Die Tiefeninformation bei Verwendung einer einzelnen Kamera kann nicht direkt bestimmt wer- 
den und muss daher separat rekonstruiert werden. Das Verfahren wird im folgenden Abschnitt 
vorgestellt. 


A.2 Sichtfeldbasierte Maskierung 


Die ermittelten Kameramatrizen bzw. die daraus resultierenden linearen Abbildungen sind kon- 
tinuierliche Transformationen. Wird eine dreidimensionale Geometrie, also die Punkte des di- 
kretisierten Messobjektes, auf die Bildebene abgebildet, wie es fiir das in Abschnitt A.1 be- 
schriebene Verfahren notwendig ist, folgt aus der Abbildung keine Information tiber sichtbare 
Bereiche (Verlust der Tiefeninformation). Dies ist aber gerade bei den verwendeten Geometri- 
en mit starken Hinterschneidungen eine elementare Information bei der Zuordnung der Daten. 
Zur vollen Rekonstruktion der Daten muss bekannt sein, welche Flächenelemente aus Sicht der 
jeweiligen Kamera zuvorderst liegen. 

An jedem Punkt der diskretisierten Geometrie, der dem Detektor zugeordnet werden kann, wird 
dafür ein Strahl gebildet. Dieser Strahl ist durch den jeweiligen Punkt X; und das optische Zen- 
trum C definiert. Für diesen Strahl werden die Schnittpunkte mit dem diskretisierten Körper 
berechnet. Dies kann unter Verwendung von Plückerkoordinaten durchgeführt werden - das 
verwendete Verfahren wird von Teller (1992) vorgestellt und von Amanatides et al. (1995) wei- 
ter optimiert. 

Um eine allgemeine Anwendung zu ermöglichen, müssen konvex-konkave Oberflächen mit 
Hinterschneidungen berücksichtigt werden können. Dazu wird an jedem Schnittpunkt zusätz- 
lich die Richtung des Strahles im Vergleich zum lokalen Oberflächenvektor bestimmt. Für jeden 
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Strahl ergibt sich dann die Anzahl der Schnitte sowie die Information, ob der Strahl in den Kör- 
per eindringt (und damit valide ist), oder aus dem Körper austritt (nicht valide). Aus logischem 
Überlegen ergibt sich die Sichtbarkeitsinformation für jeden Strahl. Nur in Ausnahmefällen! 
muss auf eine aufwändige Berechnung der Schnittpunktkoordinaten zurückgegriffen werden. 
Diese wird, wenn nötig, mithilfe des Möller-Trumbore-Algorithmus (Möller et al., 1997) durch- 
geführt. So ist das hier vorgestellte Verfahren trotz der fein aufgelösten Oberflächen mit vielen 
Elementen (ng, > @(10°)) schnell. 


A.3 3D Panorama 


Das Ziel der experimentellen Untersuchungen ist eine optimale optische Abdeckung des Priif- 
lings. Da die einzelnen Bildereiche während der Auslegung nur abgeschätzt werden können, 
ergeben sich Bereiche, in denen sich mehrere Bilder überlappen. Aufgrund der angewandten 
Methodik zur Kalibrierung und Transformation kann ein Zusammensetzen der Daten von meh- 
reren Aufnahmen von verschiedenen Kameras erst im Objektraum durchgeführt werden. Be- 
kannte Methoden setzen aber die Zweidimensionalität der Daten voraus. Zwar bleibt die Di- 
mensionalität bei der Rücktransformation erhalten, zur Reduktion der implizit notwendigen 
dritten Dimension muss aber die u-v-Parametriesierung der Oberfläche global definiert (und be- 
kannt) sein? . Dies ist eine Einschränkung an die bisher erreichte Freiheit der zu analysierenden 
Oberfläche. So sind mit modernen Methoden erstellte Freiformflächen meist nur stückweise 
(oder sogar infinitisemal) kontinuierlich parametrisiert. Eine dokumentierte Methodik zum Zu- 
sammenfügen der Daten in drei Dimensionen ist dem Autor zu diesem Zeitpunkt nicht bekannt 
und wird daher im Folgenden präsentiert. Das Verfahren beschränkt sich auf die Überlappung 
von genau zwei Datensätzen. Sollten sich in seltenen Fällen mehr als zwei Datensätze an einem 
Ort überlappen, kann eine schrittweise Überblendung durchgeführt werden. 

Bevor die einzelnen Datensätze zusammengefügt werden können, muss sichergestellt sein, dass 
die Daten an stationären Zeitpunkten bei konsistenten Bedingungen aufgenommen werden. 
Aufgrund der teilweise langen Zeitspannen (mehrere Minuten) für eine komplette Untersu- 
chung an einem Betriebspunkt und verbleibenden Restschwankungen der Randbedingungen ist 
dies nur für die entdimensionierten Totalkühleffektivitäten TCE gegeben. Daher können auch 
nur diese finalen Daten sinnvoll kombiniert werden. 

Der Einfluss des Sichtwinkels wurde für diese Daten schon korrigiert. Unter Annahme einer 
fehlerfreien Korrektur verbleibt als letztes Qualitätskriterium der Daten die Distanz von der 
Bildmitte im Bildraum. Diese beeinflusst die Qualität des Bildes und damit der Daten, da die 
Fehler optischer Systeme im Allgemeinen mit steigendem Radius zunehmen. Die auf diesem 
Qualitätskriterium basierende Methodik wird im Folgenden vorgestellt: 


' Dies ergibt sich bei ungerader Anzahl von Schnittpunkten, deren gewichtete Ein-Aus-Summe nicht 0 ergibt. 
Diese treten nur bei beschriebenen Körpern mit Hinterschneidungen auf. 

2 Solche Parametrisierungen finden bei der Entzerrung mittels Zylinder- oder Kugelpanoramamethoden Anwen- 
dung. Ist diese Parametrisierung global definiert, können aus der Fotografie bekannte Methoden zur Panorama- 
Erzeugung im zweidimensionalen u-v-Raum angewandt werden. 
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1. Der geometrische Uberlappungsbereich wird im Objektraum bestimmt. An N Punkten der 
diskretisierten Oberfläche (Menge D) liegen jeweils zwei Werte der Totalkiihleffektivitat 
vor, die i. A. unterschiedlich sind. 


2. Die Bildmitten (Berechnung nach Abschnitt 4.3) der jeweiligen Kameras werden in den 
Objektraum transformiert. Die Richtung des Verbindungsvektors dieser Punkte wird be- 
stimmt. Für jeden der zwei Datensätze wird der Abstand zu den jeweiligen Bildmitten in 
Richtung dieser Verbindungsrichtung bestimmt. 


3. Im Allgemeinen ist keine Kante des Überlappungsbereiches D senkrecht auf den Ver- 
bindungsvektor der Bildmitten. Dies kann zu Sprüngen in den Gewichtungswerten an 
den Rändern von D führen. Daher wird das volumenminimierte, gerichtete kartesische 
Hüllvolumen (oriented bounding box, OBB) der Menge D bestimmt, das eine genaue 
Richtungsbestimmung der Ränder erlaubt. 


4. Die Winkel zwischen den Achsen der OBB und dem Verbindungsvektor wird bestimmt. 
Die Hauptachse mit geringstem Winkel wird als neue Überblendungsrichtung definiert. 
Die vorher bestimmten Abstände werden auf diese Richtung projiziert. Nach Normierung 
ergibt sich für jeden Datensatz an jedem Punkt ein Gewichtungsfaktor wi(X ). Es gilt 
YL wi(X = const.) = 1. 


5. Die Daten werden gewichtet aufsummiert. 
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A.4 Kühleffektivitäten Nennbetriebspunkt, absolut 


0,55 


0,45 


0,40 


P4 SP S9 
Position, param. Lauflänge [-] 


Abbildung A.1: Verlauf der gemittelten Totalkühleffektivität, Vergleich der Geometrien bei 
Nennbetriebspunkt, absolute Darstellung 
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A.5 Kihleffektivitaten Nennbetriebspunkt, relativ 


TCE rel. [-] 


Position, param. Lauflange [-] 


Abbildung A.2: Verlauf der gemittelten Totalkühleffektivität, Vergleich der Geometrien bei 
Nennbetriebspunkt, relative Darstellung 
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A.6 Ablaufplan des genetischen Algorithmus 


Auswahl von n Auswahl von n 
Objektpunkten X Bildpunkten X 


Population erzeugen 
X; =X 
Xj = xX + Apx 
u = { (Xixi) Vi € u} 


Optimierer starten 


u Individuen 


(u +A) Individuen | Population evaluieren 
Eignung bestimmen 


A neue Ind. 

Mutation und Statistik 
Kreuzung / auswerten / 

Pm, Pmi, Pex 


Finale Kamera- 
matrix bestimmen 
Popt = f(Xi xi) 


Auswahl der u 
bestgeeigneten Ind. 
NSGA-II-Algorithmus 


Abbildung A.3: Ablaufplan des genetischen Algorithmus 
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werden, vor allem die Vorhersage der Materialtemperatur der Schaufeln 
muss hoch präzise erfolgen. 


Die vorliegende Arbeit präsentiert neue Verfahren zur Kühlung von 
Rotorschaufeln. Es wird speziell auf die Drallkühlung eingegangen, 
die sich in grundlegenden Untersuchungen und in Verbindung mit 
den sich rapide verbessernden Möglichkeiten der additiven Fertigung 
als vielversprechendes Verfahren herausgestellt hat. Zur Beurteilung 
der neuen Konzepte werden Untersuchungen durchgeführt, bei denen 
die wichtigen strömungsmechanischen Ahnlichkeitskennzahlen der 
Maschine eingehalten werden. So können die Ergebnisse unmittelbar 
auf den Einsatz im Triebwerk übertragen werden. Die Untersu- 
chungen bedeuten eine Herausforderung für die als Messtechnik 
eingesetzte Infrarotthermographie. Neue Kalibrierverfahren werden 
entwickelt und überprüft. Durch Anwendung dieser Verfahren kann die 
Oberflachentemperatur der Prüflinge hoch genau bestimmt werden. 
Anhand der mit den neuen Verfahren bestimmten Verteilung der Ober- 
flachentemperatur und weiteren aerodynamischen Messungen werden 
abschließend Empfehlungen formuliert, wie ein auf Drallkühlung 
basierendes Kühlsystem für Turbinenrotorschaufeln leistungsfähiger 
werden kann und so die geforderte Lebensdauer dieser Bauteile 
sicherstellen kann. 
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